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In  this  dissertation  we  develop  an  analytic  model  for  minority- 
carrier  transport  in  polysilicon  devices,  and  provide  experimental 
corroboration  for  the  model.  The  model  is  used  to  facilitate  the 
development  of  experimental  techniques,  compatible  with  conventional 
device  processing,  to  control  the  effects  of  grain  boundaries. 
Techniques  are  investigated  to  reduce  bulk  recombination  current  by 
gettering  intragrain   impurities  in   polysilicon. 

Key  assumptions  are  made,  with  justification,  to  simplify  the 
three-dimensional,  nonlinear  boundary-value  problem  that  defines 
minority-carrier  transport,  including  recombination,  in  polysilicon 
devices.  These  assumptions  enable  the  separation  of  the  grain-boundary 
recombination  analysis,  which  is  based  on  quasi-equil  ibri  um,  from  the 
intragrain  transport  analysis,  which  is  done  by  partitioning  the  grain 
into  subregions  in  which  the  minority-carrier  flow  is  predominantly  one- 
dimensional.  The  analyses  are  coupled  through  the  effective  minority- 
carrier  recombination  velocity  at  the  grain  boundary,  which  generally  is 


dependent  on  the  minority-carrier  density  in  the  adjacent  quasi-neutral 
grain. 

The  dependence  of  the  recombination  velocity  on  the  carrier  density 
(excitation)  at  the  grain  boundary  is  experimentally  demonstrated  using 
electron-beam-induced-current      (EBIC)      measurements.  To      facilitate 

quantitative  interpretation  of  the  EBIC  measurements,  we  solve  the 
underlying  carrier  transport  problem  subject  to  the  nonlinear  boundary 
condition  at  the  grain  boundary  using  computer-aided  numerical 
analysis.  By  comparing  the  theoretical  and  experimental  EBIC  responses, 
we  predict  values  for  the  surface- state  density  at  typical  grain 
boundaries  and  the  minority-carrier  diffusion  length  in  the  grains. 

Experimental  results  are  presented  that  imply  potential 
improvements  afforded  by  aluminum  diffusion  in  both  bulk  and  thin- film 
polysilicon  solar  cells.  For  bulk  cells,  a  high- temperature  aluminum 
diffusion  (alloying)  is  shown  to  increase  the  minority-carrier  diffusion 
length  by  gettering  intragrain  impurities.  For  thin- film  cells,  a  low- 
temperature  aluminin  diffusion  is  shown  to  substantially  passivate  grain 
boundaries  and  hence  decrease  the  recombination  velocity.  The  decrease 
is  evaluated  using  EBIC  measurements,  the  interpretation  of  which  is 
aided  by  the  numerical   carrier  transport  model    developed. 
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CHAPTER   1 
INTRODUCTION 

Polycrystalline  silicon  (polysil  icon)  is  being  widely  used  in  the 
semiconductor  industry  to  fabricate  unipolar  as  well  as  bipolar  devices 
[1].  The  emergence  of  polysil  icon  devices  has  motivated  studies  of 
majority-carrier  transport  through  the  grain  boundaries  [2-6],  which 
defines  the  resistivity  of  the  polysil  icon,  as  well  as  studies  of 
minority-carrier  recombination  at  the  grain  boundaries  [7-11],  which 
defines  the  bipolar  characteristics  of  the   polysilicon. 

Compared  to  silicon  (majority-  and  minority-carrier)  devices, 
polysilicon  devices  perform  poorly.  The  grain  boundaries,  in  general, 
are  responsible  for  the  poor  performance  of  the  polysilicon  devices. 
For  example,  in  polysilicon  MOSFETs  the  presence  of  grain  boundaries 
causes  a  turn-on  characteristic  that  is  beyond  the  strong-inversion 
threshold  [5],  while  in  thin- film  polysilicon  solar  cells  [12]  the 
conversion  efficiency  is  lowered  due  to  recombination  at  the  grain 
boundaries.  The  understanding  of  how  grain  boundaries  influence  the 
transport  of  carriers  in  polysilicon  is  hence  of  engineering 
signi  ficance. 

To  enable  the  design  and  development  of  polysilicon  solar  cells 
that  are  cost-competitive  with  silicon  cells,  we  must  understand  how  the 
grain  boundaries  influence  the  minority-carrier  transport  and  then 
devise  fabrication  techniques  to  minimize  this  influence.  Although 
several  approximate  models  [7-10,  13,  14]  for  the  minority-carrier 
transport  in   polysilicon   have  been  developed,   they  are   inadequate  in  the 
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general  case.  The  reason  for  this  is  that  the  models  [7-10,  13]  are 
based  on  restricting  assumptions  that  limit  their  ranges  of  validity; 
furthermore,  the  solution  of  the  minority-carrier  transport  equation  is 
obtained  by  truncating  an  infinite  series  [8,  14],  and  hence  the 
accuracy  of  the  models  is  questionable.  Many  processing  (passivating) 
techniques  [15-19]  have  also  been  developed  to  control  the  grain- 
boundary  recombination.  However  none  of  these  techniques  is  easily 
adapted  to   solar  cell    processing  and/or  is  totally  effective  [15-19]. 

The  purpose  of  this  research  is  to  develop  a  general  minority- 
carrier  transport  model  for  polysilicon  that  is  based  on  a  good  physical 
understanding  of  grain  boundaries,  and  to  use  it  to  develop  experimental 
techniques,  compatible  with  conventional  device  processing,  to  control 
the  effects  of  grain  boundaries;  for  example,  to  passivate  grain 
boundaries  and  thereby  improve  the  efficiency  of  polysilicon  solar 
cells.  Another  limitation  on  the  efficiency  of  polysilicon  solar  cells 
is  due  to  the  base  minority-carrier  lifetime,  which  is  typically  shorter 
than  that  in  silicon  cells  due  to  high  intragrain  defect  densities. 
Hence  another  purpose  of  this  research  is  to  develop  experimental 
techniques  to  getter  intragrain  impurities.  The  primary  tasks  of  the 
research  are 

(1)  to  develop  an  analytic  model  for  the  minority-carrier  transport 
in   polysilicon; 

(2)  to  provide  experimental  corroboration  for  the  transport  model 
and  to  determine  values  for  the  pertinent  grain-boundary 
parameters; 

(3)  to  develop  experimental  techniques  to  improve  the  performance 
of  bulk  and  thin- film  polysilicon  solar  cells  and  other  bipolar 
devices. 


Polycrystall  ine  silicon,  as  the  nomenclature  implies,  has  more  than 
one  crystalline  orientation  (grain),  each  orientation  is  separated  by  a 
plane  of  dislocations  (grain  boundary).  At  the  grain  boundaries  are 
defects  [1],  e.g.,  dangling  bonds,  which  produce  localized  electron 
states  (traps)  throughout  the  energy  gap.  The  net  charge  on  these 
states,  which  results  from  majority-carrier  trapping,  gives  rise  to  a 
potential  barrier  which  influences  the  conduction  properties  of  the 
polysilicon  [6].  In  minority-carrier  devices,  such  as  solar  cells,  the 
grain-boundary  states  serve  as  recombination  centers  fbr  the  minority 
carriers  and  hence  constitute  a  source  of  recombination  current.  Thus 
the  development  of  the  minority-carrier  transport  model  fbr  polysilicon 
is  complicated  by  the  presence  of  randomly  oriented  surfaces  (grain 
boundaries)   at  which  significant  recombination  can  occur. 

Impedance  measurements  [20],  current-voltage  measurements  [15,  16], 
and  electron-beam-induced-current  (EBIC)  measurements  [21,  22]  have  been 
used  to  characterize  the  barrier,  or  the  defect  density  at  the  grain 
boundary,  either  of  which  is  a  measure  of  the  grain-boundary 
recombination.  We  use  the  EBIC  measurements  to  infer  the  grain-boundary 
parameters  since  it  facilitates  isolation  of  the  grain  boundary  from  the 
bulk  grain.  Such  an  isolation  is  beneficial  in  order  to  study  the 
grain-boundary  characteristics,  e.g.,  the  influence  of  the  excitation  on 
grain-boundary  recombination.  Furthermore  the  EBIC  measurements  do  not 
require  tedious  sample  preparation  as  do  current-voltage  measurements 
[15,  16]  to  obtain  the  grain-boundary  parameters.  An  additional 
advantage  is  that  we  can  infer  the  lifetime  in  the  grain  from  EBIC 
measurements. 


The  presence  of  defect  states  results  in  current  loss  at  the  grain 
boundary  and  hence  deteriorates  the  performance  of  polysilicon 
devices.  To  control  the  detrimental  influence  of  the  grain  boundaries, 
it  is  necessary  to  reduce  the  defect  states  at  the  grain  boundary.  This 
reduction  can  be  achieved  by  selective  (preferential)  diffusion  of 
certain  impurities  [16-19]  down  the  grain  boundaries  or  by  fabrication 
of  special  structures  to  eliminate  them  from  the  active  portions  of  the 
device  [15].  Furthermore,  polysilicon  contains  high  intragrain  defect 
densities  which  cause  short  minority-carrier  lifetimes.  Hence  it  is 
necessary  to  consider  processes  that  can  getter  the  intragrain 
impurities  and  hence  increase  the  minority-carrier  lifetime. 

In  Chapter  2,  we  develop  an  analytic  model  for  the  minority-carrier 
transport  in  polysilicon.  The  general  transport  problem,  which  is 
three-dimensional  and  nonlinear  [23],  is  simplified  by  making  key 
assumptions  with  physical  justification.  The  main  feature  of  the  model 
is  the  separation  of  the  grain-boundary  analysis,  which  yields  an 
expression  for  the  effective  recombination  velocity,  from  the  grain 
analysis,  which  yields  an  expression  for  the  base  recombination 
current.  The  analyses  are  coupled  through  the  effective  recombination 
velocity  at  the  grain  boundary. 

Our  results  indicate  that  the  effective  recombination  velocity  is, 
in  general,  dependent  on  the  excitation  (minority-carrier  density)  at 
the  grain  boundary,  and  in  fact  decreases  with  increasing  excitation. 
This  result  renders  most  transport  analyses  [13,  14],  which  are  based  on 
the  common  assumption  that  the  effective  recombination  velocity  is  a 
constant,  not  valid  in  general.  Our  results  also  indicate  that  at  low 
and    intermediate    forward   voltages,    the   base    recombination   current   in    a 


polysilicon  diode  comprises  a  component  independent  of  the  grain- 
boundary  parameters  and  a  component  strongly  dependent  on  the  grain- 
boundary  parameters.  Contrarily,  at  high  forward  voltages,  the  grain- 
boundaries  have  negligible  effect  on  the  base  recombination  current,  and 
the  polysilicon  virtually  behaves  like   single-crystal    silicon. 

In  Chapter  3,  we  characterize  the  nonlinear  effective  recombination 
velocity  at  the  grain  boundary  using  EBIC  measurements  [21,  22].  The 
electron  beam  is  incident  from  the  top  surface  of  a  polysilicon  cell  and 
generates  electron-hole  pairs,  predominantly  in  the  base  region  of  the 
cell.  The  generation  (short-circuit)  current  is  measured  directly  using 
an  ammeter.  The  EBIC  response  of  a  grain  in  the  vicinity  of  a  grain 
boundary  is  obtained  by  measuring  the  generation  current  as  a  function 
of  the  distance  of  the  center  of  the  beam   from  the  grain  boundary. 

To  facilitate  quantitative  interpretation  of  the  EBIC  response,  we 
develop  a  computer-aided  nunerical  analysis  of  the  underlying  minority- 
carrier  transport  problem.  Several  numerical  [24]  and  analytical 
[25,  26]  solutions  have  been  derived  for  the  transport  problem.  However 
their  ranges  of  validity  are  limited  owing  to  simplifying  assumptions  of 
a  point  [26]  or  a  spherical  [24,  25]  generation  source,  and  a  constant 
grain-boundary  recombination  velocity  [24-26].  Our  analysis  removes 
these  restricting  assumptions,  and  hence  is  more  general.  The  computer- 
aided  numerical  analysis,  in  this  respect,  is  advantageous.  The 
experimentally  obtained  EBIC  response  shows  excellent  agreement  with  our 
numerical  analysis  and  demonstrates  clearly  the  dependence  of  the 
effective  recombination  velocity  on  the  excitation  at  the  grain 
boundary,  which  provides  corroboration  for  the  transport  model  in 
Chapter     2.  By     comparing     the     theoretical      and     experimental      EBIC 


responses,  we  predict  values  for  the  defect  density  at  typical  grain 
boundaries  and  the  minority-carrier  diffusion  length  in   the  grains. 

Having  demonstrated  the  detrimental  influences  of  the  grain- 
boundary  recombination  on  the  minority-carrier  transport,  we  next  focus 
our  attention  on  experimental  methods  to  improve  the  performance  of 
polysilicon  solar  cells.  In  Chapter  4,  we  describe  an  experimental 
technique,  using  aluminum  diffusion,  to  passivate  the  grain  boundaries, 
i.e.,  to  reduce  the  defect  density  at  the  grain  boundaries.  Aluminum 
was  chosen  owing  to  its  compatibility  with  solar  cell  processing,  and  to 
its  large  diffusion  coefficient  down  the  grain  boundaries.  It  can  also 
concomitantly  getter  the  intragrain  impurities,  i.e.,  reduce  the  defect 
density  in  the  grains,  which  results  in  an  increase  in  the  minority- 
carrier  lifetime.  The  minority-carrier  lifetime  in  polysilicon  is 
typically  lower  than  that  in  silicon  due  to  high  intragrain  defect 
densities,  and  hence  it  is  useful  to  consider  processes  that  improve 
lifetime. 

For  the  development  of  bulk  polysilicon  solar  cells,  in  which  the 
intragrain  recombination  is  dominant,  v/e  propose  a  high- temperature 
(above  eutectic,  577°C)  al  uminum-di  f fusion  process  from  the  back 
surface.  For  the  development  of  thin- film  cells,  in  which  the  grain- 
boundary  recombination  is  dominant,  we  propose  a  low- temperature  (below 
eutectic)  aluminum-diffusion  process  from  the  front  surface  of  the 
polysilicon.  We      use      forward      current-voltage      and      reverse-bias 

capacitance  measurements  to  demonstrate  the  intragrain  gettering,  and 
EBIC  measurements  to  show  the  grain-boundary  passivation.  Our  results 
indicate  that  the  high- temperature  aluminum-diffusion  process 
effectively    getters    intragrain    impurities,    whereas   the    low- temperature 


aluminum-di  f  fusion  process  produces  significant  grain-boundary 
passivation.  We  also  find  that  the  grain  boundaries  emit  impurities 
during  the  high- temperature  process  and  that  these  emitted  impurities 
are  effectively  gettered  by  the  aluminum. 

Although  we  attribute  the  benefits  produced  by  the  high- temperature 
aluminum-di  f fusion  process  in  Chapter  4  to  intragrain  gettering,  the 
observed  improvement  in  polysilicon  device  performance  is  also 
commensurate  with  the  fbmiation  of  a  back-surface  field  (BSF)  [27].  The 
BSF,  which  is  nonuniform  and  nonplanar,  is  produced  as  a  result  of  the 
preferential  dissolution  of  silicon  during  the  high- temperature 
aluminum-diffusion  process  [28].  In  Chapter  5,  we  examine  the 
mechanisms  of  gettering  and  BSF  formation  in  an  attempt  to  isolate  the 
actual  mechanism  that  is  instrumental  in  improving  the  polysilicon 
device  performance.  The  studies  are  based  on  forward  current-voltage 
measurements  of  al  umin  urn- di  f  fused  Wacker  (cast)  polysilicon  cells,  made 
before  and  after  lapping  off  the  back  p'^p  junction.  Our  results 
indicate  that  the  observed  improvement  in  the  polysilicon  device 
performance  is  due  to  effective  gettering  of  intragrain  impurities  that 
occurs  during  the  high- temperature  aluminum-diffusion  process,  and  that 
the  nonuniform  BSF  is  ineffective.  No  improvement  in  the  perfbrmance  of 
semiconductor-grade  and  solar-grade  [29]  silicon  cells  is  observed,  and 
hence  we  conclude  that  the  benefits  resulting  from  the  aluminum- 
diffusion  process  are  peculiar  to  the  (cast)  polysilicon.  The  benefits 
result  possibly  because  of  unique  impurities  prevalent  in  the 
polysilicon  that  are  readily  gettered  by  the  aluminum.  We  surmise  that 
aluminim  getters  oxygen  and/or  oxygen  related  defects  which  are  present 
in  higher  densities  in  Wacker  (cast)  polysilicon  than  in  semiconductor 
grade  and   solar-grade   silicon. 


We  sunman'ze  our  studies  and  review  the  main  conclusions  and 
accomplishments  of  this  dissertation  in  Chapter  6.  We  also  discuss  the 
scope  and  limitations  of  this  work  and  provide  suggestions  for  future 
research. 


CHAPTER   2 

ANALYSIS   OF  MINORITY-CARRIER  TRANSPORT    IN 

POLYSILICON   DEVICES 


2.1     Introduction 


The  performance  of  bipolar  pol ycry stall ine  silicon  devices  can  be 
limited  because  of  minority-carrier  recombination  losses  at  the  grain 
boundaries  [7-11].  For  example,  the  efficiency  of  thin- film  polysilicon 
solar  cells  [12],  which  can  potentially  be  used  to  produce  cost- 
competitive  terrestrial  photovoltaic  energy-conversion  systems,  is  lower 
than  that  of  silicon  cells  because  of  grain-boundary  recombination. 
Hence  a  complete  theoretical  understanding  of  the  influence  of  grain 
boundaries  on  minority-carrier  transport  is  essential  in  order  to 
optimize  the  design  and  performance  of  polysilicon  solar  cells  and  other 
bipolar  polysilicon  devices. 

In  the  general  case  this  influence  is  described  by  a  complicated, 
three-dimensional  boundary-value  problem  having  nonl  inear  boundary 
conditions  [23].  For  example,  in  the  fbrward-biased  polysilicon  n'^'p 
junction  with  columnar  grains  illustrated  in  Fig.  2.1,  the  electron 
current  injected  into  the  quasi-neutral  fvtype  base  is  governed,  for 
low-injection  conditions,  by  the  three-dimensional  electron  continuity 
equation  subject  to  nonlinear  boundary  conditions  at  the  grain 
boundaries  [7-11].  The  three-dimensionality  results  basically  because 
the  carriers  are  injected  at  a  surface  (i.e.,  the  pn  junction)  that  is 
not     parallel    to  the  grain-boundary     surface.       The  boundary  conditions. 
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which  characterize  carrier  recombination  through  energy-gap  states  at 
the  grain  boundaries,  can  be  expressed  by  an  effective  recombination 
velocity  that,  in  general,  is  dependent  on  the  intragrain  electron 
density  and  on  position.  Approximate  solutions  [7-11,  13,  14,  30]  for 
this  and  related  problems  have  been  derived,  but  their  ranges  of 
validity   have  not  been  checked. 

To  help  advance  our  understanding  of  minority-carrier  transport  in 
poly  silicon,  and  to  aid  the  engineering  design  of  optimal  poly  silicon 
bipolar  devices,  including  the  solar  cell,  we  describe  in  this  chapter 
the  development  of  an  analytic  model  for  the  minority  electron  current 
injected  into  the  base  of  the  n"''p  junction  in  Fig.  2.1.  The  model, 
which  is  subject  to  the  uncertainties  and  inhomogeneities  in  the 
morphology  and  in  the  electrical  properties  of  polysilicon,  is  not 
exact,  but  is,  we  believe,  sufficiently  representative  of  the  physics  to 
ensure  its  utility.  Although  we  emphasize  polysilicon,  the  analytical 
methods  used  in  the  model  development  are  applicable  to  other 
polycrystall  ine   semiconductors. 

o 

The    electron    current  Ig  injected    into    the    base    of  each    grain    in 

n 

Fig.    2.1    derives    from    recombination    at    the    grain    boundaries  (Ipn)  and 

within   the  quasi-neutral    p-region  (Iqu)-     A  key    feature  of  the  model    is 

the     separation     of    the    grain-boundary    recombination     analysis,     which 

defines    the    effective    recombination    velocity  S  ,    ^«  ,    from    the    quasi- 

n  ( e  t  f ) 

neutral  grain  transport  analysis,  which,  with  S  ,  ..<  ,  defines  iL  The 
model  differs  from  previous  ones  [7-11,  13,  14,  30]  in  the  following 
ways:  (a)  it  is  based  on  a  computer-aided  numerical  solution  of 
Poisson's  equation  within  the  grain-boundary  space-charge  region,  which 
provides    the    physical     insight    necessary    to    analytically    characterize 
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PR 
^n(eff)''      ^'^^     ^^    unifies    previous    analyses    based    on    the    restricting 

assumptions    of   carrier    depletion    [7,    9,    10]    and    of  quasi-equil  ibrium 

[7,  8]     in     the    grain-boundary     space-charge    region,    as    well     as    those 

analyses    based     on     the     common     assumption     that      s'^?    ...    =  «  [13],   by 

incorporating    into    the    quasi-equil  ibrium    analysis    the    physical     upper 

PR  c 

limit  on  S^^^^^^,  i.e.,  the  kinetic-limit  velocity  S^^(=  5  x  lo  cm/ sec 
at  300°  K)  imposed  by  the  random  thermal  motion  of  the  conduction-band 
electrons  [31];  (c)  it  avoids  possible  errors  involved  in  truncating 
infinite-series  solutions  ibr  the  three-dimensional  electron  continuity 
equation  in  the  quasi-neutral  grain  [8,  14,  23,  30]  by  reducing  the 
problem  to  coupled  one-dimensional  ordinary  differential  equations  that 
are  ibmiulated  based  on  "gradual-case"  approximations  [32];  ( d)  it  is 
the  first  analytic  model  for  minority-carrier  transport  in  polysilicon 
derived  from  the  physical  insight  provided  by  computer-aided  numerical 
sol  utions. 

The  model  development  begins  with  the  computer  solution  of 
Poisson's  equation  in  the  grain-boundary  space-charge  region  adjoining 
quasi-neutral  p-type  grains,  assumed  to  be  in  lovv  injection.  The 
numerical  analysis  is  based  on  the  quasi-equil  ibri  un  assumption  [7,  8,] 
but  the  inherent  limitations  are  effectively  removed  by  accounting  for 
the  kinetic-limit  electron  velocity.  The  calculation  of  the  grain- 
boundary  potential  barrier  height  is  based  on  the  assumption  of  an 
"effective"  monoenergetic  density  of  donor-type  traps  located  near  the 
middle  of  the  energy  gap.  This  assumption  is  commensurate  with  the 
experimental  observation  [20]  that  the  "neutral  Fermi  level"  of 
polysilicon  grain  boundaries  is  near  midgap.  Once  the  barrier  height  is 
determined,     the    composite    recombination    rate    of    hole-electron     pairs 
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through  both  donor- type  and  acceptor- type  traps  is  calculated  accounting 
for  the  commonly  large  difference  between  the  carrier  capture  parameters 
for  neutral    and  ionized  traps  [9].     This  recombination  rate  is  then  used 

PR 

to     derive  S„,    ^.x,  which     is     a     function     of    the     electron     density. 
n(eff)  "^ 

PR 

Empirical  dependencies,  implied  by  the  nimerical  analysis,  of  S  ,^^s  on 
the  grain-boundary  surface-state  densities,  on  the  grain  doping  density, 
and  on  the  "excitation"  of  the  grain  boundary,  i.e.,  the  separation  of 
the    electron    and    hole    quasi-Fermi    levels,    are    used    to     formulate    an 

PR 

expression  for  S  ,  -^^ .  The  expression  is  supported  by  qualitative 
discussion  of  the  underlying  physics,  including  mention  of  the  effects 
of  assuming  different  trap  energy  levels. 

To  complete  the  model,  the  electron  continuity  equation  in  the 
quasi-neutral  p-type  base  of  a  representative  grain,  assumed  to  be  a 
right-circular  cylinder,   is  analytically   solved   subject  to  the  nonlinear 

PR 

boundary  condition  [S  ,  ^^xl  at  the  grain  boundaries,  which  varies  with 
excitation  along  the  boundaries,  and  to  the  excess  electron  density  at 
the  edge  of  the  junction  space-charge  region  resulting  from  the  forward 
bias     V.  This     complicated     three-dimensional      problem     is     solved 

analytically  for  relatively  large  grains  by  partitioning  the  grain  into 
subregions  in  which  the  electrons  are  assumed  to  flow  predominantly  in 
one  dimension.  The  one-dimensional  solutions  are  then  coupled  to  give 
the  composite  solution,  which  is  approximate  but  nevertheless 
illustrates  well  the  effects  of  the  nonl  inearity.  The  result  is  an 
analytic  expression  for  Ip  in  terms  of  V  and  the  parameters  defining  the 
morphology  and  electrical  properties  of  the  grain  and  the  encompassing 
grain  boundaries. 
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2.2     Analysis 

The  assumed  representative  grain,  a  right  circular  cylinder,  of  the 
polysilicon  n"^p  junction  in  Fig.  2.1  is  illustrated  in  Fig.  2.2(a).  The 
grain  and  grain  boundary  are  assumed  to  be  isotropic.  Hence  rotational 
symmetry  exists  about  the  central  axis,  and  the  carrier  transport 
problem  can  be  reduced  to  two  dimensions  as  indicated  in  Fig.  2.2(b). 
Our  objective  is  to  analytically  characterize  the  injected  electron 
current  I?  =  I?„  +  I?.,     that   results    from    the  junction    forward   bias  V. 

b  bD  I^N 

(Note  that  we  are  accounting  for  recombination  at  the  portion  of  the 
grain  boundary  adjacent  to  a  quasi-neutral  region.  Recombination  at  the 
grain  boundary  within  the  junction  space-charge  region,  which  can  be 
significant  [33],  is  discussed  in  Section  2.3.)  Rigorously  this  would 
require  the  solution  of  the  complicated  two-dimensional  electron 
continuity  equation  in  the  quasi-neutral  p-region  of  the  grain  subject 
to    appropriate    boundary    conditions    at    the    grain-boundary    and   junction 

space-charge  regions  [23].     The  derived  electron   density  N(r,z)   must  be 

B 
integrated    over    the    volume    of  the    p-region    to    determine  Iqj.,  and    its 

gradient    integrated    over    the    encompassing    grain-boundary    surfaces    to 

determine  Ig„.      Such   a    rigorous    solution    would,    even    with    simplifying 

assunptions,  require  a  computer-aided  nunerical    analysis  [13]. 

To      facilitate     the     engineering     design     of     polysilicon     bipolar 

devices,    and    to    advance    our    understanding    of    the    pertinent    physics 

underlying     their     performance,     an     approximate     analytic     solution     is 

useful.       To    derive    such    a    solution,    we    separate    the    analysis   of  the 

grain-boundary  recombination    from  that  of  the  intragrain  transport,   even 

though     the      mechanisms    are     interrelated.         This      interrelationship 
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REGION 


QUASI-NEUTRAL 
p-TYPE  BASE  OF  GRAIN 


Fig.  2.2  Representative  grain  (a)  of  the  forward-biased  n  p  junction 


and  its  cross-section  (b) 
,GB 


We  assume  r^>>W, 


,GB 


hence  the 


G  "SCR^ 
z-dependence  of  W^^p  deriving  from  the  z-dependence  of  the 

grain-boundary  excitation  is  inconsequential  and  has  not  been 
shown. 
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materializes   when    we    combine    the    results   of  the    analyses    to    determine 

D 

N(r,z)    and  to  characterize  Ip(V). 


2.2.1     Grain  Boundary 

We  first  consider  the  grain-boundary  space-charge  region,  which  is 
identified      in      Fig.       2.2(b).  By      describing      the      hole-electron 

recombination  rate  Urg  along  the  grain  boundary,   assumed  to   involve  only 

band-bound      transitions,       we      can       define      the      effective       surface 

PR 
recombination  velocity  S^,^^x    for  minority  electrons  at  the  edge  of  the 

grain-boundary   space-charge  region. 

At  an  arbitrary  depth  z  in  the  base  region,  we  consider  the  energy- 
band  diagram  in  the  grain-boundary  space-charge  region.  This  diagram, 
which  follows  from  the  solution  of  Poisson's  equation  to  be  discussed, 
is  illustrated  for  non-equil  ibritm  conditions  (V  >  0)  in  Fig.  2.3.  The 
existence  of  the  potential  barrier  height  ^^  is  due  to  majority-carrier 
(hole)  charging  of  the  grain-boundary  surface  states  and  has  been 
discussed  extensively  in   previous  papers  [2,   7-10,   20]. 

To  describe  <^^  of  the  polysilicon  grain  boundary  adjacent  to  a 
p-type  grain,  we  assune  an  "effective"  monoenergetic  density  N^j/^^x  of 
donor-type  grain-boundary  surface  states,  or  traps,  which  are  either 
neutral  or  positively  ionized,  is  representative  of  the  grain-boundary 
potential  barriers  for  p-type  polysilicon  because  the  "neutral  Fermi 
level"  of  polysilicon  grain  boundaries,  which  is  related  to  the  unknown 
nature  of  the  surface  states,  is  experimentally  observed  to  be  near 
midgap  [20].  This  means  that  the  electron  energy  bands  will  bend  down 
at  the  grain  boundaries  in  p-type  polysilicon  (see  Fig.  2.3);  this  will 
be   predicted  by  an  analysis  based  on  donor- type  traps.      Indeed  McT/pff\ 
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is  the  grain-boundary  trap  density  that  would  be  inferred  from 
conductance  or  capacitance  measurements  [2,  20]  of  p-type  polysilicon. 
Although  the  assumption  of  a  monoenergetic  density  of  traps  at  Ej  = 
E^.  (0)  may  not  characterize  exactly  actual  grain  boundaries  in  p-type 
polysilicon  [34],  it  is  sufficiently  representative  to  describe  their 
basic  influence  on  the  minority  electron  transport.  The  utility  of  a 
monoenergetic- trap-density  analysis  in  characterizing  the  electrical 
behavior  of  polysilicon  has  been  previously  demonstrated  [8,  10,  20, 
23]. 

We    recognize     however    that    Ugg,    which    depends    on     ^q,    must    be 

calculated  based  on   the   possible  existence  of  both  donor-type  (Ncj)   and 

a 
acceptor-type  (N^-p)   trap  densities.      We  later   describe   this  calculation 

PR 

and  the  resulting  description  of  S  ,^«,  which  depend  critically  on  the 
large  difference  between  the  carrier  capture  parameters  of  neutral  (C") 
and  ionized  (cM  traps  [9].  We  also  later  discuss  qualitatively  the 
effects  of  moving  Ej  away   from  midgap. 

For  steady-state  conditions,  i.e.,  (j)g  time-independent,  the  energy- 
band  diagram  in  Fig.  2.3  and  Ugg,  which  defines  the  rate  of  flow  of 
electrons  and  holes  to  the  grain  boundary,  are  determined  by  the 
condition  that  the  grain-boundary  charge  that  produces  overall 
neutrality  in  the  grain-boundary  space-charge  region  be  commensurate 
with  equal  net  capture  rates  for  electrons  and  holes  at  the  grain- 
boundary  surface.  Mathematically,  the  quasi-equil  ibrium  assumption 
enables  us  to  express  this  condition  by  combining  Poisson's  equation  fur 
the  grain-boundary  space-charge  region  with  Shockl ey-Read-Hall  capture- 
emission    statistics  [35]    for  the  grain-boundary   surface    states. 
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Referring  to  Fig.  2.3,  in  which  the  grain  boundary  is  located  at 
;  =  0,  we  write  Poisson's  equation  in  the  space-charge  region 
(0  <  X  <  ^11^)   as 


d^E  (x)         2 

L_=l-[P(x)   -  N(x)   -  N     ] 


(2.1) 


where  P  and  N  are  the  hole  and  electron  densities,  N^^  is  the  grain 
(acceptor)  doping  density,  q  is  the  electron  charge,  e  is  the  dielectric 
constant  of  silicon,  and  E^-  is  the  intrinsic  Fermi  level,  whose 
derivative  with  respect  to  x  defines  the  electric  field.  The  electric 
field  at  x  =  0  i  s  determined  by  the  grain-boundary  surface- state  charge 
through  Gauss'  law  and  provides  one  of  the  two  boundary  conditions 
necessary  to   solve  (2.1): 


^ST(eff) 


x=0 


2e 


(2.2) 


where  f  is  the  steady-state  probability  of  electron  occupancy  for  the 
surface  states.  The  factor  of  1/2  in  (2.2)  reflects  our  implicit 
assumption  of  symmetrical  energy  bands  on  either  side  of  the  grain 
boundary,   i.e.,   identical    adjacent  grains. 

Assuming    a    unity    degeneracy    factor    for   the   grain-boundary    surface 
states,  we  get    from  the  Shockley-Read-Hall    theory  [35] 


:2.3] 


which   simplifies  under  our   assumption   that  Ej  =   E^-(O).      In   (2.3),   n^-   is 


C^N(O)  +  C  n.exp 

rE.(0)-E,l 

f  - 

kl 

C^N(O) 

+  C^n.exp 

VS-io)i 

'E    (0)-E   - 
.  CpP(O)  .  Cpn.exp^-i-^ 

L     ^"^ 
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for  electrons,  which  we  assume  is  C''  =  10"'^cm-^/sec  (at  300°K)  because 
the  trap  is  (positively)  ionized  prior  to  the  electron  capture  [36];  and 
C       is     the     capture     coefficient      for      holes,      which     we     assume      is 


C 


^-9_3 


10"^cm-^/sec  (at  300°K)  because  the  trap  is  neutral  prior  to  the 
hole  capture  [36], 

To  write  the  second  boundary  condition,  we  assume  low  injection  in 
the  uni  formly  doped  quasi-neutral    p-region: 


dE. 


X  >  w 


GB  " 
SCR 


(2.4) 


This  assumption  also  means  that 


P^^scR^  ~-  h^ 


(2.5) 


Within  the  grain-boundary   space-charge  region,   P  and  N  are  related  to  E^ 
through  the  respective  qua  si -Fermi   levels  Epp  and  E^^: 


P(x)  =  n^.exp 


rEi(x)-Epp 


TT 


(2.6) 


N(x)  =  n^.exp 


Ep,-E.(x) 


TT 


(2.7) 


The  difference  between  Ep^^  and  Epp  reflects  the  "excitation"  of  the 
grain  boundary  that  results  from  the  forward  bias  on  the  n'^p  junction. 
To  facilitate  the  separation  of  the  grain-boundary  recombination 
analysis  from  the  intragrain  transport  analysis,  we  assume  that 
conditions  of  quasi-equil  ibrium  exist  in  the  grain-boundary   space-charge 
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region.  This  assumption  implies  that  Epj^  and  Epp  are  nearly  flat  (i.e., 
that  their  variations  are  less  than  kT)  in  this  region,  and  will  be 
valid  if  N  and  P  are  sufficiently  large  and  the  electron  and  hole 
currents  sufficiently  small,  as  we  discuss  later.  This  will  be  the  case 
for   sufficiently   high  excitations  (see  Fig.   2.3), 


qVg  A  Ep^  -   Epp        .  (2.8) 


and  sufficiently  low  potential   barriers, 


q<t.g   A   E.(W^^j^)    -  E.(0)        .  (2.9) 


Note  that  Vg  defined  in  (2.8)  decreases  with  z  in  accordance  with 
N(rQ,z),  and  that  necessarily  Vg(z)  is  everywhere  less  than  V  because  of 
the  electron  density  gradient  in  the  grain. 

As  we  will  see,  the  quasi-equil  ibri  in  assumption  enables  us  to 
characterize  Ugg  without  having  to  describe  the  mechanism,  e.g., 
thermionic  emission  [10],  by  which  majority  carriers  (holes)  are 
transported  from  the  quasi-neutral  region  to  the  grain-boundary 
surface.  This  simplification  is  made  without  loss  of  generality  because 
when  the  quasi-equil  ibrium  assumption  is  invalid,  for  example,  when  the 
grain-boundary  space-charge  region  is  virtually  depleted  of  free 
carriers  [7,  9,  10],  the  minority  electrons  are  approaching  the  grain- 
boundary  surface  with  velocities  about  equal  to  their  kinetic  limit 
S^,_  (=  5  X  10^  cm/ sec  at  300°K)  [31].  Hence,  as  we  show  later,  when 
(2.8)  becomes  inapplicable,  the  grain  boundary  can  be  adequately  modeled 
as  a    surface    having   a   recombination   velocity   equal    to  Si/,  .     A  detailed 
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analysis  of  grain-boundary  recombination  based  on  the  depletion 
approximation   thus  is  unnecessary  and,   in    fact,   is  invalid. 

The  solution  of  (2.1)-(2.7),  with  arbitrary  excitation  Vg  in  (2.8), 
is  impossible  to  derive  analytically.  Thus  to  aid  the  development  of  an 
analytic  model  for  S^^^^^  in  terms  of  Vg  and  the  polysilicon 
properties,  we  use  computer-aided  nimerical  techniques.  A  Harwell 
subroutine  [37],  which  represents  (2.1)  using  finite-di  fference 
approximations  for  the  derivatives,  has  been  slightly  modified  to  solve 
(2.1)-(2.9).  Details  of  the  numerical  algorithn  and  the  listing  of  the 
computer  program  are  given  in  Appendix  A.  We  now  discuss  results  of 
this  numerical    analysis  of  the  grain-boundary  region. 

Consider     first    the    thermal -equi  1  ibri urn    case    (V    =    0)     for    which 

Vg  =  0     and     Ep,^    and     Epp     are    coincident    with    the     Fermi     level     Ep. 

Equilibrium    barrier    heights    <figQ    at    polysilicon    grain    boundaries    have 

been   previously  calculated  [7]   and  measured  [20,   34,   38],   and  comparison 

of  our  nunerical    results  with  these  earlier  ones   provides  an   assessment 

of  the  assumptions  underlying  the  analyses.     For  example,   it  implies  the 

validity  of  the   depletion  approximation    for  equilibrium  conditions.      In 

Fig.   2.4  we   have   plotted  <\>^q  versus  N^^  for  representative  values  of  the 

effective  grain-boundary   surface-state  density  N°  ,    .«  .     When  N°  /    rrs 

ST(eff)  ST(eff) 

10         2 
is  \iery   low  (<  10       cm"  ),  4,^^  is  extremely   anal! ,   and  electrically  the 

polysilicon  resembles  single- crystal  silicon.  For  intermediate  values 
°^ '^ST(eff)  ^~  ^^  ^^  ^'  *B0  decreases  monotonically  with  increasing 
n^f^.  This  is  explained  by  the  fact  that  ^^q  is  gnall  enough  that  Ep  is 
below  Ej  (see  Fig.  2.3),  and  hence  virtually  all  of  the  (donor)  surface 
states  are  empty  (positively  ionized)  for  all  N^^^.  Thus  as  N^^ 
increases,      the      space-charge     region     width,     which     contains     charge 
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(including  ionized  acceptor  impurities)  that  neutralizes  the  grain- 
boundary   charge,    diminishes   and    hence   (t)gQ   decreases.      When  N^,,    -:«   is 

12       -2 
high  (>  10      cm     ) ,      ^^q   is    large   enough    that   Ep   is   essentially    fixed 

at,   or   slightly  above,    Ey   (see   Fig.    2.3).      Thus  as  N^^^  increases  and  Ep 

moves  closer  to   the  valence-band  edge   in   the   quasi-neutral    grain  [i.e., 

farther  away   from  E^-    as  described  by   (2.5)    and   (2.6)],    Ej  follows  Ep  at 

X  =  0,   and  hence  (j)gQ  increases  slightly. 

We  note  that  the  calculated  values  of  (jjgQ  in  Fig.  2.4  are  generally 
higher  than  those  derived  from  measurements  [34]  of  the  activation 
energy  for  the  zero-bias  conductance  of  grain  boundaries  in  p-type 
polysilicon.  This  discrepancy  results  because  of  our  assumption 
Ej  =  E^. (0),  which  is  not,  as  we  mentioned  earlier,  strictly  valid. 
Assuming  Ej  <  E^-(O)  brings  the  <^^q  predictions  and  measurements  into 
better  agreement,  but  does  not  significantly  augment  the  physical 
insight  provided  by  our  analysis.  It  is  interesting  to  note  that  for 
n-type  polysilicon,  analogous  calculations  based  on  an  assumed 
monoenergetic  density  of  acceptor-type  traps  at  midgap  yield  values  of 
(t)gQ  that  are  roughly  consi  stent  with  those  derived  experimentally  [38]. 

Consider  now  non-equil  ibriun  cases  (V  >  0)  fbr  which  Vg  >  0  in 
accordance  with  the  quasi-equil  ibrium  assumption  in  the  grain-boundary 
space-charge     region.  Non-equil  ibriim     conditions      fbr     which     this 

assumption  fails,  alluded  to  earlier,  will  be  treated  explicitly 
later.  The  relationship  between  Vg(z)  and  V  is  complex  and  follows  from 
the  electron  continuity  equation  in  the  quasi-neutral  p-region  of  the 
grain,  which  we   solve  in  the  next  subsection. 
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2.4  Calculated  thermal-equilibrium  barrier  height  versus  grain  doping 
density  (p-type)  for  representative  values  of  effective  grain- 
boundary  surface-state  density  (T  =  300°K). 
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Figs.  2.5  and  2.6  show  plots  of  calculated  barrier  heights  <^n 
versus  Vg  for  representative  values  of  N<-j/  ^^  and  for  N,.  =  10^^  cm"'^ 
and  10^'  cm"-^  respectively.  These  plots  are  restricted  to  conditions  of 
low  injection   in  the  quasi-neutral    p-region.     They  imply  that 

3vr=-i  (2.10) 


for  Vg  greater  than  a  few  thermal  voltages  (kT/q  =  0.026  V),  a  condition 
that     validates     the      quasi-equil  ibrium      assumption.  The     monotonic 

reduction  of  <^q  with  increasing  Vg  is  consistent  with  the  results  in  [8] 
and  is  qualitatively  explained  as  follows.  The  electron  capture  rate  at 
the  grain  boundary,  which  is  proportional  to  N(0),  increases  with 
increasing  excitation.  Consequently  the  net  charge  at  the  grain 
boundary  and  hence  the  grain-boundary  barrier  height,  which  is  directly 
proportional  to  the  net  charge,  decrease  with  increasing  excitation. 
The   slope  of  -1/2  in   (2.10)    has  a   physical   basis  as  we  now  describe. 

For  steady-state  conditions  with  sufficient  excitation,  the 
electron  and  hole  capture  rates  at  the  grain-boundary  surface  states, 
including  both  donor  and  acceptor  types,  are  nearly  equal.  Thus,  for 
either  type  of  state,  we  can  write 


af-CC^d   -   f)N(O)]  =  g|-[CpfP(0)]       .  (2.11) 


From  (2.3),  for  the  assumed  conditions,  we  see  that  f  =  C^M(0)/[CpN(0)  + 
CpP(O)],  which  is  always  between  zero  and  one.  Because  the  change  in  f 
induced  by  a  small  variation  in  Vg  will  be  much  smaller  than  the 
corresponding  changes  in  N(0)   and  P(0),   (2.11)    implies 
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Fig.  2.5  Calculated  non-equilibrium  barrier  height  versus  grain- 
boundary  excitation  for  representative  values  of  effective 
grain-boundarv  surface-state  density  and  for  a  grain  doping 
density  of  10   cm"^  (T  =  300°K). 
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1      3N(0)  1      8P(0)  ,,    i,x 


The  relation  in  (2.12)  is  a  general  result  of  invoking  the  steady-state 
condition  (d<j)g/dt  =  0)  at  a  sufficiently  excited  grain-boundary 
surface.     Combining   it  with 


P(0)N(0)  =  n^expj-^l      .  (2.13) 


which   follows   from  {2.6)-(2.8),  reveals  that 


1     3P(0)         q 


Then   substituting 


(2.14) 


P(0)  =  N,,  exp(:J      '  (2.15) 


which  results  from  (2.5),  (2.6),  and  (2.9)  with  the  assumption  that  Epp 
is  nearly  flat,  into  (2.14)  yields  (2.10).  (Note  that  this  analytic 
description  of  the  grain-boundary  physics,  as  well  as  others  to  follow, 
are   facilitated  by  the  numerical    solutions.) 

The  important  conclusion  to  be  drawn    from  (2.10)   and  the  associated 
relationships  mentioned  above  is  that 


P(0),   N(0)   cc  exp[-^|      ,  (2.16) 
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and  consequently,  as  we  detail  later,  Ugg  «  exp(qVg/21<T) .  We  now 
discuss  the  dependencies  of  (fig,  P(0),  and  N{0)  on  n^j,  rry.  and  on  N^^  as 
derived   from  the  numerical    solutions  of  (2.1). 

These    solutions,    those    plotted    in    Figs.    2.5    and    2.6    and    others, 
imply    further  that 


for    sufficient  excitation   (Vg)    as  defined  earlier.      The   implications  of 
(2.17)   regarding  P(0)    and  N(0)    follow   from  (2.13)    and  (2.15): 


~  r  n  ,-'log(e)  -1/2 

P(0^  '  [^ST(eff)l  =   KT(eff)l  ^^.IS) 


and 


~   r  n  Jog(e)       ^  ^  1/2 

^(0^   "   KT(eff)l  ~-   KT(eff)  '  ^^.19) 


which  are  also   seen  directly   from  the  ntnerical    results. 

A  physical  explanation  for  (2.18)  and  (2.19)  can  be  given,  again 
based  on  the  condition  of  nearly  equal  electron  and  hole  capture  rates 
through  the  (effective)   grain-boundary  traps: 


^n^l   -   ^^ST(eff)^(°^   =  S^ST(eff)P(0^        '  ^^.20) 


For  sufficient  Vg,  partly  because  of  the  large  discrepancy  between 
C^  (  C"")  and  Cp  (=  c"),  f  =  1;  i.e.,  most  of  the  traps  are  filled  with 
electrons.       Therefore    the     factor  (1  -    ■'')N5y/g^^^    in    (2.20)    must    be    a 
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subl  inear     function    o^Ncy/pff);   i.e.,     f    must    increase     slightly    with 


increasing  NsT(eff)*     ^^^^  (2.20)    suggests  that 
N(0)  ~  ^D 


Vm      "ST(eff) 


(2.21) 


Then     (2.13)     and    (2.21)     give    approximately    the    observed    dependencies 
(2.18)   and  (2.19). 

Finally,   the  nunerical    solutions  of  (2.1)    for   sufficient  excitation 
reveal    that 


d(|) 


B  3kT 


8LlogN^^]  "     2q 


(2.22) 


and  that 


^     [1  -i.log(e)]         ./^ 


(2.23) 


and 


N(0)  ^  H 


-Cl-|log(e)]         _^/3 


AA 


'AA 


(2.24) 


These   approximate   dependencies  can  be   supported  by   physical    argument  as 
follows.     From  (2.15),   we  get 


1  - 


ap(0) 


(2.25) 


For   <|,g   >   0,    the    second   term  within   the  brackets  of  (2.25)    is  less  than 
one,    and,    typically,    considerably    so   as  evidenced   by    (2.15),    Fig.    2.3, 
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and  the  exponential  relationship  between  P  and  Epp,  which  defines  the 
modification  in  the  energy-band  diagram  resulting  from  a  change  in 
N^^.     Thus  (2.25)    implies  that 

which  is  approximately  equivalent  to  the  observation  (2.22). 

Our  computer-aided  detemiination  of  the  dependencies  of  (^^  and  of 
P(0)  and  N(0)  on  Vg,  NcT(pff).  and  N^^  enables  the  characterization  of 
the  steady- state  hole-electron  recombination  rate  Ugg  through  the  grain- 
boundary  surface  states,  provided  the  active  surface- state  densities  are 
specified.  This  calculation  is  based  on  the  Shockl  ey-Read-Hall  capture- 
emission  model  [35]  fbr  recombination-generation  through  localized 
states  in  the  energy  gap. 

Because    the    exact   nature   of  the    grain-boundary    surface    states   in 

polysilicon    is   not   known,    we    allow    fbr    the    possible   existence  of  both 

donor- type  (Ncj)  and    acceptor- type  (^^T)   states    in    our    calculation    of 

Ugg.       We    assume    that    all     states    are    located    at    E-r    =    E^-(O).       Note 

that  M^,  and  N_,  are      different       from  Nct/    r^-i  >  which      was      used      to 
ST  ST  ST(eff) 

calculate     (j)g,      and     P(0)      and     N(0).  Physically,     based     on     this 

D  A 

model,  N<--p  and  H^j  produce    grain-boundary    charge     that    is    effectively 

characterized  by  Ncj/girrv- 

If  we  neglect  bound-bound  carrier  transitions  between  the  donor  and 

D  A 

acceptor   traps,    the  composite   recombination   rate   through  N<~j  and  N^j  is 

simply 


'gB   =   ^GB^^ST^   '  'lAl^  (2.27) 
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D  A 

where  U^g   and  U^g   are     Shock! ey-Read-Hall      representation     [35]     of    the 


recombination  rates  through  the  donor  and  acceptor  traps  respective! 


P                 P(0)N(0)-n? 
^GB  "  P(0)  +  n. M(0)   +  n.        '  ^2.28) 

C  ^3^  C   N3^ 


^                P(0)N(0)-n? 
^GB  =   P(0)   +  n, m)   +  n.      '  (2.29) 


'   + 


cX7~     C^"<T 


In  (2.28)  and  (2.29)  we  have  used  the  previously  assumed  capture 
parameter  values  C^  =  10"^  cm-^/sec  and  c"  =  10"^  cm^/sec  for  the  ionized 
(positively  or  negatively)    and  neutral    states  of  the  traps. 

The  dependencies  of  Ugg  on  the  excitation  Vg  and  on  the  grain- 
boundary  surface-state  densities,  ^JjKeff)'  ^ST'  ^"^  ^ST'  ^"^  °"  ^^^ 
grain  doping  density  N^^  are  implicitly  given  by  (2.27)-(2.29)  and  the 
explicit  dependencies  of  P(0)  and  N(0)  derived  from  the  solutions  of 
(2.1).     For  example,    (2.16)   and  (2.27)-(2.29)    show 


qVr 
^GB^^^pI^I  (2.30) 


for  the  sufficient  excitations  needed  to  formulate  the  analytic 
relationships  discussed  earlier.  The  result  (2.30)  agrees  with  [7]  and 
[8],  but  is  not  equivocal  like  [7]  and  [8]  because  of  questionable 
assumptions,  which,  in  fact,  are  shown  to  be  invalid  by  our  numerical 
results.      Furthermore   (2.30)    has  been   derived   here    for   the   general    case 
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in  which  both  donor-  and  acceptor-type  traps  are  present  at  the  grain- 
boundary  surface,  and  the  constant  of  proportionality  for  (2.30)  is 
related  correctly  to  the  assumed  model  for  the  surface  states,  in 
contrast  to  [7]   and  [8]. 

If  we  neglect  recombination  in  the  grain-boundary  space-charge 
region  (see  Fig.  2.3),  Ugg  defines  the  effective  recombination  velocity 
for  minority  electrons  at  the  edge  of  the   space-charge  region  [7,   8]: 


^(eff) 

__           ^GB 

2^(WSCR^ 

where 

,    from   (2.5)-(2.8), 

mil,) 

"i            ^^B 

(2.31) 


:2.32: 


,GB 


Vg      and  N(W^.pn),  both      dependent     on 


z,      define      the      grain-boundary 
excitation   and  are  determined    from   the  coupled   solution   to   the  electron 


continuity   equation   to  be   discussed   in    the    following    subsection.      Note 
the  nonlinearity  associated  with  S^,    ..x   as  defined  by  (2. 30) -(2. 32): 


Sn{eff)   =^o[^ST(eff)'^ST'^ST'^AAl^^p(- 


qVc 


WTj 


hMlp? 


1/2 


(2.33) 


where  the  pre-exponential  coefficient  Sq  reflects  the  functional 
dependencies  derived  from  the  numerical  solutions.  An  expression  for 
Sq,  derived  empirically   from  the  numerical    solutions,  is  given  by 
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C"C^N 


So  = 


AA 


'ST 


2n. 


^ 


°92   /5xlO 

K,0 


^0^2       .-/lois^iogz  ,     , 


'ST 


'AA 


D    \m 

ST 


5x10 


10 


ST 


l5\1og2 


^ST 


T7? 


J^^] 


(2.34) 
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Values    0^  S  ,    -:Jr^    have    been    numerically    calculated    for    two    cases 

In   the    first  case,  we  let 


involving  assumed  trap  densities  N^^  and  N 


ST 


ST' 


N^-j.  =  N^-p,    ^«   and  N^^  =  0,   which     corresponds     to     the     simple     common 
assumption    that    only    donor-type    traps    are    effective    on    a    polysilicon 

PR 

grain    boundary   between    r>-type   grains   [23].      Plots  of  SjJ?    ^^x   versus   Vg 
for    this   case    are    shown    in    Figs.    2.7    and   2.8    for   N^^^   =    10^-   cm 


^15   .^-3 


and 


10^'    cm"-^    respectively.       In    the    second   case,  we    let  N^,  =  2N'    ,    .^s 

A  D 

and  Nj^  =  H^^,    ^«    to     approximately     represent     a     possible     situation 

involving    both    donor-    and    acceptor- type    traps  at    the    grain    boundary. 


Figs.     2.9     and     2.10     show  S 


GB 
n(eff) 


for     this    case,     and    emphasize    the 


possibility  of  certain  traps  influencing  U^g  but  not  significantly 
affecting  (}ig  or  the  conductance  or  capacitance  associated  with  the  grain 
boundary. 

Superimposed  on  the  plots  in  Figs.  2.7-2.10  is  an  indication  of  the 
kinetic-limit  velocity  \^  (=  5  x  10^  cm/ sec  at  300°K) ,  which  is  defined 
by  the  random  thermal  motion  of  electrons  in  the  conduction  band  [31]. 
The  velocity  S,^i_  is  the  average  directional  thermal  velocity  defined  by 
Maxwell -Boltzmann  statistics  for  conduction-band  (free)  electrons,  and 
is  the  physical  upper  limit  for  S^^^^^^.  We  therefore  recognize  that 
our    plots  of  S^^g^^^   in   Figs.    2.7-2.10   must   be   truncated   at   S.^^-      T^^' ^ 
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recognition  effectively  resnoves  any  uncertainty  or  restriction  of  our 
model  due  to  the  possible  invalidity  of  the  quasi-equil  ibrium 
assumption.     We  now  demonstrate   thi  s. 

The  quasi-equil  ibrium  assumption  is  valid  if  the  variations  in  the 
quasi-Fermi  levels  across  the  grain-boundary  space-charge  region  (see 
Fig.  2.3)  are  less  than  kT  [8,  9].  A  quantitative  self-consistency 
check  for  the  validity  of  this  assumption  is  described  in  [8].  Using 
this  check  and  our  results  shown  in  Figs.  2.7-2.10,  we  find  that 
generally  the  quasi-equil  ibri  in  assumption  is  justified  provided 
^n(eff)  ^'^  ^^^^  ^^^^  ^KL'  Since  S  ,  ^^k  cannot  physically  exceed  S^^, 
our  model  is  then  generally  valid  if,  as  indicated  in  Figs.  2.7-2.10,  we 
stipulate  that 


"«Keff)]    =   ^KL       •  <2-35) 


Hence  for  all  conditions  under  which  the  grain-boundary  space-charge 
region    is    depleted    of    free    carriers,    which    obviously    preclude    quasi- 

PR 

equilibriin,  S  ,  -,v  =  Sj,.  ,  and  a  depletion-approximation  analysis 
[7,  9,  10]  is  unnecessary.  In  fact,  a  recent  such  analysis  [10]  yielded 
results  that  comply  with  these  conclusions. 

Note    in    Figs.    2.7-2.10    that    for   low  values  of  N?-r,    ^«  ,   depending 

TR 
on    N^^,  S^,g^^v   as  calculated    from   the   ninerical    solutions  of  (2.1)    is 

insensitive    to    low    excitations    Vg.        Our    analytic    model     (2.33)     is 

inapplicable      fbr  these     conditions  because  of  the   underlying  assumption 

of      "sufficient      excitation."  This      insensitivity      reflects      the 

insensitivity  of  <|)g  on  Vg   fbr  low  MJ-|-,g^^^      illustrated  in  Figs.   2.5   and 

2.6.       Fig.   2.4     reveals     that     (|)g  =  <^^n^     for  these     cases  and     explains 
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Calculated  effective  recombination  velocity  for  minority  electrons 
at  the  edge  of  the  grain-boundary  space-charge  region  versus  grain- 
boundary  excitation  and  surface-state  density.  The  kinetic-limit 
velocity,  which  is  the  physical  uoper  limit  for  the  recombination 
velocity,  is  indicated  here  and  in  Figs.  2.8  -  2  10 
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Fig.  2.8  Calculated  effective  recombination  velocity  for  minority  electrons 
at  the  edge  of  the  grain-boundary  soace-charge  region  versus  grain- 
boundary  excitation  and  surface-state  density. 
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Fig.  2,9  Calculated  effective  recombination  velocity  for  minoritv  electrons 
at  the  edge  of  the  grain-boundary  space-charge  region  versus  arain- 
boundary  excitation  and  surface-state  density 


Fig,  2,10  Calculated  effective  recombination  velocity  for  minority  electrons 
at  the  edge  of  the  grain  boundary  space-charge  region  versus  grain- 
boundary  excitation  and  surface-state  density. 
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PR 
why  S^^g^^    is     independent     of    the     grain-boundary     excitation.         The 

excitation  is  insufficient  to  create  a  high  enough  electron  density  to 
alter  the  barrier  height  from  its  low  equilibriun  value  [N(0)  << 
P(0)].  Thus,      when      the     electron      emission      rate      is     negligible, 

^n(eff)  -  '^n^ST     ^"^  ^^^'^  "^^^• 


2.2.2     Grain 

We     now    derive     an     approximate     solution     for    the     two-dimensional 

electron    continuity    equation    in    the    quasi-neutral     p-base   of  the    grain 

illustrated     in      Fig.      2.2(b).  This     solution,     which     is     governed 

PR 
by  ^^iQff)   given       in       (2.33)       and       (2.35),       yields       the       electron 

op  p 

current  I^  =  I^g  +  l^^  injected   into    the   base  of  the    forward-biased  n^p 
junction. 

Because  of  the  nonlinear  boundary  condition  at  the  grain 
boundaries,  a  rigorous  derivation  of  this  solution  would  require  a 
computer.  To  simplify  the  problem  such  that  an  analytic  solution, 
useful  for  engineering  design  and  for  demonstration  of  the  pertinent 
physics,  can  be  obtained,  we  partition  the  i>region  as  shown  in 
Fig.  2.11.  The  implicit  assumption  made  in  doing  this  is  that  in 
certain  subregions  of  the  base,  the  electron  flow  is  predominantly  one- 
dimensional;  that  is,  the  divergence  of  the  current  density  has  one 
predominant  term.  Thus  the  resulting  one-dimensional  forms  of  the 
continuity  equation  can  be  solved,  and  the  solutions  then  coupled  to 
give  N(r,z),  which  defines  I^.  This  simplification  is  similar  to 
Shockley's  "gradual-case"  approximation  in  his  analysis  of  the  unipolar 
transistor  [32]  and  was  qualitatively  mentioned  in  [23].  Its  utility 
was  established  in  a  numerical  study  [39]  of  the  two-dimensional 
minority-carrier    flow  in  a    forward-biased   planar  diffused  pn  junction. 


41 


As  indicated  in  Fig.  2.11,  we  assume  that  in  the  central  portion  of 
the  grain  (r  <  r-^) ,  the  electron  flow  is  not  influenced  by  the  grain- 
boundary  recombination,  and  is  hence  in  the  downward  z-direction. 
Contrarily,  we  assume  that  in  the  vicinity  of  the  grain  boundary  (r,  < 
r  <  rQ),  the  electron  flow  is  strongly  influenced  by  the  grain-boundary 
recombination,  and  is  thus  primarily  in  the  lateral  r-direction.  These 
assumptions  are  commensurate  with  letting 


(2.36) 


where   L^    is   the   electron    diffusion   length.      Note    that   (2.36)    restricts 

our  quasi-two-dimensional    analysis  to  cases    for  which  r^  >  L^.     Although 

B  R 

this    inequality    implies   that  Ig  is   not   dominated   by  Igg,  the    solutions 

for  these  cases  nevertheless  serve  to  illustrate  well    the  effects  of  the 

PR 

nonlinearity  introduced  by  S  ,  \|:« .  If  rg  <  L^,  which  implies  a 
predominance  of  I^g,  the  electron  flow  is  truly  three-dimensional 
everywhere,  and  this  case  requires  a  computer  solution.  However,  from  a 
practical  viewpoint,  this  case  may  be  undesirable  because  the  electron 
current  is  dominated  by  grain-boundary  recombination  (I^  >  I  );  it 
could    be    avoided    perhaps    by    increasing    N^y^    to    reduce    L^,    while    still 

D 

decreasing   Ip. 

In    the   central    region   of  the   grain,    the   electron   density   N(z)    for 
low- injection  conditions  is  defined  by 


^-Jiif.O  (2.37) 


dz' 


with  the  boundary  conditions 
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Fig.  2.11  Partitioned  p-region  of  the  representative  grain  showing  predom- 
inantly one-dimensional  electron  flow  in  the  subregions  {r^>l   ) 
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«0).^exp(^) 


(2.38) 


and 


N(W„)   =  0 


(2.39) 


the  latter  being  due  to  the  olnic  contact  on  the  back  of  the  junction. 

The  solution  of  (2.37)-(2.39)  which  is  obvious,  then  provides  one 
boundary  condition,  i.e.,  H{r-^,z) ,  for  the  coupled  problem  in  the 
vicinity  of  the  grain  boundary: 


r  ar 


9N(r,z) 

3r 


N(r,z] 
~T2~ 


(2.40; 


The  general  solution  of  (2.40)  is  an  infinite  series.  To  obtain  an 
approximate  solution  that  is  consistent  with  our  parti tioned-base  model, 
we  initially  neglect  the  second  term,  but  later  account  for  the 
recombination  current  it  defines  by  integrating  qN(r,z)/Tp  (  t^  is  the 
electron  lifietime)  over  the  qua  si -neutral  base  of  the  grain.  This 
results  in  a  ln(r)   dependence    for  N. 

The     second    boundary    condition     for    (2.40)     is    that    defining    the 
grain-boundary  recombination  discussed  in  the   previous   subsection: 


q     N     G 


-D. 


9N(r,z] 
3r 


ilfn'^'G^ 


(2.4i: 
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where  D^  is  the  electron  diffusion  coefficient  and  where  S^J^^  >,  which 
depends  on  z,  is  given  either  by  (2.33),  in  which  N(W^^  )  -  N(r  ,z) ,  or 
by  (2.35).  The  determination  of  which  expression  to  use  involves  an 
iterative  process.  We  initally  use  (2.33),  solve  (2.37)-(2.41) ,  compare 
(2.33)  and  (2.35)  everywhere  on  the  grain  boundaries,  and  then  use 
(2.35)      where     necessary.  This     complicated      process     is     typically 

simplified  because  unrealistic  values  of  S^^^^^^  >  S^^  i-^Ply,  through 
the  solution  of  the  electron  continuity  equation  in  the  grain,  that 
Vg  =  0  at  the  point  being  considered;  this  is  also  implied  by  the 
realistic  value  of  S^^^^^^  =  s^^  ^'^'^^'^  hi  >  V^n »  which  is  common. 
The    electron    density   N(r,z)    derived    from    (2.37)-(2.41)    can    now  be 

B 
G" 


used    to   calculate   the    injected   electron   current  I^^.     The  component  due 
to  recombination  at  the  grain-boundary   surfaces  is 


4   =   -2.qDr,     ^B  iNllr^ 


dz       ,  (2.42! 


and  that  due  to  recombination  in  the  quasi-neutral    p-base  i 


p-base  is 


tB     _   2TTq       ^B     "^G     ,,,        ,        ^     , 
^ON  "     T         J       /       N(r,z)    r  dr  dz 
n     0       0 


(2.43) 


in     which     the     integration     must     be     done     in     accordance     with     the 
partitioning  of  the   p-region. 
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The    general     form    of    the    resulting    expression     for    the    electron 
current,    (I^g  +  l^^) ,  is 


In   (2.44)   the  exp(qV/kT)  component  is  the  electron  current  that  would  be 

injected   in   the   absence  of  grain   boundaries.      For  example,   i  f  L^  <   Wg, 

2  2 

then   Iq^  =  ^G^^i^n^'^AA^n  ^^^""^  ^q,  ^  ^"^q  ''*  ^    ^^^    ^^^^    °^  ^^^    grain.       In 

the  exp{qV/nkT)  component,  which  in  general  comprises  both  grain- 
boundary  and  intragrain  recombination,  Iq  depends  on  Sq,  given  in 
(2.34)  in  addition  to  the  p-region  parameters  that  define  Iqt^.  The 
reciprocal  slope  factor  n(V),  which  also  depends  on  Sq  and  the  p-region 
parameters,  increases  from  one  to  two  as  V  increases,  albeit  in  the 
range  corresponding  to  low-injection  conditions  in  the   p-region. 

The  general  shape  of  the  Ig(V)  characteristic  in  (2.44)  for  typical 
values  of  the  grain  and  grain-boundary  parameters  is  illustrated  in 
Fig.  2.12.  The  shape  resembles  that  of  the  dark  current-voltage 
characteristic  of  a  BSF  solar  cell  [40],  e.g.,  a  p^nn"^  cell,  in  which 
the  inflection  in  the  curve  results  because  of  the  relative  significance 
of  different  components  of  current  having  different  reciprocal  slope 
factors. 

In    Fig.    2.12,     for   low   V,    I^  «  exp(qV/kT)  because    n(V)   =  1.      This 

occurs    because    the    electron    injection    level     is    very    low,    and    hence 

^n(eff)'  "^'"^^      '^      proportional       to  [N(W^^p)]"^/2   (=[N(rg,z)]"^/2  as 

shown    by    (2.33),    is    high    and    fixed    at   S|^|_   everywhere    (see    Figs.    2.7- 

R  R 

2.10).     The  nonlinearity  is  thereby  removed  and  both  1°     and  I^j^  vary  as 

PR 
exp(qV/kT).        As    V     increases,  N(W^^j^)       tends    to     increase    also,     and 
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thus  S^^g^^j  begins  to  fall  below  S^^  at  portions  of  the  grain 
boundaries.  The  nonlinearity  is  hence  apparent,  and  n(V)  exceeds  unity; 
the  slope  of  Ig(V)  decreases  as  shown  in  Fig.  2.12.  At  higher  V,  the 
exp(qV/nkT)  component  becomes  insignificant,  and  Ig(V)  -►  I^,exp(qV/kT) , 
resulting  in  the  inflection  in  the  curve,  above  which  the  grain 
boundaries  are  inconsequential  because  S^^^^^  is  low  where  the  electron 
density  is  appreciable. 

The  general  shape  of  the  Ig(V)  curve  predicted  by  our  analysis  and 
shown  in  Fig.  2.12  is  consistent  with  measured  current-voltage 
characteristics  of  gated  n'^p  diodes  fabricated  on  Wacker  polysilicon 
[16].  The  characteristics  of  diodes  containing  substantial  grain- 
boundary- surface  area  show  an  inflection  similar  to  that  in  Fig.  2.12 
whereas  those  of  diodes  having  little  or  no  grain-boundary- surface  area, 
e.g.,   diodes   fabricated  within  a  grain,    shown  no  inflection. 

If  Sq  ^  0,  then  Iq^  >  0,  and  the  Ig(V)  characteristic  approaches 
that  of  a  single-crystal  silicon  junction.  If  S^^^^^,  ^  S^^  everywhere, 
then  n  >  1,  and,  as  in  the  Sq  =  0  case,  Ig(V)  =c  exp(qV/kT),  but  with  a 
pre-exponential  coefficient  greater  than  Iq^  because  of  the  grain- 
boundary  recombination.  This  discrepancy  is  commonly  described  in  terms 
of  an  "effective"  minority  electron  lifetime  in  the  grain,  which  is 
smaller  than  x^  [1,  8,  13].  It  is  important  to  note  however  that  such  a 
description  pertains  only  to  I^,  and  not  necessarily  to  other  responses 
of  the  junction,  e.g.,  the  short-circuit  photocurrent,  and  is  meaningful 
only  when  n  =  1. 
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Fig.  2.12  Theoretical  representation  of  the  electron  current  injected  into 
the  base  of  the  forward-biased  polysilicon  n"''p  junction. 
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2.3     Di  scussion 


In  this  chapter  we  have  developed,  using  computer-aided  nimerical 
solutions,  an  approximate  model  for  the  electron  current  Ip  injected 
into  the  base  of  a  representative  grain  of  the  forward-biased 
polysilicon  n^p  junction  illustrated  in  Fig.  2.1.  Key  assumptions  have 
been  made  to  simplify  the  general  three-dimensional,  nonlinear  boundary- 
value  problem,  i.e.,  the  electron  continuity  equation,  and  to  make  this 
development  possible.  We  now  discuss  these  assumptions  and,  where 
necessary,   the  model    limitations  they  imply. 

We  based  the  grain-boundary  recombination  analysis  on  the 
assumption  of  quasi-equil  ibrium,  i.e.,  nearly  flat  quasi-Fermi  levels, 
in  the  grain-boundary  space-charge  region.  This  assumption  facilitates 
the  separation  of  the  grain-boundary  recombination  analysis  from  the 
intragrain  electron  transport  analysis,  the  results  of  which  define  the 
grain-boundary  excitation.  The  assumption  further  enables  a  complete 
analysis  of  the  grain-boundary  recombination  without  having  to  ascertain 
the  mechanism,  e.g.,  thermionic  emission  [10],  by  which  majority  holes 
are  transported  from  the  quasi-neutral  grain  to  the  grain-boundary 
surface. 

The  limitations  associated  with  the  quasi-equil  ibrium  assumption 
are  effectively  removed  by  recognizing  that  when  conditions  obtain  that 
negate  quasi-equil  ibrium,  the  minority  electrons  flow  to  the  grain- 
boiindary  surface  with  velocities  about  equal  to  the  kinetic-limit 
velocity  Skl  [31].  J^^^^  ^g  ^^^^-^^  truncated  in  (2.35)  our  quasi- 
equilibrium-based  prediction  (2.33)  for  the  effective  electron 
recombination  velocity  S^^^^^   at  S^l  when   (2.33)  yields  S^^^^^^    .  S^^_. 
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Even     with     the     quasi-equil  ibrium     assumption,     the     grain-boundary 

analysis,      i.e.,      the      solution      of     Poisson's      equation      (2.1),      is 

formidable.     Thus  we   resorted   to   a  computer-aided  nuTierical    solution  of 

(2.1)    to    facilitate    the    derivation    of  the    analytic   model     for  S^?    .... 

n(eff) 

Empirical      relationships     implied     by      the     numerical      solutions     were 
supported  by  qualitative  discussion  of  the   underlying   physics. 

We  also  assumed  monoenergetic  densities  of  surface  states  at  the 
grain  boundary.  First,  in  accord  with  the  position  of  the  "neutral 
Fermi  level"  of  polysilicon  grain  boundaries  being  near  the  middle  of 
the  energy  gap,  we  postulated  the  existence  of  effective  donor- type 
(adjacent  to  p-type  grains)  states  near  midgap  to  calculate  the  grain- 
boundary  potential  barrier  height  ^g.  Then,  based  on  our  assumed 
existence  of  both  donor-type  and  acceptor- type  traps  at  midgap,  <j)g  was 
used     to     calculate     the     grain-boundary     recombination     and     ultimately 

^n(eff)' 

The  utility  of  an  analysis  based  on  a  monoenergetic  density  of 
grain-boundary  surface  states,  which  is  not  realistic,  has  been 
discussed  before  [8,  10,  20,  23].  Such  a  model  provides  physical 
insight  and  possibly  simulates  well  actual  surface- state  distributions 
in  the  energy  gap  that  significantly  influence  carrier  recombination  at 
the  grain  boundary.  A  recent  analysis  [11],  which  assumes  the  surface- 
states  to  be  distributed  in  the  energy  gap,  has  yielded  results  that  are 
consistent  with  ours.  To  generalize  our  results  somewhat,  we  studied 
the  effects  of  moving  the  trap  level  Ey  away  from  midgap.  Provided  Ey 
remains  relatively  deep  in  the  energy  gap,  i.e.,  =  5kT  above  or  below 
midgap,  the  functional  dependencies  in  our  model  prevail.  That 
IS,   S^^g^^^   as    described    in     (2.33)     for     sufficient    excitation    of    the 
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grain  boundary  depends  on  the  electron  density  in  the  adjacent  grain  and 
hence  manifests  a  nonlinear  boundary  condition  for  the  electron 
continuity  equation   in   the   quasi-neutral    p-region.     When   E-r  is   shall 


ow. 


near  either  the  conduction  or  the  valence  band,  the  carrier  emission 
(exponential)  terms  in  (2.3)  tend  to  dominate  and  render  f,  and  hence  ^g 
nearly  insensitive  to  the  grain-boundary  excitation.  Thus  in  this  case, 
the  nonlinearity  does  not  occur,  and  S^^^^^  is  nearly  constant  equal 
either  to  S^^  (e.g.,  for  E^  near  the  conduction  band)  or  to  lower  values 
(e.g.,    for  Ej  near  the  valence  band). 

However  the  observation  [20,  34,  38]  that  the  neutral  Fermi  level 
of  polysilicon  grain  boundaries  is  near  midgap  implies  that  Ej  is  not 
shallow.  Thus  the  electron  transport  within  the  quasi-neutral  [>region 
is  generally  complicated  by  the  nonlinear  boundary  condition  defined  by 
^n(eff)  ^""  (2.33).  To  enable  the  derivation  of  an  analytic  solution 
for  the  electron  continuity  equation,  we  partitioned  the  [^region  into 
subregions  in  which  we  assumed  the  electron  flow  is  predominantly  one- 
dimensional.  In  the  central  portion  of  the  grain,  we  assumed  that  the 
electron  flow  was  not  significantly  influenced  by  the  grain-boundary 
recombination.  The  one-dimensional  solution  in  this  subregion  then 
provided  a  boundary  condition  fbr  the  coupled  transport  problem  in  the 
subregion  adjacent  to  the  grain  boundary  where  the  strong  influence  of 
the  grain-boundary  recombination,  simulated  by  S^'?  .«  [i.e.,  (2.33)  and 
(2.35)],  was  assumed  to  cause  the  electrons  to  flow  one-dimensionally 
toward  the  grain  boundary. 

The  general    results  of  the   intragrain   analysis,   valid    fbr  r^  >  L 

b         n ' 

was     the     expression     (2.44)      fbr   I^(V) ,  which     includes    an     exp(qV/nkT) 
component  (1  <  n  <  2)    that  occurs  because  of  the  nonlinearity  introduced 
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^^^n(eff)'  '^^^  occurence  of  this  component  complicates  the 
interpretation  of  measured  current-voltage  characteristics  of 
polysilicon  pn  junctions.  Because  the  current  component  deriving  from 
recombination  in  the  junction  space-charge  region  exhibits  the  same 
dependence  on  the  forward  bias  [41],  it  is  generally  impossible  to 
distinguish  between  the  two  components  unless  special  structures  are 
used. 

It  is  likely  however  that  the  junction  space-charge  region 
recombination  current  could  be  predominantly  due  to  recombination 
through  grain-boundary  surface  states  within  the  junction  space-charge 
region  [33].     The  voltage  dependence  of  this  current  I^J^  can  be  derived 

bD 

by  applying  to  this  part  of  the  grain  boundary  those  portions  of  our 
analysis  related  to  (2.10)  and  (2.30).  Generalizing  the  quasi- 
equilibrium  assumption  to  mean  nearly  flat  quasi-Fermi-level  planes  in 
the  space-charge  region  surrounding  the  intersection  of  the  grain 
boundary  with  the  metallurgical  junction  [8,  33],  we  have  for  this 
case  Vg  =  V,   and  hence    from  (2.30) 


tSCR       ,,SCR  I   qV\ 


where  the  reciprocal  slope  factor  is  exactly  two.  Thus,  (2.45)  possibly 
may  facilitate  the  decomposition  of  the  measured  current-voltage 
characteristic  of  a  polysilicon  pn  junction  and  thereby,  with  the 
minority-carrier  transport  model  (2.44)  developed  herein,  provide 
physical  insight  into  the  performance  of  polysilicon  bipolar  devices  and 
their  optimal    designs. 
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From    the    grain-boundary   analysis   described    in    this  chapter,    it   is 
clear    that    the   boundary   condition    at    the    grain    boundary   adjacent   to   a 

quasi-neutral    grain    is,    in   general,    nonlinear.      We   will    experimentally 

PR 
demonstrate    in    the    next   chapter    the    nonlinearity    in   S  ?    ^^v    using    EBIC 

measurements    [21,    22]    interpreted    quantitatively    via    a    computer-aided 

numerical     solution    of    the    underlying    transport    problem.        Numerical 

simulation    of    the    EBIC     is    used    because    of    the    complexity    of    the 

transport     problem,     i.e.,     the     three-dimensional     continuity     equation 

subject    to    the    nonlinear    boundary    conditions.       With    the    aid    of  the 


grain-boundary   model    developed   in    this  chapter,    we   will    obtain    typical 
val  ues     for  \ 
poly  si  1  icon. 


values     for  N^-j.^^^^^   and     Ey    at    a     grain     boundary     in     Wacker     (cast) 


CHAPTER   3 

DETERMINATION   OF   GRAIN-BOUNDARY   RECOMBINATION   VELOCITY 

FROM   ELECTRON-BEAM-INDUCED-CURRENT  MEASUREMENTS 


3.1     Introductic 


The  minority-carrier  transport  model  developed  in  Chapter  2  reveals 
that  the  grain-boundary  recombination  velocity  s''?  ...  is,  in  general, 
dependent  on  the  excitation  level  (carrier  density)  at  the  grain 
boundary.  This      nonlinearity      can      cause      unique     current-voltage 

characteristics  for  polysilicon  diodes.  To  provide  some  experimental 
support  for  the  model,  we  investigate  in  this  chapter  the  recombination 
properties  of  grain  boundaries  using  the  scanning  electron  microscope 
(SEM)  in  the  electron-beam-induced-current  (EBIC)  mode  [42].  The  EBIC 
technique  has  been  widely  used  to  measure  bulk  diffusion  length  [21,  22, 
43-45]  and  surface  properties  [21,  22,  24,  46-49]  in  semiconductor 
devices. 

There  are  nimerous  advantages  in  using  EBIC  over  other  forms  of 
excitation  to  characterize  semiconductor  devices.  The  EBIC  probe  volume 
has  been  well  investigated  and  is  well  defined.  For  example,  the 
penetration  depth  of  an  electron  beam  is  dependent  on  the  atomic  number 
of  the  semiconductor  and  is  independent  of  the  energy  gap  [21].  This  is 
in  contrast  to  optical  excitation  [25],  in  which  the  absorption 
coefficient  depends  strongly  on  the  energy  gap  and  possibly  on  the 
impurity  concentration.  The  EBIC  technique  is  well  suited  to  probe 
finite    semiconductor  regions   since   the   excitation  volume   is   small.      For 
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example,    in    polysilicon    the    EBIC    technique    facilitates    isolation    of  a 

grain    boundary     from    the    adjacent    quasi-neutral     grain    without    tedious 

sample     preparation,     such    as     fabrication    of    small-area    diodes    [16]. 

Furthermore    EBIC  is  a    potentially   powerful    tool     for  measuring   transport 

properties      near      semiconductor      surfaces     [46,      49],      especially      in 

integrated  circuits. 

PR 

In  this  chapter  we  will  determine  S\g^^v  at  a  grain  boundary  from 
an  analysis  of  the  EBIC  response  in  a  grain  adjacent  to  it.  The  grain 
boundary  is  assumed  to  be  perpendicular  to  the  collecting  (n"'"p) 
junction,  and  the  electron  beam  traverses  it  as  illustrated  in 
Fig.  3.1.  The  quantitative  interpretation  of  the  EBIC  response  requires 
the  solution  of  the  underlying  minority-carrier  (electron)  transport 
problem  in  the  p-type  base,  which,  in  general,  is  three-dimensional  and 
has  nonlinear  boundary  conditions.  The  three-dimensionality  arises 
because  the  carriers  are  generated  over  only  a  finite  region  in  the 
semiconductor,    while    the   nonlinearity   in    the   boundary   condition    is  due 

PR 

to  the  dependence  o^S^^^^^^  on  the  excitation  as  shown  in  Chapter  2. 
The  existence  of  a  finite  hoi e-electron- pair  excitation  volume,  which, 
in  general,  is  described  by  a  Gaussian  function  [47,  50,  51],  further 
complicates  the  analysis. 

The  conventional  methods  to  solve  the  transport  problem  treat  the 
grain  boundary  as  a  surface,  which  can  be  justified  using  the  "folding 
technique"  [26].  The  EBIC  response  in  the  vicinity  of  a  surface 
perpendicular  to  the  collecting  junction  was  first  derived  by 
Van  Roosbroeck  [52]  for  the  case  of  a  point  generation  source.  Later 
8erz  and  Kuiken  [25]  developed  an  analytic  model  by  assuming  a  spherical 
generation     source     (volume)     and  using  the    method     of    images.       These 
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analyses     have     been     used     in     the     determination     of    minority-carrier 

lifetime       and       grain-boundary       recombination       velocity       from       EBIC 

measurements  [21,    22,    24,    47].      However   the    point  [52]   or  the   spherical 

[25]   generation- source  models  do  not  give  a  realistic  description  of  the 

actual    electron-beam   generation  [50].      Furthermore   none  of  the   analyses 

[21,    22,    24-26,    47,    52]    account    for    the   nonlinear  S^?    ..>    at   the   grain 

n(eff)  ^ 

boundary,     the    effects    of    which    depend    critically    on    the    generation 
source. 

The  approach  we  take  is  to  solve  nunerically  the  minority-carrier 
transport  problem  underlying  the  EBIC  response  subject  to  the  nonlinear 
boundary  condition  at  the  grain  boundary.  The  three-dimensional, 
steady- state  electron  continuity  equation  is  reduced  to  two  dimensions 
by  recognizing  that  the  grains  are  sufficiently  wide  (in  the 
y-direction)  that  the  y-dependence  of  the  electron  concentration  is 
insignificant  (see  Fig.  3.2).  The  solution  of  the  two-dimensional 
continuity  equation  yields  the  electron  density,  N(x,z).  The  EBIC 
collected  at  the  junction,  I^gic  ^' ^  subsequently  calculated.  Both  the 
spherical  and  the  Gaussian  generation- source  models  have  been  considered 
in  the  numerical    analysis. 

Our  results,  obtained  by  nunerically  solving  the  two-dimensional 
continuity  equation,  are  in  agreement  with  published  analytic  solutions 
[25,  26].  .  Experimental  data  obtained  on  Wacker  (cast)  polysilicon  cells 
demonstrate  the  dependence  of  S^^^^^^^  on  the  excitation  level,  which  is 
consistent  with  the  results  in  Chapter  2.  With  the  aid  of  the  numerical 
analysis,  we  predict  values  for  the  effective  grain-boundary  trap 
density  N^^^^^^^  at  passivated  (see  Chapter  4)  and  unpassivated  grain 
boundaries  in  Wacker     polysilicon.       We  also    find  that     using   the 
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assumption  of  a  simplified  spherical  generation  source  to  predict  the 
EBIC  response  is  valid  provided  the  beam  is  farther  than  2Zq  from  the 
grain  boundary,  where  Zq  is  the  penetrating  range  of  the  primary 
electrons  [21]. 

3.2     Formulation  of  the  Problem 

The     configuration     of    an     assumed     representative     grain     of    the 

polysilicon   n'''p  junction   is   shown   in   Fig.    3.2.      The  grain  is  assumed  to 

extend     infinitely     in     the     positive    x-   and     z-directions    and     in     the 

positive    and    negative    y-directions.        The    metallurgical     junction    is 

located  at  a  distance  x^    from   the  SCR  edge   top   surface,   and  W^^r^  is  the 

width    of    the    (nearly    one-sided)    junction     space-charge    region.       The 

electron  beam  is  incident   from  the  top  surface  and  creates  electron- hole 

pairs  over  a    finite  region   in   the    semiconductor.     The  collection  of  the 

electrons,   which  are  assuned  to    flow  only  by   diffusion,   by  the  junction 

gives    rise    to    I^gjp-       Our    objective    is    to    characterize    I^gjQ,    which 

decreases    as    the    beam    moves    closer    to    the    grain    boundary    because    of 

recombination    losses    there.       From    the    EBIC    response    we    can    determine 

S  ,    ffv    at    the    grain    boundaries    as    well     as    the    electron    diffusion 
niett) 

length  L^   in  the  grains,  which  we  assu^ie  to  be  everywhere  constant. 

In  order  to  interpret  the-  EBIC  measurements  quantitatively,  we 
solve,  numerically,  the  electron  transport  problem  in  the  p-type  base 
subject  to  the  appropriate  boundary  conditions.  The  steady- state 
continuity  equation  for  electrons  when  an  incident  beam  of  electrons 
creates  g(x,y,z)   electron-hole   pairs  (cm"'^-sec'-'-)    is 


V    .  J^   -  qU  +  qg  =  0       .  (3.1) 
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In  (3.1)  J^  is  the  electron  current  density,  q  is  the  electron  charge, 
and  U  is  the  net  thermal  recombination  rate  (cm"-^-sec"M .  For  a 
homogeneously  doped  base  under  low-level -injection  conditions,  (3.1)  can 
be  written  as 


dh       dh   .    d^H         N  ^  g  x,y,z         ^ 
3x^       ay'^       dz^       l^  ^n 

where  N(x,y,z)  is  the  excess  electron  density  at  any  point  in  the  base, 
D^  is  the  electron  diffusion  coefficient,  and  L^^  =  (D^t  )^^^  for  an 
electron  1  i  fetime  Tp. 

Complex  analytical  solutions  of  (3.2)  have  been  published  for  the 
cases  of  point-  [26],  spherical-  [25],  and  Gaussian-  [53]  generation 
sources.  These  solutions  are  not  completely  general  since  they  do  not 
account    for    the    nonlinearity    in    the    boundary    condition    at    the    grain 


GB 
n(eff) 


on 


boundary.  The  nonlinearity  occurs  due  to  the  dependence  of  S 
the  excitation,  i.e.,  N  at  the  grain  boundary,  given  by  (2.33).  The 
incorporation  of  the  nonlinear  boundary  condition  complicates  the 
solution  of  (3.2),  and  hence  we  must  resort  to  a  computer-aided 
numerical  solution.  The  direct  numerical  solution  of  (3.2)  in  three- 
dimensions,   however,  would  involve  excessive  computer  times. 

Fortunately,  we  can  reduce  (3.2)  to  a  two-dimensional  equation  that 
faithfully  represents  the  actual  electron  transport  if  the  grains  are 
sufficiently  wide,  i.e.,  the  surfaces  in  the  y-direction  are  far  away 
(~»).  Ascertaining  that  the  EBIC  response  is  invariant  along  the 
y-direction  (due  to  the  symmetry),  Donolato  [53]  proved  mathematically 
the  validity  of  the  simplification  of  the  transport  problem  from  three- 
to   two   dimensions.      He  concluded   [53]    that   the   detailed   distribution  of 
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the     generation      function     g(x,y,z)      along     the     y-direction     does     not 

influence  the  resulting  value  of  the  EBIC,   and   hence  it  is   sufficient  to 

solve,    as  we  will,    the   two-dimensional    (in   x  and   z)    continuity  equation 

to  calculate   I^bic-     ^®  ^^"^^    show  qualitatively  that   I^gjQ  in   the  actual 

(three-dimensional)    case  differs   from  that  calculated    from  our  nuTierical 

(two-dimensional)     analysis    by    a    constant     factor     for    all     excitation 

conditions.  This      difference      does     not      hinder     our     quantitative 

interpretation    of   the    EBIC   measurements,    e.g.,    evaluation    of  S^?    ^.^  . 

n(e  f  f) 

since  we  do  not  use  absolute  values  of  current  in  our  interpretation. 

The  simplification  to  two  dimensions  (x  and  z)  is  valid  if  the  flux 
of  electrons  in  the  y-direction  is  inconsequential.  This  indeed  is 
generally  true  with  regard  to  relative  value  of  I^bIC  ^^^"  ^^^ 
electron  beam  is  far  away  from  the  grain  boundary,  the  electrons  that 
diffuse  (initially)  in  the  y-direction  have  virtually  the  same 
probability  of  being  collected  by  the  junction,  rather  than  recombining 
at  the  grain  boundary,  as  those  electrons  diffusing  in  the  x-direction 
toward  the  grain  boundary.  For  this  condition,  the  y-dependence  in 
(3.2)  can  be  neglected.  As  the  beam  is  moved  close  to  the  grain 
boundary,  the  nunber  of  electrons  that  diffuse  in  the  y-direction  is 
small  compared  to  the  number  that  diffuse  in  the  x-  (and  z-) 
directions.  This  is  true  because  the  grain-boundary  surface  (and  the 
junction)  is  an  effective  sink  for  the  minority  electrons,  and  hence 
they  tend  to  diffuse  directly  towards  the  grain  boundary  (or  the 
junction).  For  this  condition  also,  the  y-dependence  in  (3.2)  is  not 
important. 

Therefore  integrating   (3.2)   over  the  y-direction,  we  get 


61 


-Ydy+/      ^dy-4-        /Ndy+^      /     g(x,y,z)    =  0 
QX  -00  92  L  -00  n       _oo 


(3.3) 


where,  in  accord  with  the  above  conclusions,  we  have  neglected  the 
y-dependence  of  the  electron  concentration,  which  implies  a  two- 
dimensional  solution  for  I^gj^.  that  differs  from  the  actual  EBIC  by  a 
constant   factor.     We   then    simplify  (3.3)    to 


•y  ''''!f^^^%^-S^^7t^^^i^  =  o  (3.4) 


3^N(x,z)  a^N(x,z)       ^     N(x,z)    ,  G'(x,z) 


3x  J'        3z^  ^      Lj;  -n 


where  Ky  is  an  effective  y-direction  width  that  serves  as  a  normalizing 
constant.  In  (3.4),  N(x,z)  is  the  electron  concentration  obtained  from 
the  two-dimensional    numerical    analysis,   and 


G'(x,z)   =       /     g(x,y,z)    dy  (3.5) 

y=_c, 

is  the  two-dimensional  generation-rate  function.  Dividing  (3.4) 
throughout  by  K       we  obtain 

9^N(x.z)        9^N(x,z)       N(x,z)   ,  G(x,z)       ^  ,_, 

^7  ^  ^^-13^=0       .  (3.6) 

which  is  the  differential  equation  that  we  solve  to  obtain  the  EBIC 
response.  Detailed  expressions  for  G(x,z)  for  spherical  and  Gaussian 
distributions  will   be  given   in  the  next   section. 

We  have  thus  demonstrated  the  validity  of  the  simplification  of  the 
electron  transport  problem  (3.2)  from  three-  to  two-dimensions.  Further 
support    for   our    two-dimensional    model    is   obtained    when    we   compare    the 
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normalized  EBIC  responses  calculated  from  the  three-  [25]  and  the  two- 
dimensional  analysis.  Excellent  agreement  between  the  two  (normalized) 
responses  is  achieved. 

The  boundary  conditions  for  this  problem  are  N(z=0)  =  0  at  the  edge 
of  the  (shorted)  junction  space-charge  region;  N(z  ^.  <»)  =  o  since  L  < 
Wg  where  Wg  is  the  thickness  of  the  base;  N(x  ^-  «)  =0  since  the  grain 
is  semi-infinite  (>>  L^)  in  the  x-direction;  and  at  the  edge  of  the 
grain-boundary   space-charge  region   (x=0), 

"ni'^Meff)"       •  (3-7) 

PR 
Where  S^^^^^^    is  given  by  (2.33).     In   (3.7)   we  have  implicitly  accounted 

for    the    flux    of  electrons   that   do   not  recombine   at   the   grain   boundary 

but    flow   through   it   to    the    adjacent  grain.      Qualitatively   this  can   be 

visualized    by    "folding    over"     the    region    to    the    left    of    the    grain 

boundary    onto    that    to    the    right    [26],    and    noting    that    the    electron 

distributions   in    the    adjacent   grains   add    to   yield    a    solution   of  (3.6) 

provided   1)    the    transport   problem   1s  linear,    2)    L^   is  the    same   in  both 

grains,     and    3)     quasi-equil  ibrium    prevails    across    the    grain-boundary 

space-charge   region.      When    these  conditions  obtain,    the   grain-boundary 

recombination    is    properly    characterized    by    (3.6)     and    (3.7)     provided 

N(x,z)    is   recognized    to   be    the    folded    sim  of  the   left-    and   right-side 

solutions.       This    is    shown    to    be   mathematically    valid    in    Appendix    B. 

From    the    numerical     solution    N(x,z)     of    (3.6)     we    calculate    the    EBIC 

collected  at  the  junction: 


'eBIC  =   Vn  J     '\,,  *       •  <3.3) 
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3.3  The  Volume  Distribution  of  the  Generation  Rate 
Before  we  solve  (3.6),  it  is  important  to  model  the  volume 
distribution  of  the  electron-beam- induced  hoi  e- electron- pair  generation 
rate  G(x,z).  The  simplest  representation  is  a  point-source  excitation 
[21,  26,  52],  i.e.,  the  electron-hole  pairs  are  generated  at  a  point  in 
the  base.  The  distance  of  the  point- source  from  the  junction  is 
determined  by  the  energy  of  the  incident  electron  beam.  It  is  obvious 
that  the  point- source  generation  model  does  not  represent  the  actual 
generation. 

Improved  models  for  the  excitation  volume  have  been  proposed.  A 
spherical  generation  source  [24,  25]  has  been  suggested  as  a  reasonable 
representation      of      the      actual       generation      distribution.  More 

sophisticated  models  such  as  the  radial  Gaussian  [50]  and  the  spherical 
Gaussian  [51]  have  been  proposed  and  have  been  shown  to  be  more 
generally  reliable  than  the  spherical  generation- source  model.  In  our 
analyses  we  consider  two  cases:  1)  a  spherical  generation  source 
[24,  25];  and  2)  a  more  complicated  radial  Gaussian  generation  source 
[50].  The  reason  for  analyzing  the  spherical  generation  source  is  that 
it  is  easier  to  implement  it  numerically  than  the  Gaussian  source,  and, 
as  we  see  later,  the  results  obtained  by  assuming  a  spherical  generation 
source  are  in  agreement  with  those  obtained  by  assuming  a  Gaussian 
source,  provided  the  beam  is  far  away  from  the  grain  boundary.  Any 
uncertainties  resulting  from  assuming  a  spherical  source  model, 
particularly  related  to  the  dependence  of  S^?  .^,  on  the  excitation, 
will    be  removed  by  the  Gaussian    source  model. 
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3.3.1     Spherical    Generation   Source 

The  electron  beam  incident  from  the  top  surface,  is  assumed  to 
generate  electron- hole  pairs  in  a  spherical  volume  [25].  In  our 
ninerical  analysis,  in  which  we  have  reduced  the  transport  problem  to  a 
two-dimensional  differential  equation  (in  the  xz-plane),  the  source  is  a 
circle  as  shown  in  Fig.  3.3.  The  radius  of  the  circle,  rg,  is  much 
smaller  than  the  electron  diffusion  length  in  the  base,  and  the 
generation  rate  at  any  point  inside  the  circle  is  assumed  to  be 
independent  of  position  [25]: 


S  =  V^-^gS  (3.9) 


where  Gq  is  the  total  nimber  of  electron- hole  pairs  per  unit  time 
generated  by  the  electron  beam,  and  K^  is  the  normalizing  constant  in 
(3.4). 

The  center  of  the  circle  is  located  at  a  distance  Zq/3    from  the  top 
surface,  where  [21] 


^0  -  ^  ^0  (3.10) 


is  the  primary  electron  range,  i.e.,  the  penetration  of  the  primary 
electrons  in  the  z-direction.  In  (3.10)  E^  is  the  energy  of  the  beam  in 
kilo-electron-volts  (keV),  and  C  is  a  constant  which  is  dependent  on 
properties  of  the   semicondcutor  [21]: 


C=^^10"^  (3.11) 
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where  p  is  the  (mass)  density  of  the  semiconductor.  For  silicon 
C  =  0.0171  X  10"^,  which  yields  Zq  =  8.6  pm  for  E^  =  35  keV.  The 
physical  reason  for  using  a  circular  generation  source  centered  at  a 
distance  2^/3  from  the  top  surface  will  be  explained  in  the  next 
subsection. 


3.3.2     Gaussian  Generation  Source 

When  high-energy  electrons  penetrate  the  semiconductor,  the 
electron-hole  pairs  are  generated  over  a  characteristic  pear-shaped 
volume  rather  than  over  a  spherical  volume  [50].  The  reason  for  this  is 
that  the  energy  dissipation  of  the  incident  electrons  has  been 
experimentally  shown  to  be  given  by  the  Gruen' s  function  [54]  in  the 
depth  coordinate,  i.e.,  in  the  z-direction.  The  Gruen' s  function  is  a 
third-degree  polynomial  in  z  and  has  a  maximum  at  a  depth  determined  by 
the  energy  of  the  incident  beam.  These  results  have  also  been  verified 
by  Monte  Carlo  simulations  [55].  Hence,  in  our  analysis,  we  will  assume 
that  the  distribution  in  the  z-direction  is  adequately  described  by  the 
Gruen'  s   function. 

In  the  radial  (x-)  direction,  we  will  assume  that  the  carrier 
generation      follows     a     Gaussian     distribution     [50].  The     Gaussian 

distribution  has  been  used  to  reliably  model  carrier  generation  due  to 
X-rays  produced  by  an  electron  beam  [56]  and  successfully  in  electron 
lithography  studies  [50].  The  amplitude  A  and  the  width  a  of  the 
Gaussian  function  will  depend  on  z  and  Zq.  The  generation  rate  an  any 
point  can  then  be  expressed  as  [50] 


A(z,z^) 

G(x,z)  =  ^ ^exp 

Kv2[2T:a  (Z,z^)]'^^ 


U-./  ^ 


2aU,z^] 


(3.12) 
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where  Ky  is  the  normalizing  constant  in  (3.4).  The  expression  for  G  in 
(3.12)  has  been  proposed  by  Donolato  [50].  Expressions  for  A  and  a  will 
be  given  shortly.  Note  in  (3.12)  that  G  is  independent  of  the  y- 
direction  since  the  transport  problem  has  been  reduced  to  two 
dimensions.  The  actual  generation  in  the  y-direction  is  accounted  fbr 
implicitly  by  matching  the  experimentally  measured  and  the  theoretically 
calculated  values  of  the  maximum  EBIC  (at  a  point  well  inside  the  grain) 
to  obtain  the  theoretical  EBIC  response,  Iebic(x). 
A  property  of  (3.12)    is 


x=- 


/       G{x,z)dx  =  A(z,z   )/K  .  (3.13) 


0"   y 


From  (3.10)  and  (3.13)  we  can  infer  that  A(z,Zo)  f'epresents  the  relation 
between  the  penetration  depth  of  the  electron  beam  and  its  energy 
[54].  An  expression  for  A(z,Zq)  has  been  experimentally  obtained  [50, 
54], 


■0  y 


^^^'^0^  =Or^^^"/'o^^      •  (3.14) 


where  A(z/Zq)    is  the  Gruen' s   function  [50,   54]; 


A(z/z  )  =  0.6  +  6.21(z7z^)  -  12.4(z7z  ) 


+  5.69   (z7z„)^  for  0  <  z'    <  1.1  z 

0  0 

(3.15) 
=  0  for  z"   >  1.1  z       ; 

0 
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in  (3.15),  z'  -  z  +  (x.  +  W^^|^).  Since  this  fbrm  for  A  yields  a  unity 
integral  between  0  and  1.1  Zq,  integration  of  (3.15)  in  the  range 
0  <  z'  <  »  correctly  yields  Gq/K^.  In  Fig.  3.4  we  have  plotted  A(z,Zq) 
for  Eq  =  35  keV.  We  see  that  the  maximum  value  of  A(z,Zq)  occurs  at 
Zq/3  from  the  top  surface  (z'  =  Zq/3)  ,  where  z^  =8.6  ym.  This 
observation  qualitatively  justifies  the  earlier  assumption  that  the 
center  of  the  approximate  circular  generation  volume  is  at  z'    =  z  /3. 

The   width  of  the   Gaussian    for   a    finite   value  of  the  beam   diameter 
is  given  by  [50] 


a^(z,z„)  =  0.36  d^  +  0. 


11   (z')"/z. 


(3.16) 


where  d  is  the  beam  diameter.  For  EBIC  measurements,  typically 
d  =  500-1000  A.  The  complete  expression  for  G  adopted  in  our  numerical 
analysis  is  (0  <  z  <   1.1  z  ) 

0 


G„/K 
G(x,z)  =-°^ 


A(Z/Z^) 


0     [2Tra^(z,zJ]^/2   ^^P 


77(77)j 


(3.17) 


where  A(z,Zq)  and  o^{z,z^)  are  described  by  (3.15)  and  (3.16).  In 
Fig.  3.5  we  have  plotted  G  as  a  function  x  for  various  values  of  z. 
These   plots  have  been  drawn    for  x     =  30  pm  and  E.  =  35  keV. 


3.4     Results 

The    computer-aided    numerical     solution    of    (3.6),     subject    to    the 

given    boundary    conditions,    was    obtained    by    using    a     finite-difference 

approximation  technique  [57,   58].     The  grain   is  divided  into   rectangular 

grids  and  the  differential   equation   (3.6)   is  approximated  by  its   finite- 
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Fig.  3.4  The  pre-exponential  coefficient  A(z,z  )  of  the  Gaussian  generation 
distribution  versus  the  depth  coordinate.  The  shaded  area 
indicates  the  total  generation  in  the  emitter  and  the  junction 
space-charge  region.  For  our  cells,  (x.  +  W^^^)  =  0.8  ym,  within 
which  only  about  6%  of  the  total  generation  occurs. 
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Fig.  3.5  The  generation  rate  (assumed  to  be  a  radial  Gaussian)  as  a  function 
of  the  distance  x  from  the  grain  boundary  at  different  depths  z. 
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difference    form    at    each    grid    point.       The    resultant    set   of  algebraic 

equations     is     combined     with     that     at     the     boundaries     and      solved 

simultaneously  using  the  Gauss-Seidel    elimination  method  with   successive 

overrelaxation    [57,     58].        The    EBIC    current     (3.8)     is    calculated    by 

numerically    differentiating    N(x,0)    and    integrating    it  over    the   area   of 

the  junction.     Salient   features  of  the  numerical    algorithn  are  described 

in  Appendix  C,  wherein  we  have  also  listed  the  Fortran  computer  program. 

Typically     EBIC     measurements     are     performed     on     shallow-junction 

devices.       In    such    devices    we    can    assume    that    the    generation    in    the 

quasi-neutral     emitter    and    the    junction     space-charge    region    is    much 

smaller   than   the   total    generation.     Hence   the   EBIC   due   to  generation  of 

carriers    in    the    quasi-neutral     emitter    and    the    junction     space-charge 

region,    which    is   a    constant,    can    be    neglected.       For    example,    in    the 

Wacker     polysilicon     cells     used     in     our     experiments       x.   =0.4pand 

^SCR  ""  °*^  ^^'      '^^^    ^^^"^    °^  ^j    ^^s    inferred    from    sheet    resistivity 

measurements,   while  the  value  of  1^3^,^  was  determined    from   the   zero-bias 

capacitance  measurements  performed  on  the  completed  cell.     From  Fig.   3.4 

we  note  that  the  generation  in  the  emitter  and  the  junction   space-charge 

region    is   about  6%   of  the    total    generation   (in   the  entire  cell),   which 

is  negligible. 

The  calculated  I^gj^  as  a  function  of  the  beam  position  is  plotted 
in  Figs.  3.6  and  3.7  fbr  two  different  values  of  the  background  EBIC 
excitation  level,  i.e.,  maximum  \qiq  produced  when  the  beam  is  far 
removed   from  the  grain  boundary: 


^EBIC  =   ^^  /T     [G(x,z)    -  iil^]    dx   dz  (3.18) 
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where   G  is  given   by   (3.17).      These    plots  have  been  drawn    subject  to   the 

PR 
condition    that  S^^^^^^    is    dependent   on    the    excitation    as    described   by 

(2.33).      Note    that   the    theoretical    EBIC    responses   in   Figs.    3.6   and  3.7 

as  wen    as  the  experimental    data   pertain   to   the    same  grain. 

Superimposed  on   the   theoretical    plots  in   Figs.    3.6   and  3.7   are  the 

measured   EBIC   responses  of  a   typical    grain    in   a    standard   ( unpassivated) 

polysilicon   n'^p  cell,    the    ftbrication  of  which  will    be  described   in  the 

next    chapter.       The     two     parameters     L^  and     N^^^^^^^     are     detemined 

simultaneously    by    fitting    the    experimental    data    at  one   current   level; 

the    fits  at  other  current  levels  provide  a  check  on  the  model,   i.e.,   the 

PR 

nonlinear  S^^^^^^   and  the  constant  1^.     By  comparing  the  theoretical    and 

experimental    responses,    we    infer    that   L_   =  120   pm   and   that  N°  ,    ,^,    = 

ST(eff) 
1  x  10  ^  cm"'^   for  the  grain  measured. 

The  error  bars  shown  in  Figs.  3.6  and  3.7  as  well  as  in  the  other 
figures  to  follow  reflect  the  uncertainty  in  the  measurement  of 
thedi  stance,  x^,  of  the  generation  volume  from  the  center  of  the  grain 
boundary.  This  uncertainty  is  due  to  the  difficulty  in  being  able  to 
locate  the  center  of  the  grain  boundary  in  our  measurements. 

In     Figs.     3.6     and    3.7,     we     have     also     plotted     theoretical     EBIC 

responses     based     on     a     constant  S^^^^^ .  The     assumption     of     a 

constant  S^^^^^    seems     to     be     a     reasonable     one     when     the     beam     is 

sufficiently    far    away    from    the    grain    boundary.      However,    as   is  clear 

from.  Figs.   3.6  and  3.7,   the  EBIC  responses   for  the  constant  S^?    ,,,   case 

n(  eff) 
deviate    from    the    experimental     data    as    the    beam    approaches    the    grain 

boundary.        The     physical     reason     for    this    deviation    is    evident     from 

Fig.   3.8,  where  we   have   plotted  S^^^^^^^     versus  x^    for  the  two  different 

excitation  levels  in  Figs.   3.6  and  3.7.     As  the  beam  nears  the  grain 


73 


1.0 


0.9 


X  o 
<  m 

2_uj 

o 

CO 


0.8 


0.7 


Nonlinear  // 

GB 


/  Constant  S^S 

'  nie 


ffl 


0.61 


10 


15 


20 

Xq  (pm) 


25 


30 


35        40 


Fiq.    3.6 


v?r?n  f  /   "measured  (•)  EBIC  (normalized)  from  a  qrain  in  the 
llVr  l^^  f   ^"  unoassivated  grain  boundary  (G  ^  lo'^/sec)   The 
EBIC  has  been  calculated  for  constant  (--)  aRd  nonlinear  (—) 

^n(eff)- 


74 


1.0 


0.9 


xo 

<m 

o 
ffi 


0.8 


0.7 


0.6 


Nonlinear  S 


10  /sec 


/  Constant  S?,? 
/ 


nleffl 


10 


15 


20 
Xr,  (Mm) 


25 


30 


35 


40 


Fig.  3.7  Calculated  and  measured  (4)  EBIC  (normalized)  from  a  qrain  in  the 
yicin-ity  of  an  unoassivated  grain  boundary  (G  -  lO^Vsec)   Note 

the  significance  of  the  nonlinearity  when  the°electron  beam  is  near 
the  grain  boundary. 


75 


boundary  we   find  that  S^^^^^^   decreases  due  to  the  increase  in 

excitation,     i.e.,     N(0,z),     at    the    grain    boundary.        Consequently    the 

actual     EBIC    collected    by    the    junction    at    any    beam    position    is    higher 

than    that    predicted    by    the    unrealistic  constant-S^^^^^,   curves.      This 

demonstrates   the    significant   nonlinearity  in  S^^^^^^    in   accord  with  our 

theoretical    results  in  Chapter  2. 

We    also    note    from    Figs.    3.6    and    3.7,    which    pertain    to    the    same 

grain,    that    the   value   of  S^^^^^^    inferred  by   comparing   the   theoretical 

PR 
(constant  S^^^^^^)     and  experimental    responses,   is  not  unique;   different 

values    are    implied    at    the    same    grain    boundary    at    different    current 
-GB 
'n(eff) 


PR 

levels     when  S^^^^^   is     assumed     to     be     a     constant.         This  further 

demonstrates  the  nonlinearity  and   shows  that  a  constant-s'^?    ..,  model    is 

physically  unreasonable.     We  note  that  when   the  nonlinearity  in  S^? 

^  n(eff) 

is    accounted    for,    comparison    of   theoretical     and    experimental  results 

ST(eff)'   ^^^     nonlinear     model  is     thus 


yields     a     unique     value     for  N^^^^^^^;   the     nonlinear    model      is     thu 
physically  reasonable. 

PR 
The    nonlinearity    in  S^^^^^^    is    further    confirmed    when    we    compare 

the    experimental    data    measured    at   the    two    different   excitation    levels 

(Figs.    3.6    and    3.7).       The    increased    excitation    (Gq)    was    obtained    by 

increasing   the  beam  current  and  keeping   the  beam  energy  constant.     (The 

increase      in      the     beam      diameter     with     excitation     level      [42]      is 

inconsequential     since    the    second    term    on    the    right    side    of  (3.16)    is 

always    dominant.)       If  S^^^^^^   were   a   constant   then    the    percentage    EBIC 

collected  by   the  junction   at  any   given   x^    should  be   independent  of  the 

excitation   level    [25,    47].      However   the   experimental    data    in   Figs.    3.6 

and  3.7    show  a   strong   dependence  on   the   excitation  level.     For  the   same 

value    of   Xq    we    observe    that    the    percentage    EBIC    lost    at    the    grain 
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boundary  at  the  lower  excitation  level  (Fig.  3.6)  is  higher  than  that  at 
the  higher  excitation  level  (Fig.  3.7).  The  reason  for  this  is  obvious 
from  Fig.  3.8  wherein  we  have  plotted  S^^^^^  as  a  function  of  x^  at  the 
two  excitation  levels  for  the  grain  boundary  in  the  standard  cell.  From 
Fig.  3.8  we  note  that  S^^^^^^  at  the  grain  boundary  in  the  standard 
cell  is  lowered  by  about  an  order  of  magnitude  at  the  higher  excitation 
level.  Hence  we  conclude  that  the  nonlinearity  in  S^^^^^  influences 
the  EBIC  response  when  the  beam  is  close  to  the  grain  boundary  and 
should  be  accounted  for  to  ensure  accurate  quantitative  interpretation 
of  EBIC  measurements. 

In  Figs.  3.9  and  3.10  we  compare  the  theoretical  plots  with  the 
experimental  data  obtained  in  a  grain  of  an  al  uminum-passivated  (active) 
n"^p  cell.  The  passivation  technique  will  be  described  in  the  next 
chapter.  By  comparing  the  experimental  data  with  the  theoretical  plots, 
which  were  obtained  by  assuming  G  to  be  a  Gaussian  given  by  (3.17),  we 
infer  that  N^^^^^^^  =  10  cm"  at  a  passivated  grain  boundary.  Once 
again  we  observe  that  S^^^^^  is  nonlinear.  However  the  nonlinearity  is 
less  pronounced  than  at  the  unpassivated  grain  boundaries  in  a  standard 
cell,  particularly  at  the  low  excitation  level.  An  explanation  for  this 
can    be    provided    using    the    grain    boundary    analysis   in    Chapter    2.       For 

^°^  ^ST(eff)'   *B  =  *B0'  "^^^'^    ^'"^PT^'es    that  S^^^^^^         (.  C.N^^^^^^)   is 
virtually  constant. 

The   theoretical    EBIC   response   obtained    in   a   grain    in   the   vicinity 

of  an    unpassivated  grain   boundary   is  compared   to   the   experimental    data 

in        Fig.   3.11.  We       infer       that       at       this       grain       boundary, 

^ST(eff)   "  ^0      cm"   ,  which     suggests    that    this     particular    grain 

boundary  is  a  high-angle  boundary.     This  value  of  N^^^^^^^      is  unusually 
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Fig.  3.8  The  grain-boundary  effective  recombination  velocity  as  a  function 
Of  Xq  at  unpassivated  and  passivated  grain  boundaries.  These 
plots  show  how  S^^g^^^  changes  as  the  electron  beam  moves  toward 
(excites)  the  grain  boundary. 
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high,  possibly  because  of  our  assumption  thatN^^,^^^.  i  s  monoenergetic 
situated  at  midgap.  The  large  value  ofN^y/g^^N  causes  ^g  to  be 
high.  Consequently  the  value  ofS^^^^^^^  does  not  vary  significantly 
with  the  excitation  level  and  is  at  the  kinetic  limit  as  described  in 
Chapter  2. 

In  Fig.  3.12  we  have  plotted  the  theoretical  EBIC  response  obtained 
by  assuming  a  circular  generation  source.  Comparing  it  to  the  Gaussian- 
generation- source  case,  we  observe  that  the  circular-generation 
assumption  is  good  provided  we  interpret  data  far  away  (~  2Zq)  from  the 
grain  boundary.  As  the  beam  approaches  the  grain  boundary  the  predicted 
EBIC  responses  in  the  two  cases  differ  from  each  other.  If  we  interpret 
the  experimental  data  in  Fig.  3.6  using  the  circular  generation  source 
assumption  we  find  that  N^^^^^^^  =  2  x  10^^  cm"^,  which  is  a  factor  of 
two   higher  than  that  obtained  by  assuming  a  Gaussian  generation   source. 


3.5     Discussion 


In    this    chapter    we    determined  sj^^^^   at    typical     passivated    and 

unpassivated    grain    boundaries    in    Wacker    (cast)     polysilicon    from    EBIC 

measurements.  We     demonstrated     the     dependence     of  s''?    ^^,   on     the 

n{  eff) 

excitation    level     as    predicted    by    the     transport    model     developed    in 
Chapter  2. 

In  order  to  facilitate  quantitative  interpretation  of  EBIC 
measurements  we  solved,  using  computer-aided  numerical  analysis,  the 
underlying  transport  problem.  The  general  transport  problem  in  three 
dimensions  was  reduced  to  one  in  two-dimensions  by  neglecting  the 
variations  in  the  (y-)  direction  normal  to  the  plane  defined  by  the 
moving  electron  beam.     This  assumption  is  justified   since  the   surfaces 
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in  the  y-direction  are  far  removed  from  the  excitation  volume.  We  found 
that  the  discrepancy  between  the  actual  solution  and  ours  is  a  constant 
factor  in  I^gjQ  provided  the  y-dependence  in  (3.2)  can  be  neglected, 
which  was  justified  qualitatively.  Our  quantitative  interpretation 
ofS  ,  -:«  is  not  hindered  due  to  this  constant  difference  since  we 
compare  only  relative  values  of  EBIC  current.  The  results  of  the  two- 
dimensional  analysis  agree  well  with  published  analytic  solutions  [25, 
26,  52].  The  simplification  of  the  transport  problem  from  three-  to 
two-dimensions  has  also  been   shown  to  be  mathematically  valid  [53]. 

In  our  analyses,  we  have  incorporated  two  different  models  fbr  the 
hoi e-electron- pair  excitation  volume:  1)  the  spherical  generation- 
source  model  [24,  25];  and  2)  the  Gaussian  generation- source  model 
[50].  The  utility  of  the  spherical- source  model  results  from  its  ease 
of  implementation  in  the  numerical  analysis.  The  EBIC  response  obtained 
by  using  this  model  far  away  from  the  grain  boundary  agrees  well  with 
that  obtained  by  using  the  sophisticated  and  hence  more  complicated 
Gaussian- source   model.       Thus   the    spherical- source  model    is    useful     for 

interpreting   EBIC   data    far  away    from   the  grain  boundary  when  S  ,    o«    is 

PR 
a  constant  (e.g.,  when  S  -    ^x\   =  Sj,.    or  at  passivated  grain  boundaries). 

In  the  spherical- source  model  the  incident  electron  beam  generates 
hole-electron  pairs  in  the  base  in  a  sphere  (circle  in  the  two- 
dimensional  analysis)  whose  radius  is  much  smaller  than  the  electron 
diffusion  length  and  within  which  the  generation  rate  is  a  constant 
[25].  To  remove  any  uncertainties  associated  with  the  spherical -source 
assumption  [47,  50,  51],  we  also  modeled  the  excitation  volume  using  a 
radial  Gaussian  distribution  [50],  i.e.,  at  any  depth  z,  G  is  given  by  a 
Gaussian    function.      In   the  z-direction   the   distribution   is  given  by  the 
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Gruen's  function  [50,  54].  This  description  of  the  generation  source  is 
more  generalized  and  better  represents  the  actual  generation  than  the 
simpl  i  fied   spherical -source  model  . 

The  analytic  solution  of  the  two-dimensional  continuity  equation 
(3.6)  in  the  general  case  of  nonlinear  S^J^^^  and  a  Gaussian  generation 
source  is  formidable.  Hence  we  resorted  to  a  computer-aided  numerical 
solution.  The  solution  was  obtained  by  approximating  (3.6)  using  finite 
diffisrences  [57,  58].  The  resultant  set  of  algebraic  equations  was 
solved  by  the  Gauss- Seidel  elimination  method  [57,  58].  The  EBIC 
current  collected  at  the  junction  was  calculated  as  a  function  of  x^, 
S^(g^^^,  andN^y^g^^.  In  the  analyses  we  have  treated  the  grain 
boundary  as  a  surface,  which  can  be  justified  by  using  the  folding 
technique  [26]. 

The  EBIC  response  of  (al  uminum-)  passivated  and  unpassivated  Wacker 
(cast)  polysilicon  cells  were  measured  experimentally.  The  cell 
fabrication  and  the  passivation  technique  will  be  described  in  the  next 
chapter.  When  the  beam  is  situated  close  to  the  grain  boundary  we 
observed  that  I^gj^  measured  experimentally  was  higher  than  predicted  by 
the  theoretical  plots  for  constant  S^^^^^^.  This  clearly  demonstrates, 
subject  to  our  assumption  of  constant  1^  which  was  indeed  supported  by 
our  results,  that  the  dependence  of  S^^^^^^  on  the  excitation  level  as 
predicted  in  Chapter  2.  We  also  observed  that  the  current  loss  at  the 
grain  boundary  for  a  given  x^  decreased  with  increasing  excitation 
level.  This      provides      further      evidence      for      the      nonl  inearity. 

Physically,    N(0,z)    increases  as  the  beam   approaches   the   grain   boundary 
or     as    the     incident    excitation     level     is    increased.        Consequently, 


^n(eff)'  ^^"^^^    ^'^     proportional     to    [N(0,Z)]l/2^     decreases    and     hence 
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PR 
^EBIC     "' ^     higher     than     predicted     by     constant       S  , ^^   curves     (see 

Figs.    3.6,    3.7,    3.9,    3.10). 

By    comparing    the    experimental    data    with    the    theoretical    plots   we 

D  12       -2 

inferred     that  ^^jjfgff)   ==  10       cm       at     a     typical      unpassivated     grain 

boundary    and  N3J(gf*^   =   10       cm"     at   a    passivated   grain    boundary.      The 

comparison    of    the    theoretical     and    experimental     EBIC    responses    also 

yielded   values    for   the   diffusion   lengths   in   the   grains    simultaneously. 

For  the  grain  of  the  unpassivated  cell,  L     =  120  um,  while   for  the  grain 

of   the    passivated    cell,  L^  =  110  um.      The    determinations    o^Ncj/pff) 

were  based  on  the  assunption  that  the  trap  level,   Ej,   is  at  midgap  (Ey  = 

E^-).       If  Ej   were   below    E^-    [34],    the   value   o^  Ncj(eff)   would  be    higher 

than    that    predicted    by    our    ninerical     analysis.        The    dependence    of 


.GB 
'n(eff) 


S„/„rin   on    the    excitation    level,     however,    remains    unaltered    and    is 


given   by    (2.33),    provided   Ej   is   located    relatively    deep  in   the   energy 
gap,   i.e.,  =  5  kT  below  E^-   (see  Chapter  2). 

Our  analyses  predict  that  the  spherical  generation- source 
assumption  adequately  describes  the  EBIC  response  provided  x  >  2z  . 
For  Xq  <  2Zq,  the  response  for  a  spherical  source  differs  from  that  fbr 
a  Gaussian  source  and  could  cause  at  least  a  factor  of  two  difference  in 
the  determination  o^'^sT(eff)'  ^'^^  reason  fbr  the  discrepancy  is  that 
the  spherical  source  cannot  accurately  describe  the  nonlinearity 
in  S^.  ^^,  which  becomes  important  as  the  beam  approaches  the  grain 
boundary.  For  x  >  2  z  ,  the  influence  of  the  grain  boundaries  on  the 
EBIC  response  is  minimal  and  hence  the  response  is  independent  of  the 
choice  of  the  generation  source.  The  use  of  the  spherical- source  model, 
which  is  easy  to  implement  and  is  widely  used  in  most  analyses  [24,  25, 
47],    is    recommended    to    interpret    EBIC    data     far    away     from    the    grain 
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boundary.     The  model   yields  reliable  results  when  S^^    .^    is  a  constant, 

e.g.,   at  a   passivated  grain  boundary  or  when  S^    ^,,    =  S^,  . 

n(  e  t  f)         KL 

In     Chapter     2     we     studied     grain-boundary     recombination     and     its 

influence  on   the  base   recombination   current  of  a    polysilicon   diode.      In 

this     chapter     we     described     an     experimental      technique     to     determine 
PR 
^n(eff)'  '^'^^^'^  ^'s  3  measure  of  the  grain-boundary  recombination.      It  is 

obvious  that  grain-boundary  recombination  deteriorates  the  performance 
of  bipolar  polysilicon  devices.  For  example  in  a  polysilicon  solar 
cell,  the  grain-boundary  recombination  both  in  the  junction  space-charge 
region  as  well  as  adjacent  to  the  quasi-neutral  base  can  limit  the  cell 
efficiency  [8].  In  order  to  improve  polysilicon  device  performance  it 
is  necessary  to  suppress  grain-boundary  recombination.  Such  suppression 
will  also  likely  reduce  carrier  scattering  by  grain  boundaries  which 
lowers  <j)g  at  the  grain  boundaries.  This  results  in  an  improvement  of 
the  carrier  mobility  in  polysilicon  [5],  for  example  that  used  in 
integrated  circuits  [3]. 

In  the  next  chapter  we  suggest  an  experimental  method  to  reduce 
grain-boundary  recombination  by  passivating  grain  boundaries  using 
preferential  diffusion  of  aluminum.  The  passivation  is  assessed  by 
comparing  the  EBIC  response  in  the  vicinity  of  an  unpassivated  and  a 
passivated  grain  boundary.  The  value  of  S^^^^^^  is  determined  from  the 
EBIC  measurements  with  the  aid  of  the  nimerical  solution  of  the 
transport  problem  described  in  this  chapter.  Possible  additional 
benefits  of  aluminun,   such  as  gettering,  are  also  examined. 


CHAPTER   4 

POTENTIAL    IMPROVEMENT  OF   POLYSILICON   SOLAR  CELLS 

iY  GRAIN-BOUNDARY  AND   INTRAGRAIN   DIFFUSION  OF   ALUMINUM 


4.1  Introduction 
The  conversion  eflnciency  of  polycrystall  ine  silicon  (polysil  icon) 
solar  cells,  as  pointed  out  in  Chapter  2,  is  limited  by  high  defect 
densities  within  the  grains  as  well  as  at  the  grain  boundaries.  On  the 
one  hand,  the  efficiency  of  bulk  polysilicon  cells  having  relatively 
large  grains  is  generally  limited  by  high  densities  of  intragrain 
defects  that  cause  short  minority-carrier  lifetimes;  the  grain-boundary 
recombination  is  hence  rendered  virtually  unimportant.  The  efficiency 
of  thin- film  polysilicon  cells  having  relatively  small  grains  may,  on 
the  other  hand,  be  generally  limited  by  high  densities  of  grain-boundary 
traps  that  produce  high  recombination  velocities.  It  is  therefore  of 
interest  to  consider  processes  that  could  possibly  passivate  the  grain 
boundaries,  i.e.,  reduce  the  grain-boundary  trap  densities,  and  also 
getter  the  grains,  i.e.,  remove  defects  (impurities)  from  the  active 
portions  of  the  grains.  Such  a  process  is  perhaps  the  diffusion  of  a 
specific  impurity  whose  grain-boundary  diffusion  coefficient  is  higher 
than  its  bulk  coefficient,  which  implies  preferential  grain-boundary 
diffusion  [59],  and  whose  presence  at  the  silicon  surface  increases  the 
solubility  of  common  defects,  which  implies  bulk-to- surface  transfer  of 
defects. 
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Phosphorous,  which  is  an  effective  gettering  agent  in  silicon  [60], 
also  has  the  potential  to  be  an  effective  grain-boundary  passivating 
agent      [16].  However      preferential       grain-boundary      diffusion      of 

phosphorous  much  deeper  than  ~  10  ym  is  not   practical    because  it  must  be 
done    at   a    temperature    (>    1000°C)    considerably    higher    than    the    normal 
solar  cell   junction- formation   temperature   (<   900°C) .      It   has  been    shown 
that    lithitn    tends    to    suppress    grain-boundary    recombination,    but    its 
stability    and    effectiveness    are    dubious   [17].      Monatomic    hydrogen    has 
been    shown    to    be    extremely    effective    in    passivating    grain    boundaries 
[18],     but     the     necessity    of    inordinate     sample     preparation     and    the 
unorthodox    nature    of  the    hydrogenation    process  make    it   impractical    in 
solar  cell    fabrication.      Passivation   of  grain   boundaries   has  also   been 
achieved      using      deuterium      [19]      but,      as      is      hydrogenation,      the 
incorporation    of   deuteriin    into    the    polysilicon    is    incompatible    with 
conventional    solar  cell    processing.     Hence   it  would  be  useful    to  design 
a      process     that     1)      is     compatible     with     conventional      solar     cell 
fabrication;   2)    passivates  most  of  the  grain  boundaries  in   a  relatively 
short  time;   and  3)    produces  concomitant  gettering  of  intragrain  defects. 
There    are    indications    that    aluminum,    an    acceptor- type    dopant    in 
silicon,    could   be    an    ideal    agent    for    improving    polysilicon    solar   cell 
efficiency.        Evidence    of    anomalous     diffusion    of    and    gettering    by 
aluminin    in     silicon     dates    back     to    the     initial     development    of    the 
alrnii  nun- alloyed     n''pp+     back- sur  face- fi  el  d     (BSF)      silicon     solar     cell 
[61].  More      recently      unusually       high      grain-boundary      diffusion 

coefficients  for  aluninum  in  polysilicon  have  been  measured  [62],  and 
unexpected  improvements  in  the  spectral  response  of  "bicrystall  ine" 
silicon  cells  having  aluminum  "BSFs"    have  been  observed  [63]. 
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The  passivation  of  a  grain  boundary  by  preferential  aluminum 
diffusion  might  occur  via  1)  reduction  of  the  grain-boundary  (deep- 
level)      trap     density      resulting      from      removal      of     dangling      bonds; 

2)  fionnation  of  a  high-low  junction  at  the  grain  boundary,  which  reduces 
the    effective    recombination    velocity  S^^^^^      (see    Chapter    2);    and/or 

3)  formation  of  a  pn  junction  at  the  grain  boundary,  which  effectively 
increases  the  collecting  junction  area.  It  is  possible  then  that 
aluminun  diffusion,  which  is  compatible  with  solar  cell  fabrication, 
could  improve  polysilicon  solar  cells  by  passivating  grain  boundaries  as 
well  as  by  gettering  intragrain  impurities  and  by  creating  an  effective 
BSF. 

In  this  chapter  we  present  results  of  experiments  designed  to  study 
the  effects  of  alumintm  diffusion  on  the  efficiency  of  both  bulk  and 
thin- film  polysilicon  solar  cells.  In  the  case  of  bulk  polysilicon, 
which  possibly  may  soon  be  a  cheap  substitute  for  silicon  in  solar  cell 
technology,  we  focus  our  attention  on  the  conventional  Al-BSF  process 
[61,  64],  which  should  ultimately  yield  the  highest  attainable 
efficiencies  as  it  has  in  silicon  cells  [65].  Possible  additional 
benefits  due  to  the  back  aluminun  diffusion  (alloying),  e.g.,  gettering 
and  grain-boundary  passivation,  are  investigated.  With  regard  to  thin- 
film  polysilicon,  we  examine  gettering  and  grain-boundary  passivation 
resulting  from  aluminum  diffusion  at  low  temperatures  (below  the 
eutectic- formation  temperature,  577°C)  from  the  front  surface.  Such 
low- temperature  diffusion  has  been  shown  to  yield  preferential  aluminum 
spikes  down  dislocations  [66]. 

To  assess  the  grain-boundary   passivation,  we  use  the  electron-beam- 
induced-current     (EBIC)     technique.         A     qualitative     analysis    is    done 
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inUially  to  detect  grain-boundary  passivation.  Quantitative  estimates 
of  grain-boundary  parameters,  e.g..  S^J^^^^,  based  on  EBIC  data  are 
derived  using  the  computer-aided  ninerical  simulation,  details  of  which 
were  disclosed  in  Chapter  3,  of  the  minority-carrier  transport 
underlying  the  EBIC  response  of  a  pn  junction  in  the  vicinity  of  a  grain 
boundary.  Current-voltage  characteristics  and  reverse-bias  capacitance 
measurements,  which  yield  the  base  doping  density  of  the  cell,  are  used 
to  determine  the  minority-carrier  diffusion  length  in  the  base,  which 
reflects  the  benefit  produced  by  the  Al-BSF   process. 

Our  experimental  results  show  that  the  conventional  Al-BSF  process 
results  in  longer  diffusion  lengths  because  of  (bulk-to- surface) 
gettering  of  intragrain  impurities,  irrespective  of  indications  that  the 
grain  boundaries  emit  impurities  during  the  high- temperature  process. 
Our  results  also  suggest  that  the  minority-carrier  lifetime  is  longer 
when  measured  using  light-bias  than  when  measured  in  the  dark.  This 
suggestion  is  consistent  with  the  theory  [67]  that  defect  clusters  in 
polysilicon  are  effectively  negated  by  high  carrier-injection  levels. 
No  grain-boundary  passivation  near  the  front  junction  is  apparent  after 
the  BSF  process.  However  we  find  significant  passivation  near  the  front 
junction  resulting  from  low- temperature  aluminim  diffusion  from  the 
front  surface,  but  no  observable  gettering.  Consequently,  it  is 
possible  that  aluminum  diffusion,  when  optimally  done,  could 
substantially  increase  the  efficiencies  of  both  bulk  and  thin-film 
polysilicon    solar  cells. 
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4.2     Cell    Fabrication 


The  cells  used  in  this  study  were  fabricated  on  f>-type  Wacker 
(cast)  polysilicon  substrates.  After  standard  cleaning  and  chemical 
polishing  of  the  front  surface,  the  substrate  thickness,  Wg,  is  about 
350  ym.  The  doping  density,  N^^,  varies  from  SxlO^^cm"^  to  IxlO^^cm"^, 
as  determined  by  reverse-bias  capacitance  measurements  on  the  completed 
cells. 

Three  different  kinds  of  cells  were  fabricated  for  comparison  in 
both  the  bulk  and  thin- film  polysilicon  studies.  Active  cells  were 
subjected  to  aluminum  diffusion  while  control  cell  s  were  not.  However, 
the  control  cells  were  subjected  to  the  same  thermal  stress  as  were 
their  aluminum-di  f  fused  counterparts.  Several  standard  cells,  which 
were  neither  al  umin urn- di  f fused  nor  subjected  to  any  abnormal  thermal 
stress,  were  also    fabricated. 

4.2.1     Bulk   Process 


For  the  studies  applicable  to  the  development  of  bulk  polysilicon 
solar  cells,  we  used  the  conventional  Al-BSF  process  for  silicon  [64]  to 
fabricate  the  active  cells.  About  1.5  m  of  aluminum  was  evaporated 
onto  the  back  surface  of  the  cell  after  the  front  phosphorous  (n^) 
diffusion.  The  aluminum  was  then  diffused  into  the  active  cell  during 
the  alloying  process  at  8QQ°C  for  1  hour  in  forming  gas.  The  control 
cell  underwent  an  identical  heat  treatment  but  without  aluminum.  Both 
large-area  (0.4  -  0.45  cm^)  and  small-area  (30-mil -diameter)  cells  were 
fabricated. 

Fabricating  aiiall-area  cells  enables  control  of  the  grain-boundary 
density   in   each  cell.      The    surface   and  edge  leakage  currents,   which  are 
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predominant  in  large-area  cells,  can  be  effectively  controlled  by  means 
of  an  MOS  guard-ring  overlapping  the  edge  of  the  emitter  in  the  small- 
area  cells  [16].  In  addition,  by  comparing  the  diffusion  lengths, 
determined  from  the  current-voltage  characteristics,  of  the  large-area 
(measured  under  light)  and  the  gnall-area  (measured  in  dark)  cells,  we 
can  deduce  the  effect  of  the  light-bias  on  minority-carrier  lifetime. 

The  n  p  junction  in  the  large-area  cells  was  formed  by  diffusing  in 
phosphorous  at  875°C  for  30  minutes  using  POCI3  as  the  source.  The 
resulting  junction  depth  is  about  0.4  pm.  Subsequent  to  the  aluminum 
diffusion,  which  followed  the  phosphorous  (n^)  diffusion,  about  0.5  un 
of  aluminum  was  evaporated  onto  the  back  surface  to  serve  as  a 
contact.  The  control  and  the  standard  cells  were  then  sintered  at  450°C 
for  15  minutes  in  forming  gas.  Following  a  mesa  etch  that  isolated  the 
cells,  the  front  contact  was  formed  by  evaporating  a  300-A-thick 
aluminum  "dot".  Details  of  these  processing  steps  as  well  as  those  of 
the    small-area  and  the  thin- film  cells  are  given   in  Appendix  D. 

The  junction  in  the  small-area  cells  was  formed  by  a  phosphorous 
predeposition  at  950°C  for  40  minutes  using  P0C13  as  the  source  followed 
by  a  drive-in  diffusion  at  1000°C  for  27  minutes.  This  resulted  in  a 
junction  depth  of  approximately  0.7  ym.  Aluminun  was  then  diffused  into 
the  active  cells  from  the  back  surface.  The  front-  and  back-contact 
aluminum  evaporations,  5000-A  and  6000-A-thick  respectively,  were 
performed  subsequently.  As  in  the  case  of  the  large-area  cells,  the 
control  and  the  standard  cells  were  sintered  at  450°C  for  15  minutes 
after  the  back-contact  evaporation. 

In  order  to  suppress  surface  and  edge  leakage  currents  in  the 
small-area  cells,  a   silicon-dioxide  layer   (~  1700-A-thick)   was  thermally 
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grown  on  the  top  surface  during  the  drive-in  step,  and  a  metal  gate  was 
deposited  overlapping  the  edge  of  the  diffused  region  [16].  A  negative 
gate  bias  (-20V  to  ensure  surface  accumulation  in  the  base)  was  used 
while  the  current-voltage  characteristic  of  the  cell    was  being  measured. 

4.2.2     Thin-Film  Process 

For  the  studies  applicable  to  the  development  of  thin- film 
polysilicon  solar  cells,  the  n^p  junction  was  formed  by  diffusing 
phosphorus  at  875°C  for  30  minutes  using  POCI3  as  the  source.  The 
resulting  junction  depth  is  approximately  0.4  pm.  A  mesa  etch  was  used 
to  isolate  the  cells.  The  active  and  control  cells  were  fabricated  on 
sequential    wafers   from  the  Wacker  ingot. 

We  use  bulk  (large-grain)  polysilicon  to  imply  the  effects  of 
aluminim  diffusion  in  thin- film  (small-grain)  cells.  This  material 
facilitates  electrical  isolation  of  a  grain  boundary  (using  EBIC)  and 
enables  direct  quantitative  assessment  of  the  passivation  yielded  by  the 
aluminum  diffusion.  We  anticipate  that  the  low- temperature  diffusion 
will  similarly  passivate  grain  boundaries  in  thin- film  polysilicon. 
However  because  the  grain  boundaries  may  differ  significantly  from  those 
in  Wacker  polysilicon,  the  extent  of  the  passivation  might  be  different. 
We  diffused  aluminum  into  the  active  cells  from  the  front  surface 
at  temperatures  below  the  Si-Al  eutectic  temperature  (577°C)  to  avoid 
irreparable  damage  to  the  junction  region.  IC-grade  aluminum,  1.5-pm- 
thick.  was  evaporated  onto  the  front  surface  before  the  phosphorous  (n^) 
diffusion.  The  aluminum  diffusion  was  then  done  at  450°C  for  3-6  hours 
in  forming  gas.  Following  the  removal  of  excess  aluminum,  a  thin  layer 
of  silicon   was  etched  off  to   avoid    shunting   of  the   n+p  junction  that  is 
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subsequently  formed.  The  control  cells,  sequential  to  the  active  cells, 
were  also  etched  in  addition  to  undergoing  an  identical  heat 
treatment.  After  the  subsequent  phosphorous  (n"*")  diffusion,  about 
6000  A  of  aluminun  was  evaporated  onto  the  lapped  back  of  the  cell  to 
serve  as  the  back  contact.  The  front  contact  was  formed  by  evaporating 
a  300-A-thick  aluminim  "dot"   after  the  isolating  mesa  etch. 

The  standard  cells  did  not  undergo  any  thermal  treatment  except  for 
the  phosphorous  (n"^)  diffusion  and  a  sinter  for  15  minutes  at  450°C 
following  the  back-contact  evaporation. 

4.3.     Results 
4.3.1     Bulk  Process 


Significant  effects  of  the  Al-BSF  process  are  revealed  in  Tables 
4.1  and  4.2,  where  we  give  representative  values  of  the  electron 
diffusion  lengths  L^  (averaged  over  each  cell)  in,  respectively,  the 
small -area  and  the  large-area  active,  control,  and  standard  cells.  The 
values  of  L^  were  derived  from  the  quasi-neutral-base  recombination 
current,  which  was  experimentally  determined  by  measuring  the  dark 
current-voltage  characteristics  [68].  We  have  assumed  that  the  emitter 
recombination  current  in  both  the  large-  and  gnall-area  cells  is 
negligible.  The  large  values  of  the  base  recombination  current 
determined  from  the  current-voltage  characteristics  justify  this 
assumption  [68]. 

The  forward  current-voltage  characteristics  of  the  snail -area  cells 
were  obtained  by  measuring  the  dark  current,  I^,  through  the  cell  for 
different  values  of  the  forward  bias,  Vq,  across  the  cell.  This 
measurement  as  well   as  the  extrapolation  of  the  Iq-Vq  characteristics  to 
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Table  4.1         Measured    electron    diffusion    lengths    in    bulk-cell -related 
small-area  active,  control,   and   standard  n'^'p  cells. 


Device 

N;^A  (cm-3) 

Lp  (m) 

STANDARD 

6  X  10^5 

74 

6  X  10^5 

33 

7  X  10^5 

27 

7  X  10^^ 

55 

7  X  10^5 

49 

7  X  10^5 

47 

7  X  10^5 

57 

7  X  10^5 

74 

7  X  10^^ 

63 

ACTIVE 

8  X  10^5 

168 

7  X  10^5 

276 

9  X  10^5 

173 

9  X  10^5 

92 

7  X  10^5 

146 

1  X  10^^ 

92 

6  X  10^5 

244 

CONTROL 

9  X  10^5 

20 

1  X  10^6 

20 

1  X  10l6 

41 

8  X  10^5 

14 

7  X  10^5 

34 

9  X  10^5 

24 
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Table  4.2  Measured  electron  diffu^sion  lengths  in  large-area  active, 
control,  and  standard  n"'"p  cells.  Each  row  refers  to  cells 
fabricated  on   sequential    wafers. 


DEVICE 

^AA  (cm-3) 

L,  (pm) 

STANDARD 

5  X  10^5 

91 

5  X  10^5 

96 

ACTIVE 

5  X  10^5 

160 

5  X  10^5 

130 

4  X  10^^ 

122 

CONTROL 

5  X  10^5 

34 

6  X  10^5  . 

57 
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obtain  the  base  recombination  current  [68]  was  achieved  using  an  HP 
9845A  desktop  computer.  The  cells  were  maintained  at  a  constant 
temperature,  T  =  25°C,  during  the  measurement  by  placing  them  in  a 
Stanam  oven.  In  Fig.  4.1  a  typical  Iq-Vq  characteristic  obtained  on  a 
small -area  active  cell    is  shown. 

The  base  recombination  current  of  the  large-area  cells  was 
experimentally  determined  by  measuring  the  short-circuit  current,  Icp, 
versus  the  open-circuit  voltage,  Vg^,  characteristic  resulting  from 
different  illumination  levels  [68].  The  reason  for  using  the  Isc-Vqc 
characteristic  to  determine  the  base  recombination  current  is  that  the 
series  resistance  associated  with  the  large-area  cells  does  not 
influence  the  Vg^  measurement  and  can  be  rendered  insignificant  in  the 
Isc  measurement  by  applying  an  appropriate  reverse  bias  across  the 
cells.  In  this  measurement  the  cells  were  maintained  at  T=25°C  by 
placing  them  on  a  Temptronic   thermal    chuck. 

To  determine  N^^,  the  capacitance  C  of  the  completed  cell  was 
measured  fbr  different  values  of  the  reverse  bias  Vr.  Assuming  a  one- 
sided (asymmetrical)  junction,  we  derived  the  base  doping  density  N,, 
from  the  slope  of  the  1/C^  versus  V,^  plots  [69].  These  measurements 
indicated  slight  variations  of  N^^  among  the  cells,  given  in  Tables  4.1 
and  4.2,   that  were  accounted   fbr  in  the  determinations  of  L  . 

We  see  in  Table  4.1  that  the  anall-area  active  cells,  which  are 
actually  n'^pp^  structures,  have  longer  diffusion  lengths  than  do  both 
the  control  and  standard  cells.  On  the  average,  L^  in  the  active  cells 
is  about  a  factor  of  three  longer  than  in  the  standard  cells.  The 
average  L^  in  the  control  cells  is  about  a  factor  of  six  shorter  than  in 
the  active  cells.     These  results  suggest  that  the  Al-BSF   process  getters 
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intragrain  impurities,  reducing  the  defect  density  by  more  than  an 
order-of-magnitude.  We  discuss  this  and  other  possible  explanations  fbr 
the  results  later. 

Measurements  of  large-area  cells,  the  results  of  which  are  given  in 
Table  4.2,  also  indicate  that  the  conventional  Al-BSF  process  is 
beneficial.  However,  the  increase  in  L^  in  the  large-area  cells  is 
seemingly  smaller  than  that  in  the  small-area  cells  indicated  in 
Table  4.1.  Typically  L^  is  about  50%  longer  in  the  active  cells  than  in 
the  standard  cells.  The  discrepancy  is  possibly  due  to  the  fact  that 
large-area  cells  were  measured  under  light-bias,  i.e.,  the  base 
recombination  current  was  extracted  from  Isc-Vqc  measurements,  whereas 
the  3Tiall-area  cells  were  measured  in  the  dark,  i.e.,  the  dark  I-V 
characteristic  was  used.  We  discuss  this  discrepancy  later  in 
conjunction  with  the  theory  [57]  that  light-bias  effectively  negates 
defect  clusters  in  polysilicon  and  increases  the  actual  minority-carrier 
di  f fusion  length. 

We  also  observe  in  Table  4.1  that  the  average  L^  in  the  control 
cells  is  about  50%  shorter  than  in  the  standard  cells.  Since  the  front 
and  back  surfaces  of  the  control  cell  were  masked  by  an  oxide  (1700A) 
during  the  800°C  heat  treatment,  the  shorter  L^  in  the  control  cells 
could  be  due  to  the  enission  of  impurities  from  the  grain  boundaries 
into  the  intragranular  regions.  Such  emission,  which  has  been  suggested 
[70],  would  also  occur  in  the  active  cells,  but  apparently  the  emitted 
impurities  as  well  as  the  intragrain  defects  are  gettered  to  the  Si-Al 
interface. 

With  regard  to  grain-boundary  passivation  resulting  from  the  Al-BSF 
process,    qualitative   results    from   EBIC  measurements   performed  on  large- 
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area  cells  show  no  significant  reduction  of  S^?  .^  .  This  conclusion  is 
based  on  a  comparison  of  the  percentages  of  background  E8IC  current  lost 
at  the  grain  boundary,  averaged  over  ten  randomly  selected  grain 
boundaries  in  the  active  and  in  the  control  cells.  This  result  is  not 
unexpected     since    we    are    using    thick    wafers  (W„   ^  350  pm) ,   and    it    is 

D 

unlikely  that  the  aluminum  could  diffuse  sufficiently  up  grain 
boundaries  to  significantly  passivate  than  in  the  vicinity  of  the  front 
n"^p  junction.  Furthermore  the  formation  of  the  Al-Si  alloy  tends  to 
restrict  the  amount  of  aluminim  entering  the  silicon  [63].  Hence  we 
conclude  that,  for  bulk  polysilicon,  aluminum  produces  no  significant 
grain-boundary  passivation  when  diffused  from  the  back  at  a  high 
temperature.  However  it  does  seem  to  getter  intragrain  impurities, 
thereby     increasing     L^.  Hence     an     optimal      Al-BSF     process     might 

substantially  increase  bulk  polysilicon  solar  cell  efficiencies,  perhaps 
by   producing  an  effective  BSF   [40]. 

4.3.2     Thin-Film  Process 

Since  the  aluminum  diffusion  from  the  front  is  done  at  a  low 
temperature  prior  to  the  phosphorus  (n"^)  diffusion,  we  observe  no 
gettering  attributable  to  it.  This  could  possibly  be  because  most  of 
the  impurities  are  immobile  at  temperatures  around  450°C.  However, 
significant  grain-boundary  passivation  is  achieved  via  this  process  as 
inferred  from  EBIC  measurements.  These  measurements  were  made  on  a 
scanning  electron  microscope  using  a  beam  energy  of  35  KeV,  which  yields 
a  maximun  electron  penetration  depth  of  about  8  pm  [21].  Hence  the 
measurements  reflect  recombination  properties  of  the  grain  boundaries 
within   about  20  pm    from  the    front   surface. 
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In  the  EBIC  images  of  the  top  surface,  we  have  observed  random 
clustering  of  colored  dislocations  in  some  grains.  We  attribute  this  to 
preferential  aluminum  diffusion  down  dislocations  that  might  possibly  be 
aluminum-induced.  The  junction,  however,  is  not  shorted  because  it  is 
formed  after  the  low- temperature  aluminum  diffusion. 

To  illustrate  the  grain-boundary  passivation  effected  by  aluminim 
diffusion,  we  compare,  in  Fig.  4.2,  the  percentages  of  background  EBIC 
lost  at  the  grain  boundaries  in  active  and  control  cells  (fabricated  on 
sequential  wafers)  and  in  the  standard  cell.  The  percentage  of  EBIC 
lost  at  a  grain  boundary  is  experimentally  determined  as   follows. 

The  maximum  EBIC,  I'^gj^,  at  two  points  far  away  from  the  grain 
boundary  and  on  either  side  of  it  is  measured  directly.  Note  that  the 
points  should  be  chosen  far  enough  away  from  the  grain  boundary  so  that 
the  grain-boundary  recombination  does  not  limit  I^JIJ^.     The  EBIC,    Irorr 

Ld  iL  to  iL  ' 

when  the  beam  is  positioned  at  the  center  of  the  grain  boundary  is  also 
measured.  The  percentage  current  lost  at  the  grain  boundary  is  the 
ratio  of  the  current  lost  at  the  grain  boundary  (given  by  the  difference 
between  the  average  of  the  maximun  EBICs  and  the  EBIC  measured  at  the 
grain  boundary)    and  the  average  of  the  maximum  EBICs. 

In  Fig.  4.2  we  see  a  clear  indication  that  the  aluminim  diffusion 
reduces  the  current  loss  at  the  graip  boundaries,  which  implies  that  it 
decreases  S^^^^^  near  the  n"^p  junction.  Furthermore  it  seems  in  some 
cases  to  reduce  the  scatter  in  the  recombination  properties  of  the  grain 
boundaries.  Additional  support  for  these  qualitative  conclusions  is 
provided  in  Table  4.3  where  we  present  similar  results  of  EBIC 
measurements  performed  on  active,  control,  and  standard  cells  fabricated 
on  nonsequential    wafers.     We  emphasize   that  these   qualitative  EBIC 
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Fig.    4.2     Statistical    distribution  of  the   percentage   background   EBIC 
lost  at  the   grain   boundaries   in   the  active,   control,   and 
standard  cells    (beam  energy=35   keV).      Each   box   refers    to  cells 
from  sequential   wafers.      Thick  arrows   indicate   the  average 
percentage   EBIC   lost  at  the   grain   boundaries   in  each  cell.      The 
aluminum  diffusion    (annealing)   time  and  subsequent  silicon-etch 
are   given   for  active    (and  control)   cells. 
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results  merely   indicate   a   trend    statistically.      The   actual    decrease    in 

PR 
^n(eff)    P'*°^"'^®^    ^"<^    ^^^    significance    with   regard    to    an    actual    thin- 
film  cell    perfbrmance  are  now  considered. 

In    order    to    interpret    the    EBIC    measurements    quantitatively    and 

PR 
evaluate  S  ,    ^^   at  passivated  and  unpassivated  grain  boundaries,  we  use 

the   solution  to   the  minority-carrier  (electron)    transport  problem   in   the 

p-type    base    outlined    in    Chapter    3.       We    will     use    the    results   of  the 

numerical    solution    for  the   simple  case  of  a  circular  generation    source, 

the    utility    of    which    was    discussed,    to    demonstrate    the    efficacy    of 

aluminun   passivation. 

The  EBIC  response  as  a  function  of  the  distance  of  the  electron 
beam  from  the  grain  boundary  is  measured  experimentally  using  the  setup 
shown  in  Fig.  4.3a.  The  electron  beam  traverses  the  grain  boundary  at  a 
known  speed.  The  Keithley  427  amplifier  converts  the  current  to  a 
voltage  which  is  recorded  on  the  x-y  plotter.  A  typical  plot  obtained 
using  this  setup  is  shown  in  Fig.  4.3b.  Note  that  it  is  important  to 
calibrate  the  sweep-rate  of  the  x-y  plotter  and  the  scanning-rate  of  the 
electron  beam  precisely  in  order  to  relate  the  x-axis  of  the  plot  in 
Fig.   4.3b  to  actual    distances  in   the   sample. 

The  theoretically  calculated  values  of  EBIC  current  as  a  function 
of  the    beam    position    are    plotted    in    Fig.    4.4    for   different  values   of 

PR 

S  /    -«    and    for  the  case  of  a  circular   (spherical    in   three  dimensions) 
generation    source.      The   value   of  L^   assumed   corresponds   to   that  of  an 

PR 

active     cell.  We      have     assumed     that  S  ,    ^^.    is     constant,      i.e., 

n  ( e  f  f )  

independent  of  the  excitation  at  the  grain  boundary,  which  is  not  valid 
in  general  as  shown  in  Chapter  2.  Superimposed  on  the  theoretical  plots 
in  Fig.   4.4  are     the  measured     EBIC  responses     obtained     on   two     typical 
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Table  4.3  Measured  EBIC  losses  at  grain  boundaries  in  thin- film- 
cell-related  active,  control,  and  standard  n'^p  cells 
fabricated     on      non- sequential      wafers.  The     aluminum 

diffusion   (annealing)    time  and   subsequent   silicon  etch  are 
given    for  the  active  (and  control)   cells. 


Cell 


ACTIVE  (6  hr;  2  ym) 
ACTIVE  (6  hr;  1  ym) 
ACTIVE  (3hr;  0.5  um) 


CONTROL  (6  hr;  2  ym) 
CONTROL  (6  hr;  1  ym) 
CONTROL  (3  hr;  0.5  ym) 

STANDARD 
STANDARD 


EBIC  Current  Loss  at  Grain  Boundary  {%) 


28 
28 
25 

29 
26 
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grain    boundaries   in    an    active   cell.      By   comparing   the    theoretical    and 

PR  /I 

experimental    responses,    we    infer    that  S^°^^^.   =     5  x  10^  cm/ sec    at   the 

passivated   grain   boundaries  near    the    n'^p  junction.      This  value   is  more 

than    an    order-of-magnitude    lower    than    the    value    obtained    for    the    same 

two  grain  boundaries  (in  a  sequential    control    cell)   without  passivation, 

which    follows    from   Fig.    4.5.      In    Fig.   4.5    calculated   and  measured   EBIC 

responses     from     the    control     cell     are     plotted.        These     plots     imply 
PR  fi 

^n(eff)   ^    ^^      cm/sec,    near    the    kinetic    limit    velocity    S^^    [31],    at 

typical     unpassivated    grain    boundaries.       Similar    analyses    reveal     that 
PR 
^n(eff)   ^^  ^  typical   grain  boundary  in  a   standard  cell    is  comparable  to 

that  in  a  control    cell  . 

We   now   support  the   assumptions  of  constant  s'"^    .^   and   a   spherical 

generation    source,    in    conjunction    with    the    results    in    Chapter   3.       In 

Chapter     3     it     was     shown     that     the     nonlinearity     in  S^^    ^«    is     not 

n(eff) 

significant  until  the  beam  is  sufficiently  close  to  the  grain 
boundary.     This  is  particularly  true  in  the  vicinity  of  passivated  grain 

PR 

boundaries,  where  S^^^^^  is  low.  The  reason  for  this  is  that  the 
minority  electron  density  at  the  grain  boundary  is  not  high  enough  to 
significantly  change  the  barrier.  Consequently  S^^  ..,  is  nearly 
constant  in  the  EBIC  measurements  made  and  analyzed  in  this  chapter. 
Furthermore  we  concluded  in  Chapter  3  that  the  spherical  generation 
source  is  a  good  representation  of  the  actual  electron-beam-induced 
generation  provided  we  interpret  data  (as  we  did  here)  taken  about  2Zq 
away  from  the  grain  boundary  (see  Fig.  3.13),  where  z^  is  the  maximum 
penetration  depth  of  the   primary  electrons  (z     =  8  pm  at  E^  =  35  keV). 

In    order     to    estimate     the    reduction     in    the    grain-boundary    trap 
density     N^-j-^^^^^     corresponding  to  the     decrease   in     S^^    ... 
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resultingfrom  the  aluminLim  diffusion,  we  use  the  analytic  model 
developed  in  Chapter  2.  It  was  shown  that  for  a  monoenergetic  density 
of  trap  states  located  at  the  center  of  the  energy  gap,  s'"?  ..>  is 
insensitive  to  the  grain-boundary  excitation  for  sufficiently  low 
"^SKeff)   ^"^     sufficiently     low     excitation.         For     these     conditions, 

^n(eff)  '  ^ST(eff)-  ^^  ^^"^  '''^''  ^^^^  ^ST(eff)  ^  ^-5  x  10^^  cm'^  at 
typical  unpassivated  grain  boundaries,  and  it  is  more  than  an  order-of- 
magnitude  lower  at  typical  passivated  grain  boundaries.  The  low  value 
°^'^ST(eff)  ^~-^°  ^"^'  ^  ^^  ^^  passivated  grain  boundary  reveals  that 
the  aluminum  diffusion  reduces  the  grain-boundary  trap  density  toward  a 
physical  lower  limit,  and  that  further  reductions  in  s''^  .  >,  which  is 
defined  by  the  grain-boundary  potential  barrier,  can  be  achieved  only  by 
increasing  N^^  (see  Chapter  2). 


4.4     Pi  scussion 

In  this  chapter  we  have  demonstrated  experimentally  that  aluninum 
diffusion  can  be  beneficial  to  both  bulk  and  thin- film  polysilicon  solar 
cells. 

With  regard  to  bulk  polysilicon  cells,  we  found  that  the 
conventional  Al-BSF  (800°C)  process  [64]  produced  concomitant  gettering 
of  prevalent  intragrain  impurities,  in  accord  with  previously  reported 
results  [63].  The  gettering,  which  increased  the  average  minority 
electron  lifetime  (T^^L^/D^lin  p-type  Wacker  (cast)  polysilicon 
significantly,  was  detected  by  dark  current-voltage  and  reverse-bias 
capacitance  measurements  of  n^p  solar  cells.  We  also  found  evidence 
that  the  grain  boundaries  emit  impurities  into  the  intragranular  regions 
during     the     high- temperature     heat     treatment,     but     that     the     emitted 


no 


impurities  as  well  as  the  intragrain  defects  are  effectively  gettered  to 
the  Si-AI  interface.  EBIC  measurements  revealed  that  the  A1-BSF  did  not 
passivate  grain  boundaries  near  the  front  junction  because  of  the 
thickness  of  the  bulk   polysilicon. 

It  has  been  theorized  that  L^  increases  with  the  illumination  level 
because  defect  clusters  are  negated  by  high  carrier-injection  levels. 
Comparison  of  L^  in  standard  and  control  cells  measured  under  light-bias 
(large-area  cells  in  Table  4.2)  and  in  the  dark  (small-area  cells  in 
Table  4.1)  support  this  theory.  For  example,  the  average  L^  in  the 
standard  cell  is  about  a  factor  of  two  higher  when  measured  using  light- 
bias  than  when  measured  in  the  dark.  Since  the  defect  densities  in  the 
active  cells  are  reduced  due  to  the  gettering,  no  such  increase  in  L^ 
with  illumination  is  observed  (compare  results  in  Table  4.1  with  those 
in  Table  4.2). 

With  regard  to  thin- film  polysilicon  cells,  we  found  that  a  low- 
temperature  (450°C)  aluminun  diffusion  from  the  front  surface  of  the 
cell  produced  significant  passivation  of  grain  boundaries.  The 
passivation  was  experimentally  assessed  using  EBIC  measurements  which 
showed  a  reduction  in  the  current  lost  at  an  al  uminum-di  f fused  grain 
boundary  in  comparison  with  one  without  aluminum.  Typically,  the  grain- 
boundary  recombination  velocity,  which  is  unaltered  by  EBIC  excitations 
far  away  from  the  grain  boundary,  is  reduced  from  greater  than 
106  cm/ sec  to  less  than  IQS  cm/ sec  by  the  aluminum  diffusion.  These 
values  of  S^^g^^^  near  the  n"^p  junction  were  determined  by  interpreting 
measured  EBIC  responses  using  a  computer-aided  nimerical  simulation  of 
the  underlying  electron  transport  in  the  quasi-neutral  base  of  the  n^p 
cell      developed     in     Chapter     3.         Using     the     analytic     grain-boundary 


Ill 


recombination  model  developed  in  Chapter  2,  we  found  that  the  grain- 
boundary  trap  density  was  also  lowered  by  more  than  an  order-of- 
magnitude  by  the  aluminun  diffusion.  The  value  of  N^y,  ^^.  at  a  typical 
Al-passivated  grain  boundary  is  -10^°  cm"^,  which  is  probably  low 
enough  to  enable  the  virtual  suppression  of  grain-boundary  recombination 
losses  in   thin- film  cells  by   proper  choice  of  N^^^   (see   Chapter  2).     For 

example,     from    Fig.    2.9    we    infer    that    if   U..    were    10^^  cm"^,   S^?    ^^, 

^'^  n(eff) 

would  be  near  10^  cm/ sec,  which  would  render  the  grain  boundary 
virtually  inactive  electrically.  We  note  that  the  grain-boundary 
passivation  produced  by  the  alumimjm  diffusion  is  distinct  from  any 
passivation  due  to  the  phosphorous  diffusion,  which  would  occur  in  the 
standard  and  control  cells  as  well  as  the  active  ones,  but  which  is 
probably  negligible  [16]. 

To  exemplify  the  real  benefit  of  the  aluminum  diffusion  to  the 
efficiency  of  a  thin- film  polysilicon  solar  cell,  we  estimate  the 
increase  in  Vqq  that  would  result  from  the  inferred  reduction  in 
NsT(eff)  "oted  above.  It  has  been  shown  [33]  that  the  predominate 
recombination  current  component  in  a  thin- film  cell  may  well  be  that 
from  the  portion  of  the  grain  boundary  within  the  junction  space-charge 
region,  I^g  ,  which  is  proportional  to  MsT(eff)  ^^^-  ^^^^  in  this 
case,     the    large     reduction     in  N^.^,    ...   would     increase     y^^r-    by    about 


'ST(eff)     "^"'^        M,^,,ca:,c       VQ(> 

rSCf 


(kT/q)ln(50)   =  100     mV,      provided  I^^"^  were     still      predominant     in     the 
grain-boundary-pa ssivated  cell . 

In  contrast  to  other  grain-boundary  passivation  techniques  [17-19], 
we  expect  the  Al-diffused  grain  boundaries  to  be  quite  stable.  The  low- 
temperature  aluminum  diffusion  was  done  prior  to  the  high- temperature 
phosphorus   (n"^)    diffusion,    which  therefore    stabilizes  the  redistributed 
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aluminun.  The  grain  boundaries  are  effectively  passivated  even  after 
the  phosphorous  diffusion,  which  implies  that  no  significant  out- 
diffusion  of  aluminum  from  the  grain  boundaries  occurs,  possibly  because 
of  a   high   segregation  coefficient   for  aluminun  at  the  grain  boundaries. 

We  emphasize  that  the  aluminum  processes  described  in  Section  4.2 
are  not  necessarily  optimal.  Additional  experiments  are  needed  to 
determine  optimal  temperatures  and  times  that  maximize  the  gettering 
efficacy  of  aluminin  in  bulk  cells  and  the  grain-boundary  passivation  in 
thin- film  cells.  Furthermore  in  the  thin- film  process,  the  amount  of 
aluminin  evaporated  onto  the  front  surface  prior  to  the  diffusion  and 
the  subsequent  silicon-etch  must  be  properly  specified.  The  main 
objective  of  this  research,  however,  is  not  to  define  optimized 
processes,  but  to  indicate  the  potential  improvements  afforded  by 
aluminim  diffusion  in  both  bulk  and  thin- film   polysilicon    solar  cells. 

Although  we  have  suggested  that  gettering  is  the  most  likely 
mechanisTi  that  produces  the  observed  increases  in  the  diffusion  length, 
i.e.,  the  decreases  in  the  base  recombination  current  J^q,  we  cannot 
rule  out  the  possible  presence  of  the  BSF  effect  [40].  Preferential 
dissolution  of  silicon  can  occur  when  aluminim  is  evaporated  onto 
silicon  and  heated  at  temperatures  greater  than  700°C  [28].  This 
dissolution,  which  occurs  most  rapidly  at  dislocations  due  to  we^k 
bonding,  can  penetrate  deeply  into  the  silicon  [28].  As  a  result,  a 
nonuniform  and  nonplanar  BSF  is  produced,  which  could  cause  an  apparent 
increase  in  the  minority-carrier  diffusion  length  as  implied  by  Jgg. 
Hence,  although  we  surmise  the  Al-BSF  getters  intragrain  impurities,  we 
stress  that  our  results  could  be  due  to  the  formation  of  a  nonuniform 
BSF. 
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In  the  next  chapter  we  examine  the  mechanisms  of  gettering  and 
nonunifbm  BSF  due  to  preferential  dissolution  of  silicon  in  greater 
detail  in  an  effbrt  to  isolate  the  actual  mechanism  that  is  responsible 
for  the  observed  improvements  in  the  bulk   polysilicon  cell    performance. 


CHAPTER   5 
ALUMINUM  GETTERIN6   IM    (CAST)    POLYSILICON 


5.1  Introduction 
In  Chapter  4  we  presented  experimental  evidence  fbr  benefits  to  the 
performance  of  poly  silicon  solar  cells  afforded  by  alumintm  diffusion. 
We  found  that  the  conventional  Al-BSF  process  [51]  resulted  in  a 
reduction  in  the  base  saturation  current  density,  JgQ,  which  we 
attributed  to  gettering  of  impurities  from  the  base  region  by  the 
aluminun  alloying.  Another  possible  explanation  for  the  reduction  in 
Jbo  could  be  the  formation  of  an  effective,  but  nonuniform,  nonplanar 
back- surface  field  (BSF)  [40]  due  to  the  preferential  dissolution  of 
silicon  during  the  high- temperature  process  [28].  Although  the  results 
in  Tables  4.1  and  4.2  imply  that  the  BSF  is  inconsequential  in  bulk 
cells  since  L^  <  Wg,  we  consider  it  because  the  preferential  dissolution 
of  silicon  can  cause  deep  penetrations  of  aluminun  into  the  base  of  the 
cell,  resulting  in  an  effective  BSF  in  some  portions  of  the  cell,  and 
because  the  L^  that  we  measured  is  an  average  value  over  several  grains 
in  which  the  local   value  may  vary   significantly. 

If  the  aluminuTi  penetrations  are  substantially  deep  (~  Wg),  and/or 
if  L^  is  sufficiently  long  (~  Wg),  then  the  BSF  effect  can  be  prevalent 
over  a  significant  area  of  the  cell,  and  the  reduction  in  JgQ  observed 
in  Chapter  4  could  be  due  to  this  nonuniform  BSF  rather  than  to 
gettering.  The  main  objective  of  this  chapter  is  to  experimentally 
identify  the  actual  mechanisn  that  i  s  responsible  for  the  reduction  in 
JBO- 
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Most  of  the  previous  research  regarding  the  a1  uminum-sil  icon  system 
has  been  focused  on  low- temperature  processing,  e.g.,  sintering,  because 
of  its  extensive  use  in  the  fabrication  of  silicon  integrated 
circuits.  Gettering  of  impurities  by  aluminLm  has  been  observed  at 
temperatures  below  the  eutectic  temperature  (577°C)  [71,  72].  Thompson 
and  Tu  [72]  found  that  annealing  a  silicon  wafer,  which  has  a  layer  of 
aluminum  on  its  surface,  at  540°C  for  100  hours  results  in  the  removal 
of  copper  and  gold  from  the  silicon.  As  intimated  in  Chapter  4, 
aluminum  also  preferentially  diffuses  down  dislocations  (aluminum 
"spiking")  when  heated  at  below-eutectic  temperatures,  resulting  in  deep 
(>  10  um)  penetrations  into  the  substrate  [65,  73,  74].  If  the  aluminum 
"spikes"  down  dislocations  in  shallow-junction  devices,  it  can  increase 
leakage  currents  and  could  even  result  in   shorted  junctions  [73]. 

The  gettering  as  well  as  the  penetrations  into  the  substrate  by 
aluminum  occur  at  above-eutectic  temperatures  also  [28,  63].  However 
the  kinetics  of  these  mechanisms  are  different  at  low  and  high 
temperatures;  solid-phase  diffusion  occurs  at  low  temperatures,  and 
liquid-phase      diffusion      occurs      at      high      temperatures.  Another 

distinguishing    feature   is  that  the  liquid-phase   diffusion    proceeds  at  a 
much   faster  rate  than  the   solid-phase  diffusion  [62]. 

In  this  chapter  we  will  study  the  impurity  gettering  as  well  as  the 
dissolution  of  silicon  into  the  aluminun  during  the  high- temperature 
processing  of  silicon  with  a  layer  of  almiinun  on  its  surface.  The 
removal  of  impurities  from  the  base  region  of  a  p-n  junction  diode, 
i.e.,  gettering,  results  in  an  increase  in  the  minority-carrier 
lifetime,  and  a  consequent  reduction  in  JgQ.  The  dissolution  of  silicon 
into    the    aluminun,    which    occurs    preferentially    [28],    may    cause    deep 
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penetrations  of  aluninum  into  the  silicon  resulting  in  an  effective 
back- surface  field.  If  the  BSF  effect  is  significant,  then  it  reduces 
JgQ  as  does  gettering.  The  formation  of  an  effective  BSF  hence  may  be 
wrongly  interpreted  as  an  increase  in  the  minority-carrier  diffusion 
length.  A  good  understanding  of  the  Al-Si  system  when  subjected  to 
above-eutectic  temperatures  is  therefore  essential  to  distinguish 
between  these  two  mechanisms. 

We  provide,  in  this  chapter,  a  qualitative  description  of  the 
mechanisms  of  gettering  and  the  formation  of  a  nonuniform  BSF  due  to  the 
dissolution  of  silicon.  In  order  to  identify  the  mechanisn  (gettering 
or  nonuniform  BSF)  that  is  responsible  for  the  reduction  in  JgQ,  we 
determine  the  minority  electron  diffusion  length  L^  in  al  urn  in  urn- di  f  fused 
(n  pp  )  polysilicon  cells  before  and  after  lapping  off  the  back  p^p 
junction.  The  values  of  L^  are  determined  from  forward  current-voltage 
and  reverse-bias  capacitance  measurements  on  large-area  Wacker  (cast) 
polysilicon  cells.  We  also  investigate  the  benefits  of  the  Al-BSF 
process  in  cells  made  from  semiconductor-grade  silicon  substrates  as 
well  as  from  solar-grade  [29]  silicon  substrates.  The  motivation  for 
doing  this  is  to  determine  if  the  benefits  of  the  Al-BSF  process  are 
unique  to  the  (cast)    polysilicon. 

Our  results  indicate  that  the  reduction  in  JgQ  which  results  from 
the  conventional  Al-BSF  process  [61]  is  predominantly  due  to  impurity 
gettering.  The  BSF  effect  due  to  the  penetration  of  aluminum  into  the 
silicon  does  not  significantly  influence  JgQ.  No  significant  increase 
in  L^  was  observed  in  the  active  (n'-pp^)  cells  fabricated  on 
semiconductor-grade  and  solar-grade  silicon,  which  suggests  that  the 
impurity    gettering    is    unique    to    the    (cast)     fwlysil  icon.      The    results 
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imply  that  the  polysililcon  contains  unique  impurities  that  limit  L  ; 
they  are  not  prevalent  in  semiconductor-grade  nor  solar-grade  silicon. 
Data  obtained  from  Wacker  Siltronic  Corporation  [75]  indicate  that  the 
concentration  of  oxygen  in  the  cast  polysilicon  is  about  four  orders  of 
magnitude  higher  than  that  in  semiconductor-grade  and  solar-grade 
silicon.  The    presence    of   oxygen    can    cause     short    minority-carrier 

lifetimes  [76]  due  to  the  creation  of  donor- type  impurities  during  heat 
treatment.  Hence  we  surmise  that  aluminLn  getters  oxygen  and/or  oxygen- 
related  defects  during   the   high- temperature   process. 

5.2     Description  of  the  Mechanisms 
5.2.1     Gettering 

Gettering  can  be  defined  as  the  process  of  removing  impurities  and 
defect  nuclei  from  crystals  and  wafers.  The  application  of  suitable 
thermal,  mechanical,  or  chemical  forces  to  the  wafer  can  cause  the 
transport  of  metallic  impurities  and  defect  nuclei  (such  as  oxygen)  from 
the  wafer  to  an  external    source   such  as  the  wafer   surface  [60]. 

The  effectiveness  of  gettering  depends  upon  two  criteria:  1)  the 
temperature  must  be  high  enough  that  impurities  can  diffuse  readily 
through  the  silicon  and  2)  the  impurities  must  be  captured  (soluble)  at 
a  noncritical  region  of  the  wafer  and  must  not  be  released  back  into  the 
wafer  on   subsequent  heat  treatments. 

Gettering  can  be  achieved  by  the  introduction  of  stresses  on  the 
surface  of  the  wafer  followed  by  heat  treatments  at  elevated 
temperatures  [60].  Metallic  impurities  such  as  iron,  copper,  and  gold 
can  be  gettered  by  the  predeposi tion  of  phosphorous,  boron,  arsenic,  or 
galliim    [60].      As    a   result  of  the   lattice    stress  caused   by   the    dopant 
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diffusion,  dislocations  are  created  at  the  surface  of  the  wafer  [50]. 
Impurities  diffuse  towards  these  dislocations,  which  are  created  away 
from  the  active  regions  of  the  device,  and  are  captured  by  the 
dislocation   sites. 

Another  model  for  gettering  has  been  propsed  by  Goetzberger  and 
Shockley  [77].  According  to  this  model,  metal  atoms  diffuse  to  the 
surface  where  they  react  with  the  gettering  agent  by  formation  of 
chemical  compounds.  The  gettering  of  copper  by  phosphorous  has  been 
explained  by  this  model.  Copper  atoms  diffuse  from  within  the  wafer  to 
the  surface  during  phosphorous  or  boron  predeposition  to  form  borates  or 
phosphates  [77,  78].  It  has  been  shown  [79]  that  metallic  impurities 
can  also  be  gettered  by  employing  inorganic  oxides  as  gettering  agents 
at  temperatures  between  800°C-1000°C. 

It  has  also  been  found  [80]  that  copper  and  gold  tend  to  accumulate 
in  the  heavily  doped  regions  of  the  wafer.  This  has  been  attributed  to 
the  high  solubility  and  low  diffusivity  of  copper  and  gold  in  heavily 
doped  layers.  The  reason  for  the  increased  solubility  can  be  explained 
by  an  analogy  to  chemical  reactions.  The  defect  (solute)  density  in  the 
doped  semiconductor  (solvent)  increases  when  the  doping  density 
(solvent)  is  increased  in  accordance  with  the  laws  of  mass  action 
[80,  81].  For  example,  when  copper  and  gold  exhibit  acceptor  behavior 
in  silicon,  an  n+  (phosphorous)  diffusion  has  been  found  to  getter  them 
out  of  the  low-doped  regions  of  the  wafer  to  the  heavily  doped  region 
[80]. 
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5.2.2     Formation  of  a  nonuniform  back-surface    field 

In  our  experiments  described  in  Chapter  4,  pure  aluminum  was 
deposited  onto  silicon  and  annealed  at  800°C  in  a  forming-gas 
atmosphere.  During  this  process  no  melting  occurs  until  the  temperature 
of  the  Al-Si  system  rises  well  above  the  eutectic  temperature  (577°C), 
or  very  close  to  the  aluminum  melting  temperature  (660°C)  [28].  The 
reason  for  this  is  that  there  is  no  material  of  eutectic  composition 
present  initially.  Once  the  aluminum  is  melted,  an  Al-Si  liquid-phase 
of  temperature-dependent  composition  is  formed  comprising  all  the 
deposited  aluminum  and  the  required  silicon,  which  comes  from  the 
substrate.  As  the  sample  is  slowly  pulled  out  from  the  furnace,  the 
composition  of  silicon  in  the  liquid- phase  must  decrease  as  implied  by 
the  phase  diagram  shown  in  Fig.  5.1.  The  excess  silicon  segregates  at 
the  silicon-liquid  interface  and  forms  a  regrowth  layer  doped  with 
aluminum  ( p"^  layer)  [28].  The  regrowth  layer,  the  thickness  of  which  is 
directly  proportional  to  the  weight  of  the  deposited  aluminum  [64], 
forms  the  necessary  p'^p  (high-low)  junction  to  create  a  back-surface 
field.  The  segregation  of  the  silicon  stops  when  the  liquid  is  cooled 
down  to  the  eutectic  point,  and  the  liquid  solidifies  [28].  The  final 
structure  is  illustrated  in  Fig.   5.2. 

During  the  Al-Si  alloying  process,  the  main  factor  that  controls 
the  points  at  which  the  dissolution  of  silicon  begins  is  the  energy 
required  to  remove  the  silicon  atoms  from  the  lattice.  Because  of  the 
low  strain  energy  associated  with  a  dislocation,  removal  of  the  silicon 
atoms  from  the  lattice  requires  the  least  amount  of  energy  at  positions 
where  dislocations  emerge  at  the  surface.  In  other  words  the  most 
loosely     bonded     atoms  (at  dislocations)    are  removed   first  [28].       Since 
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Fig.    5.2     Structure  of  an   aluminum-alloyed    (n"''pp"'")   cell 
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the   grain   boundaries   are   dislocations,    the   removal    of  the    atoms  at   the 

grain    boundaries    takes    place    readily,    resulting    in    aluminur,    "spikes" 

along  the  grain  boundaries  as  reflected   in  Chapter  4.     The   silicon  atoms 

are    taken    into    the  molten    aluminun   until    the    proportion   of  the    silicon 

reaches  about  25%   (at  800°C)   as  inferred    from  Fig.   5.1.     Until    this  melt 

composition   is  achieved,    the    system   is  in   nonequil  ibri  um,   and   hence   the 

silicon    atoms   are    removed    from   the   lattice   at  a   rapid    pace.      Once    the 

melt    composition    corresponding    to    the    temperature    is   reached,    further 

dissolution    takes    place    in   near   equilibrium   conditions,    and    silicon    is 

dissolved  more   or   less   uniformly    from   the    surface    [28].      The   points  of 

initial    attack   (i.e.,   the   points  at  which  the  grain  boundaries  emerge  at 

the    surface)    due    to    the    preferential    dissolution   of   silicon    are   where 

the  deep  penetrations  of  aluminum  into  the   silicon  occur. 

On  cooling  the  wafer,  the  excess  silicon  in  the  Al-Si  alloy,  which 
penetrates  to  different  depths  at  different  points  in  the  silicon 
substrate,  recrystall  izes.  This  recrystall  ized  silicon,  doped  with 
aluminum,  creates  a  nonplanar  back- surface- field  layer.  The  reason  fbr 
the  nonuniformity  is  the  varying  depths  of  penetration  (penetration 
depths  up  to  75  pm  have  been  measured  [28])  of  aluminum  into  the 
silicon.  The  BSF  effect,  if  dominant,  can  reduce  Jg^  [40],  thus 
resulting  in  an  apparent  increase  in  the  minority-carrier  diffusion 
length. 


5.3     Results 


We  fbcus  our  attention  on  two  main  topics  in  this  section: 
1)  Identi  fication  of  the  actual  mechanism,  gettering  or  nonuni  fom  BSF, 
that    is    responsible    fbr    the    reduction    in    Jg^    of   (Wacker)     polysilicon 
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cell  s      fabricated      using      the      conventional      Al-BSF      process.  This 

identification  is  achieved  by  comparing  L^  in  active  (n^pp+)  cells  with 
and  without  the  back  p^p  high-low  junction.  2)  Investigation  of  the 
beneficial  effects  of  the  Al-BSF  process  in  semiconductor-grade  and 
solar-grade  [29]  silicon  cells.  The  purpose  of  this  study  is  to 
ascertain  whether  or  not  the  benefits  that  result  from  the  Al-BSF 
process  are  unique  to  the  (Wacker)  polysilicon  by  comparing  L  in  the 
active  (n"^pp^)  cells  and  in  the  standard  cells  fabricated  on 
semiconductor-grade  and   solar-grade   silicon   substrates. 

5.3.1     Polysil  icon  cell  s 

In  order  to  determine  if  the  reductions  in  Jg^  due  to  the 
conventional  Al-BSF  process  [61]  were  due  to  gettering  or  due  to  a  back- 
surface  field,  we  compare,  in  Table  5.1,  the  values  of  L^  in  large-area 
Wacker  (cast)  polysilicon  cells.  The  cells  were  fabricated  using  the 
bulk  process  outlined  in  Chapter  4,  wherein  we  also  described  the 
technique  to  determine  L^  from  the  short-circuit  current  (I^p)  versus 
open-circuit  voltage  (Vq(.)  characteristic.  The  values  of  L  in  the 
active  (n  pp^)  cells  are  about  50%  higher  than  those  in  the  standard 
cells,  which  is  consistent  with  our  earlier  results.  However  our 
purpose  here  is  not  to  compare  the  active  and  standard  cells,  but  to 
compare  L^  in  al  urn  in  urn- di  f  fused  cells  with  and  without  the  back-surface 
field.  If  the  decrease  in  JgQ  is  caused  by  a  back-surface  field,  and 
consequently,  the  inferred  increase  in  L^  is  only  an  apparent  one,  then 
removing  the  p^p  (high-low)  junction  from  the  al  umin  urn- di  f  fused  cells 
would  increase  J^q. 
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Table    5.1        Measured    elec^tron    diffusion    lengths    in    active    (with    and 
without  the   p     layer)    and   standard  n  p  polysilicon  cells. 


Device 

hA   (cm-3) 

L,  (um) 

ACTIVE 
(n^pp^) 

6  X  10^5 

190 

5  X  10^5 

150 

5  X  10^5 

165 

ACTIVE 
(Back-lapped) 

6  X  10^5 

160 

5  X  10^5 

165 

STANDARD 

6  X  10^5 

105 

5  X  10^5 

100 
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The  removal  of  the  p+p  junction  was  achieved  by  lapping  off 
25-50  vm  of  material  from  the  back  surface  of  the  active  (n^pp^)  cells, 
which  were  fabricated  using  the  conventional  Al-BSF  process.  Aluminum 
back  and  front  contacts  were  then  evaporated  as  described  in 
Chapter  4.  The     cells     were     not     sintered     after     the     back-contact 

evaporation. 

The  values  of  JgQ  in  the  back-lapped  samples  were  measured  to  be 
about  the  same  as  in  the  active  (n^pp^)  cells.  This  suggests  that  the 
back-surface  field  had  no  influence  in  the  reductions  on  JgQ  observed  in 
the  aluminum-diffused  poly  silicon  cells.  We  can  thus  conclude  that  the 
conventional  Al-BSF  process  getters  intragrain  impurities  in  Wacker 
poly  si  1  icon. 

To  further  consolidate  the  fact  that  the  BSF  effect  due  to  the 
dissolution  of  silicon  is  not  important,  we  measured  the  penetration 
depths  of  aluminum  into  the  polysilicon  in  the  active  (n+pp+) 
structures.  In  order  to  measure  these  depths,  we  etched  the  Al-Si  alloy 
layer  from  the  back  of  the  active  cells  by  immersing  them  in  HCl  for  2-3 
days.  We  observed  that  the  back  surface  of  the  n^pp+  cell,  after 
removal  of  the  Al-Si  alloy  layer,  was  pitted  due  to  preferential 
dissolution  of  silicon.  Using  an  optical  microscope,  we  estimated  the 
depths  of  the  pits  to  be  10-20  ym  on  the  average,  which  is  much  less 
than  the  substrate  thickness  (Wg  =350  pm) .  To  produce  an  effective 
BSF,  the  dissolution  must  penetrate  to  near  the  front  junction  [40], 
which  requires  the  pits  to  be  at  least  200-ijm  deep.  Since  this  is 
definitely  not  what  is  observed  experimentally,  we  conclude  that  the 
nonuniform,  non planar  BSF  due  to  the  preferential  dissolution  of  silicon 
i  s  ineffective. 


126 


5.3.2     Semiconductor-grade  and   solar-grade   silicon  cells 

Having  established  that  the  conventional  A1-BSF  process  getters 
intragrain  impurities  in  Wacker  (cast)  polysilicon,  we  next  examine  its 
gettering  efficacy  in  semiconductor-grade  as  well  as  solar-grade  [29] 
silicon  cells.  The  solar-grade  silicon  is  devoid  of  grain  boundaries 
but      has     twin      planes      in      it.  The     motivation      for      fabricating 

semiconductor-grade  and   solar-grade   silicon  cells  is  to  determine  if  the 
gettering  is  unique  to  the  Wacker   (cast)    polysilicon. 

We  fabricated  large-area  active  and  standard  cells  on 
semiconductor-grade  and  solar-grade  silicon  substrates  using  the  bulk 
process  described  in  Chapter  4.  The  semiconductor-grade  silicon 
(Czochral  ski)  wafers  were  obtained  from  Wacker,  and  their  base  doping 
densities  were  determined  to  be  about  3  x  10^^  cm"-^  from  reverse-bias 
capacitance  measurements.  The  base  doping  densities  of  the  solar-grade 
silicon  wafers  were  determined  to  be  about  8  x  10^^  cm"-^. 

The  base  diffusion  lengths,  determined  from  Isc-^oc  characteristics 
(see  Chapter  4),  for  the  active  and  standard  cells  fabricated  on 
semiconductor-grade  silicon  substrates  are  given  in  Table  5.2.  We  find 
no  significant  change  in  L^  due  to  the  Al-BSF  process.  In  Table  5.3  we 
give  values  of  L^  obtained  on  solar-grade  silicon  cells.  We  find  from 
Table  5.3  that,  as  in  the  case  of  silicon  cells,  the  active  (n^pp^) 
cells  have  about  the  same  L^  as  do  the  standard  cells.  Since  L  in  the 
active  cells  fabricated  in  semiconductor-grade  or  solar-grade  silicon  is 
not  longer,  we  conclude  that  the  conventional  Al-BSF  process  does  not 
getter  impurities  in  these  materials,  as  it  does  in  the  polysilicon.  We 
cannot  say  whether  or  not  the  grain  bounadries  having  diffused  aluninum 
facilitate     the     gettering     process,     or  that  the     impurities  in     Wacker 
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Table    5.2  Measured    electron    diffusion    lengths    in    active    (with    and 

without  the   p^  layers)    and   standard  n"'"p  silicon  cells. 


DEVICE 

^A  (^"1"^^ 

Ln  (ym) 

ACTIVE 

3  X  10^5 

200 

3  X  10^5 

265 

3  X  10^5 

280 

ACTIVE 
(Back-lapped) 

3  X  10^5 

200 

3  X  10^5 

280 

STANDARD 

3  X  10^5 

250 

3  X  10^5 

280 

3  X  10^5 

300 
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Table   5.3      Measured   electron   diffusion   lengths  in   active,   control,    and 
standard  n  p  solar-grade   silicon  cells. 


DEVICE 

^AA  (C"^'^) 

Ln  (um) 

ACTIVE 
(n+pp^) 

8  X  10^^ 

175 

8  X  10^5 

145 

CONTROL 

8  X  10^5 

180 

9  X  10^5 

200 

STANDARD 

8  X  10^5 

210 

7  X  10^5 

150 
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(cast)  polysilicon  are  unique,  but  one  or  both  of  these  is  possible. 
Most  probably,  because  L^  is  considerably  shorter  in  the  polysilicon, 
there  are  prevalent  impurities  in  it,  which  can  be  gettered  by  the 
Al-BSF   process  and  which  are  not  prevalent  in  the  other  materials. 

To  identify  the  impurities  present  in  Wacker  (cast)  polysilicon  and 
to  determine  if  they  are  unique,  we  refer  to  data  provided  by  Wacker 
Siltronic  Corporation  [75].  We  compare  the  types  and  densities  of 
impurities  prevalent  in  Wacker  semiconductor-grade  silicon  [75]  and 
(cast)  polysilicon  used  in  our  experiments.  The  impurities  prevalent  in 
Wacker  polysilicon  are  oxygen,  carbon,  titanium,  cobalt,  antimony,  and 
iron     [75].        The    concentration     of    oxygen     in     Wacker     polysilicon     is 

~  5  X  10  cm"  ,  which  is  at  least  four  orders  of  magnitude  greater 
than  that  in  Wacker  semiconductor-grade  silicon  [75].  Because  of  this 
large  difference  between  the  oxygen  content  in  the  polysilicon  and  that 
in  the  semiconductor-grade  silicon,  it  is  likely  that  the  gettering 
observed  in  the  polysilicon  is  due  to  the  removal  of  the  oxygen  atoms 
from  the  bulk. 

When  oxygen-rich  samples  are  annealed  at  temperatures  greater  than 
400°C,  the  oxygen  atoms  are  thermally  activated  into  donor- type 
impurities  [82,  83].  These  donors  cause  a  reduction  in  the  minority- 
carrier  lifetime  [76].  One  possible  reason  for  this  reduction  could  be 
that  the  oxygen  donors  form  effective  recombination  centers,  i.e., 
produce  deep-level  traps  in  the  energy  gap.  A  second  reason  could  be 
that  mobile  impurities,  which  are  gettered  by  oxygen,  are  released  to 
form  effective  recombination  centers  when  oxygen  is  activated  into  the 
donor   state  [76]. 
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We  thus  surmise  that  aluminum  getters  oxygen  and/or  the  mobile 
impurities  possibly  released  by  the  oxygen  during  the  donor-activation 
process.  Because  the  oxygen  content  in  the  polysilicon  is  much  higher 
than  that  in  semiconductor-grade  silicon,  this  conjecture  is  consistent 
with  our  resul  ts. 

We  note  that  although  we  have  suggested  that  aluminum  getters 
oxygen  in  polysilicon,  it  is  not  unlikely  that  other  impurities  present 
in  Wacker  polysilicon  such  as  iron,  titanium,  cobalt,  and  antimony, 
which  fomi  effective  recombination  centers  in  silicon  [84],  are  gettered 
by  aluminum.  Furthermore  Wacker  polysilicon  contains  a  hgher  density 
(-10  cm"-^)  of  carbon,  which  is  at  least  two  orders  of  magnitude 
higher      than      that      in      semiconductor-grade      silicon.  However      the 

recombination  effects  due  to  carbon  are  not  well  understood  [60]  and 
hence  are  not  discussed  here. 

Not  much  information  is  avail  al  be  about  the  type  of  impurities  and 
their  densities  prevalent  in  solar-grade  silicon.  However,  from 
Table  5.3,  we  observe  that  the  measured  diffusion  length  in  standard 
cells  fabricated  on  solar-grade  silicon  is  higher  than  that  in 
polysilicon  (standard)  cells  (see  Table  5.1),  and,  in  fact,  is 
comparable  to  that  measured  in  semiconductor-grade  silicon  (standard) 
cells.  This  suggests  that  the  impurity  (presumably  oxygen)  content  in 
solar-grade  silicon  is  lower  than  that  in  polysilicon,  which  provides  a 
possible  explanation  as  to  why  aluminum  gettering  is  not  observed  in 
solar-grade   silicon. 

In  Table  5.3,  note  that  the  control  cells  made  on  solar-grade 
silicon  show  no  significant  reduction  in  L  from  that  of  the  standard 
cells.      Contrarily   the  control    cells  made  on   polysilicon    substrates   had 
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a  lower  L  than  that  in  the  standard  cells  (see  Tables  4.1  and  4.2), 
which  we  attributed  to  emission  of  impurities  by  the  grain  boundaries 
into  the  grains.  Since  the  solar-grade  silicon  [29]  has  only  twin 
planes  and  no  grain  boundaries  in  it,  there  is  no  emission  of  impurities 
into  the  bulk  region,  and  consequently  there  is  no  reduction  in  L  . 
This  provides  additional  support  for  our  conclusion  that  grain 
boundaries  emit  impurities,  which  has  also  been  observed  elsewhere  [70]. 

5.4     Discussion 


In  this  chapter  we  studied  the  Al-Si  system  when  it  is  subjected  to 
elevated  temperatures.  We  established,  by  means  of  Isc"^OC  "Tiessurements 
on  solar  cells,  that  the  Al-BSF  process  getters  intragrain  impurities  in 
Wacker  (cast)  poly  silicon.  Furthermore,  we  found  that  the  gettering  is 
apparent  only  in  polysilicon  and  not  in  semiconductor-grade  or  in  solar- 
grade  [29]  silicon  cells.  We  then  surmised  that  aluminum  getters  oxygen 
and/or  oxygen-related  impurities,  which  are  present  in  higher  densities 
in   the   polysilicon  than   in   semiconductor  and   solar-grade   silicon  [75]. 

We  provided  a  qualitative  description  of  the  mechanisms  of 
gettering  and  the  formation  of  a  nonuni  form  BSF  due  to  preferential 
dissolution  of  silicon  that  could  possibly  produce  the  observed 
reduction  in  JgQ  in  al  umin  urn- di  f  fused  polysilicon  cells.  Several  models 
[60,  77,  78,  80]  have  been  proposed  in  the  literature  fbr  gettering. 
The  strain- field  model  [60,  78]  postulates  that  during  predeposition  of 
dopants,  dislocations  are  introduced  at  the  surface  as  a  result  of 
lattice  stress.  At  high  temperatures,  the  impurities  are  mobile  and 
hence  they  diffuse  towards  the  dislocations  and  are  captured  by  them. 
The    compound- formation    model     [67]     suggests    that    metallic    impurities. 
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such  as  copper,  fom  a  compound  with  the  gettering  agent  (a  phosphate  is 
formed  in  the  case  of  phosphorous  gettering).  The  impurities  readily 
diffuse  towards  the  surface  since  the  concentration  of  impurities  there 
is  very  low.  The  defect  solubility  model  [80]  postulates  that 
particular  impurities  readily  dissolve  in  heavily-doped  layers  due  to 
their  high  solubility.  For  example,  copper  and  gold,  if  present  as 
acceptor- type  impurities  in  silicon,  tend  to  accumulate  in  the  n+  region 
adjacent  to  the  surface  of  the  wafer  [80],  and  are  thus  gettered  out  of 
the  low-doped  regions  of  the  wafer. 

The  nonuniform,  nonplanar  BSF  is  fomed  as  a  result  of  preferential 
dissolution  of  silicon  during  the  aluminum-diffusion  process.  When 
aluminum  is  alloyed  with  silicon,  the  necessary  silicon  is  obtained  from 
the  substrate  preferentially  from  the  most  weakly  bonded  sites 
(dislocations)  [28].  Consequently  aluminum  penetrates  down  the  grain 
boundaries,  since  grain  boundaries  are  dislocations,  to  varying  depths 
in    the    silicon    resulting    in   a   nonuniform   BSF   layer.      The  effect  of  the 

nonuniform    BSF    layer,    if   predominant,    is    to    cause    a   reduction    in   J,, 

BO' 
which  implies  an  apparent  increase  in  L   . 

We  measured  Isc-^qc  ^"d  reverse-bias  C(V)  characteristics  of  active 
(Al-BSF)  n%p+  polysilicon  cells,  from  which  we  determined  the  electron 
diffusion  length  in  the  base.  In  order  to  prove  that  the  observed 
increase  in  L^  was  real  and  not  merely  apparent,  which  would  be  the  case 
if  the  BSF  was  the  dominant  mechanism,  we  lapped  the  back  p%  junction 
in  the  active  cells  and  determined  L,.  We  fbund  that  values  of  L^  in 
the  back-lapped  cells  were  about  the  same  as  those  in  the  active  (n%p+) 
cells,  which  proves  that  the  Al-BSF  process  getters  intragrain 
impurities  in  Wacker  (cast)    polysilicon. 
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It  should  be  pointed  out  that  the  reduction  in  JgQ  in  active 
(Al-BSF)  cells  is  also  commensurate  with  a  possible  reduction  in  the 
minority-electron  mobility  within  a  diffusion  length  of  the  junction 
[85].  The  reduced  mobility  could  be  caused  by  the  lattice  strain 
produced  during  the  al  umi nun- di  f fusion  process.  If  the  reduction  in  Jon 
were  due  to  a  lowering  of  the  electron  mobility,  then  lapping  off  the 
back  p'^p  (high-low)  junction  should  have  increased  JgQ.  The  reason  for 
this  is  that  lapping  off  the  p'^p  junction  decreases  the  lattice  strain, 
and  consequently  the  mobility  increases.  However,  we  experimentally 
observed  no  change  in  JgQ  when  the  p^p  junction  was  lapped  off  from  the 
active  (Al-BSF)  cells.  This  suggests  that  it  is  unlikely  that  mobility 
reduction  is  an  important  mechanign  that  reduces  JgQ  in  the  active  (Al- 
BSF)   cells. 

Measurements  on  semiconductor-grade  and  solar-grade  [29]  silicon 
cells  revealed  no  significant  increase  in  L^  due  to  the  conventional 
Al-BSF  process  [51].  From  this  we  concluded  that  the  benefits  due  to 
the  Al-BSF  process  were  unique  to  the  polysilicon.  A  possible  reason 
for  this  is  that  Wacker  (cast)  polysilicon  contains  much  higher 
densities  of  oxygen  compared  to  semiconductor  and  solar-grade  silicon 
[75].  The  presence  of  high  densities  of  oxygen  can  reduce  minority- 
carrier  lifetime  in  polysilicon  as  a  result  of  the  formation  of  thermal 
donors  during  heat  treatment  [82,  83].  The  reduction  in  the  minority- 
carrier  lifetime  could  be  due  to  the  formation  of  deep-level  traps, 
which  act  as  effective  recombination  centers,  created  by  the  thermal 
donors,  or  due  to  the  release  of  mobile  impurities  that  become  effective 
recombination  centers  when  oxygen  is  activated  into  the  donor  state. 
The    observed    increase   of  L^    in    the    polysilicon    could   then    possibly   be 
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due  to  the  removal    of  the  oxygen  atoms  or  the  released   impurities  by  the 
al  uminum. 

We  also  found  that  L^  in  the  control  cells  processed  on  solar-grade 
silicon,  which  has  no  grain  boundaries  in  it,  were  about  the  same  as  in 
the  standard  cells  (see  Table  5.3).  The  control  cells  processed  on 
Wacker  polysilicon  showed  a  decrease  of  the  diffusion  length  (see  Tables 
4.1  and  4.2),  from  which  we  concluded  that  grain  boundaries  in 
polysilicon  emit  impurities  at  high  temperatures,  which  is  consistent 
with  [70]   and  with  our   solar-grade-sil  icon  experiments. 

Finally  we  point  out  that  the  characterization  of  the  Al-Si  system 
presented  in  this  chapter  is  not  conclusive.  Further  experiments, 
possibly  using  an  ion  probe  or  Auger  electron  spectroscopy  (AES)  [86], 
are  recommended  to  enable  better  characterization  of  the  Al-Si  system. 
Comparison  of  the  type  and  relative  concentrations  of  impurities 
prevalent  in  the  bulk  of  the  polysilicon  and  at  the  Al-Si  interface 
would   serve  as  an   independent  method  to  detennine  gettering. 


CHAPTER  6 
SUMMARY,    DISCUSSION,   AND   RECOMMENDATIONS 

In     this    dissertation    we     have    developed    an    analytic    model     for 

minority-carrier    transport    in     polysilicon,     and    provided    experimental 

corroboration      for     the     model      using     EBIC     measurements.          The     EBIC 

measurements    were    quantified    by    means    of   a    computer-aided    numerical 

solution   of  the   relevant   transport    problem.      We    have   also    demonstrated 

the     potential     of    aluminLn    diffusion     to     improve     the     performance    of 

poly  si  1  icon   solar  cell  s. 

The  main  accomplishments  of  the  research  are   sunmarized  as   follows: 

1)  An  expression  for  the  effective  recombination  velocity  as  a 
function  of  the  grain-boundary  excitation  (excess  carrier 
density)  was  derived,  and  the  importance  of  including  the 
nonlinearity  in  the  general  solution  of  the  transport 
problem  was  di  scussed. 

2)  An  approximate  analytic  model  for  the  minority  electron 
transport  in  the  base  of  a  n"^p  polysilicon  diode  was 
developed,  and  it  was  shown  that  the  transport  is  influenced 
by  the  grain  boundaries  only  at  low  and  intermediate  forward 
voltages. 

3)  EBIC  measurements  were  used  to  demonstrate  the  dependence  of 
the  grain-boundary  recombination  velocity  on  the  excitation 
as  well  as  to  determine  the  surface- state  density  at  the 
grain  boundary  and  the  diffusion  length  in  the  grain. 
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4)  The  conventional  Al-BSF  process  [51]  was  shown  to  getter 
intragrain  impurities  (in  Wacker  polysilicon)  as  implied  by 
current-voltage  measurements  on  polysilicon  cells,  producing 
at  least  a  factor-of-two  increase  in  the  bulk  diffusion 
length,  vis-a-vis,  a  factor-of- four  increase  in  minority 
(electron)   li  fetime. 

5)  The  gettering  was  found  to  be  peculiar  to  Wacker  (cast) 
polysilicon  cells,  which  we  attributed  to  the  presence  of 
unique  impurities  in  polysilicon  that  are  not  present  in 
semiconductor-grade  nor  in  solar-grade  [29]  silicon  cells. 
The  polysilicon  contained  a  much  higher  density  of  oxygen 
than  semiconductor-grade  silicon,  and  hence  we  speculated 
that  aluminum  getters  oxygen. 

6)  A  low- temperature  aluminum  diffusion  process  was  shown  to 
effectively  passivate  the  grain  boundaries,  producing  more 
than  an  order-of-magnitude  reduction  in  the  defect  density 
at  a  typical    (Wacker   polysilicon)   grain  boundary. 

The    grain-boundary    analysis    in    Chapter    2,    which    is    based    on    an 

effective    monoenergetic     density    of    trap     states     situated    at    midgap, 
PR 

predicted  thatS^^^^^^  is  nonlinear.  The  functional  dependencies  of 
n(eff)  Pf'^dicted  by  the  grain-boundary  analysis  are  valid  even  when 
the  trap  level  is  not  at  the  intrinsic  level,  but  is  within  about  5kT 
above  or  below  it.  For  example,  the  predominant  trap  level  at  a  grain 
boundary  adjacent  to  a  |>type  grain  in  Wacker  polysilicon  has  been 
experimentally  observed  to  be  below  the  intrinsic  level  [34].  The 
effective  recombination  velocity  in  this  case  is  still  dependent  on  the 
excitation     at     the     grain     boundary.  However,      the     onset     of     the 
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nonlinearity  and  the  magnitude  of  the  recombination  velocity  are  both  a 
function  of  the   trap  level. 

In  general  the  grain-boundary  traps  are  not  monoenergetic,  but  may 
be  distributed  in  the  energy  gap  [11],  although  the  exact  nature  of  the 
distribution  is  unknown.  We  recommend  that  Poisson's  equation  in  the 
grain-boundary  space-charge  region  be  solved  subject  to  a  realistic 
distribution  of  trap  states  in  the  energy-gap  at  the  grain  boundary,  and 
thus  used  to  verify  the  accuracy  of  the  results  obtained  using  the 
monoenergetic-density-of- states  assumption. 

The  solution  of  the  minority-electron  continuity  equation  in  the 
base  discussed  in  Chapter  2  is  based  on  the  gradual  approximation  [32], 
and  is  valid  if  and  only  if  r^  >  L^,  which  restricts  the  analysis  to 
relatively  large  grains.  The  base  recombination  current,  derived  from 
the  solution  of  the  continuity  equation,  comprises  the  one-dimensional 
diode  current  and  a  component  of  current  dependent  on  the  grain-boundary 
parameters.  We  found  that  the  latter  component  was  insignificant  at 
high  forward  voltages.  In  order  to  obtain  a  general  transport  model  we 
recommend  that  a  computer-aided  numerical  analysis  of  the  minority- 
carrier  continuity  equation  subject  to  the  nonlinear  boundary  condition 
at  the  grain  boundary  be  done. 

In  Chapter  3,  we  determined  S^^^^^^  at  both  unpassivated  and 
passivated  grain  boundaries  in  Wacker  (cast)  polysilicon  using  EBIC 
measurements.  We  found  that  the  percentage  of  current  lost  at  the  grain 
boundary     is    lowered     as    the    excitation     is    increased,     which     proves 

PR 
that  S^^g^^^    is   nonlinear   as   described   in    Chapter   2.      The    quantitative 

interpretation    of  the   EBIC   measurements   was    facilitated   by   a   computer- 
aided    numerical     solution    of    the     pertinent    two-dimensional     continuity 
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equation  subject  to  the  nonlinear  boundary  condition  at  the  grain 
boundary.  By  comparing  the  theoretical  and  experimental  EBIC  responses, 
we  evaluated  the  trap  density  at  typical  passivated  and  unpassivated 
grain  boundaries,  assuming  the  trap  states  to  be  monoenergetic  at 
midgap. 

We  recommend  that,  in  order  to  determine  the  trap  level  in  the 
energy  gap  from  EBIC  measurements,  the  expression  fbr  S^?  ..«  used  in 
the  nunerical  analysis  be  modified  to  include  the  dependence  on  the  trap 
level.  Our     two-dimensional      numerical      analysis     is     restricted     to 

interpretation  of  EBIC  measurements  in  relatively  large-grain 
polysilicon.  To  enable  the  interpretation  of  EBIC  measurements  in 
small-grain  polysilicon  as  well  as  to  determine  the  bulk  diffusion 
length  directly  from  the  absolute  EBIC  response,  a  refinement  of  our 
analysis  to  include  the  variations  in  the  third-dimension  is  called 
for.  It  should  be  pointed  out  that  the  determination  of  the  diffusion 
length  from  the  absolute  EBIC  response  requires  precise  measurement  of 
the  electron  beam  current. 

In  Chapter  4,  we  demonstrated  the  beneficial  effects  of  aluminum 
diffusion  in  polysilicon.  We  found  that  high- temperature  aluminum 
diffusion  from  the  back  surface  gettered  intragrain  impurities  in  Wacker 
polysilicon.  This  finding  followed  from  measurements  of  the  forward 
current-voltage  and  the  reverse-bias  capacitance  characteristics  of  n^p 
polysilicon  cells.  A  low- temperature  aluminun  diffusion  from  the  front 
surface  of  the  polysilicon  resulted  in  effective  grain-boundary 
passivation.  The  passivation  was  detected  using  EBIC  measurements 
interpreted  via  the  computer-aided  numerical  analysis  developed  in 
Chapter  3. 
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To  optimize  the  aluninum-di  ffusion  processes,  we  recommend 
fabrication  of  bulk  and  thin- film  solar  cells  using  aluminun  diffusion 
at  different  temperatures  and  for  various  times.  Pure  aluminum  is 
recommended  for  the  thin- film  process,  since  it  results  in  deeper  spikes 
[73]  and  hence  effective  grain-boundary  passivation.  The  determination 
of  the  distribution  of  the  aluminun  down  typical  passivated  grain 
boundaries,  possibly  using  ion  probe  and/or  Auger  Electron  Spectroscopy 
(AES)  measurements,  is  called  for  to  expedite  the  optimization  of  the 
passivation   process. 

The  current-voltage  characteristics  measured  on  n+pp+  cells, 
reported  in  Chapter  4,  implied  that  the  conventional  Al-BSF  process  [61] 
does  indeed  getter  intragrain  impurities.  The  benefits  produced  by  the 
Al-BSF  process  is  also  commensurate  with  the  formation  of  an  effective 
BSF  solar  cell.  The  BSF,  which  is  nonuni  fonn  and  nonplanar,  is  fomed 
as  a  result  of  preferential  dissolution  of  silicon  during  the  high- 
temperature  process  [28].  In  Chapter  5  we  examined  the  mechanigns  of 
gettering  and  BSF  formation  using  current-voltage  characteristics  on 
Wacker  (cast)  polysilicon  cells  before  and  after  lapping  off  the  p+p 
junction.  We  established  that  the  observed  improvement  in  the  (cast) 
polysilicon  cell  performance  was  due  to  gettering  and  that  the  BSF 
effect  was  inconsequential.  We  also  demonstrated  that  the  benefits 
produced  by  the  Al-BSF  process  were  unique  to  the  polysilicon,  i.e.,  no 
improvement  in  the  performances  of  semiconductor-grade  and  solar-grade 
silicon  cells  was  observed.  This  suggested  that  the  polysilicon  has 
unique   impurities  in  it  that  are  readily  gettered  by  the  aluminum. 

Data     obtained     [75]     on     Wacker     (cast)      polysilicon     indicate     the 
prevalence  of  a    high  concentration  of  oxygen   in    polysilicon  compared   to 
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that  in  semiconductor-grade  and  solar-grade  silicon.  It  has  been  shown 
[76]  that  minority-carrier  lifetime  is  shortened  due  to  the  presence  of 
excess  oxygen.  We  thus  surmise  that  aluminum  getters  oxygen  from  the 
polysilicon.  However  we  do  not  rule  out  the  possibility  of  gettering  of 
other  impurities  by  aluminum  such  as  iron,  cobalt,  antimony,  and 
titaniun  that  are  also  present  in  (cast)  polysilicon  [75],  and  which 
form  effective  recombination  centers.  To  clarify  the  uncertainty 
involved  as  to  the  type  of  impurities  gettered  by  aluminum,  we  recommend 
that  an  AES  analysis  of  the  Al-Si  system  be  done.  Such  an  analysis 
would   serve  as  an   independent  method  to  detemine  gettering. 

We  discuss  in  conclusion  a   few  general    applications  of  the  research 
described   in    this   dissertation.      The    solution  of  Poisson's  equation    in 
the    grain-boundary    space-charge    region,    i.e.,    the    calculation    of   the 
barrier-height,     is     applicable     also     in     studies     of    majority-carrier 
transport  through  grain  boundries,   e.g.,   SOI  MOSFETs  [5]   and   polysilicon 
resistors    [4].       Since    a    real     semiconductor    surface    is   analogous   to    a 
grain     boundary,     we    can     therefore     anticipate     that,     in    general,     the 
recombination    velocity    at    a    semiconductor    surface    will    be    nonlinear. 
However,     most     analyses    assume     that     the     recombination     velocity     is 
constant  at  a    surface,    the  validity  of  which  is   hence   questionable.     To 
remove    this   uncertainty,    we    recommend   a  comprehensive   investigation   of 
the    recombination    properties   at    passivated    (by    thermal    oxides   or    even 
polysil    emitters    [87]),    and    unpassivated    surfaces    based    on    numerical 
interpretations  of  EBIC  measurements. 


APPENDIX  A 
NUMERICAL   SOLUTION  OF   POISSON'S   EQUATION 
IN  THE   GRAIN-BOUNDARY   SPACE-CHARGE   REGION 


Referring    to    Fig.    2.3,    in    which    the    grain    boundary   is  located   at 

x=0,  we  write  Poisson's  equation   (one-dimensional)    in  the  grain-boundary 

PR 
space-charge  region   (0  <  x  <  W^pn)   as 


d7=7  ■  (A.l) 


where  p  is  the  net  charge  density,  e  is  the  dielectric  constant  of 
silicon,  and  E  is  the  electric  field  in  the  grain-boundary  space-charge 
region. 

The  net  charge  density  in  the  grain-boundary   space-charge  region  is 
given  by 

P  =  q[P   -  N  -  N^^]  (A. 2) 

where  N^^  is  the  grain  (acceptor)  doping  density,  q  is  the  electron 
charge,  and  P  and  N  are  the  hole  and  electron  densities  respectively  in 
the  grain-boundary  space-charge  region.  P  and  N  are  related  to  the 
intrinsic    Fermi    level,    E^- ,    through    the    respective    quasi-Fermi    levels. 


E.(x)   -  E     I 
"i   ^^P| kT —]        '  (A. 3) 


141 


H2 


N  =   n.    expj  ^j    1  (A. 4) 

where  n^.  is  the  intrinsic  carrier  concentration.  The  electric  field  in 
the  grain-boundary  space-charge  region  is  given  by  the  derivative  of  the 
intrinsic  Fermi   level : 


1   '^^i 
^   =   q^       •  (A.5) 


Substituting  (A. 2)  -  (A. 5)  in  (A.l),  we  get 


where 


/E.(x)  -  Eppi        iZ^,,   -  E.(x)\ 


The  electric  field  at  x=0  is  determined  by  the  grain-boundary 
surface-state  charge  through  Gauss'  law  and  provides  one  of  the  two 
boundary  conditions  necessary  to   solve  (A. 6): 


dE, 


x=0  ^ '*-^' 


where  N5j(gf^)  is  an  effective  monoenergetic  density  of  grain-boundary 
states  situated  at  the  midgap  (see  Chapter  2)  and  f  is  the  steady-state 
probability    of    electron     occupancy     for     the     surface     states    given     in 
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(2.3).      The    second  boundary  condition   is  written    in   the   uniformly  doped 
quasi-neutral    p-region  by  assuming  low-injection  conditions: 


dE. 


GB        =   ^        •  (A.9) 


The  differential  equation  (A. 6)  is  not  homogeneous  and  has  a 
nonlinear  boundary  condition,  and  hence  must  be  solved  numerically.  The 
following  nimerical  algorithm,  called  the  n-method  [88],  has  been  used 
to   solve  the   problem. 

Let  y  be  an  approximate  solution  of  the  differential  equation  (A. 6) 
that  also  satisfies  the  boundary  conditions  (A. 8)  and  (A.9). 
Substituting  E.(x)   =  y(x)    in   (A. 6),  we  get 


y"   -7Q(y)   =  0       .  (A. 10) 


A  better  approximation  to  the  solution  of  the  differential  equation 
(A.  10)  is  then  given  by  (y  +  n) ,  where  n  is  a  small  variation. 
Substituting   this  in   (A.  10),  we  get 

(y  +  n)"  -  f  Q(y  +  n)   =  0       .  (A. 11) 

Expanding  (A. 11)  by  a  Taylor  series  and  neglecting  higher-order  terms, 
which  is  justified  if  n  is  small,  we  obtain 


^"  -^7§-  -^y"  -^^iyn     .  (A.12] 
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The  right  side  of  (A. 12)  is  known  from  the  previous  solution  (A.  10), 
and  —  can  be  obtained  by  differentiating  (A. 7).  Thus  (A. 12)  is  a 
second-order  differential    equation   in  n  whose  boundary  conditions  are: 

n'    -   n^^=   -Cy'-  |Q^(y)]       at  x  =  0        ,  (a. 13) 


where  Q^   =  ^N3^(g^^(l-f)    is     the     net    charge     density     at     the     grain 
boundary,  and 


^'    =  °  ^^  ^   ^  W3^CR       •  (A.14) 

The  boundary  condition  (A. 13)  was  obtained  by  substituting  E-  =  (y  +  n) 
in  (A. 8),  expanding  by  a  Taylor  series,  and  neglecting  higher-order 
terms. 

The    differential    equation    (A. 12)     for    n    i  s   a    linear,    second-order 
equation  with  linear  boundary  conditions,    and  is   solved   using  a    finite- 
difference  approach   using   the   Harwell    subroutine   DD02AD   [37].      The   grid 
size     chosen     was     about     40     A,     which     yielded     good     convergence     and 
satisfactory    results.       Starting     from    an    initial     solution    of   E^- ,    the 
numerical    analysis  yields   the   correction    term   n,    which   is   added   to   the 
initial     solution    to    obtain    the    new   value    of  E- .      The   coefficients   in 
(A. 12)    are    evaluated    using    the    new   value   of   E • ,    and    (A. 12)    is    solved 
again    to   obtain    a   new   value    of  n,    and    so    on.      Convergence    is  reached 
when    the   correction    term   n   is   less   than   a   certain   minimun   value   which 
was    chosen    to    be    5    meV    in    our    analysis.       This    nunber    was   chosen    to 
achieve    a    trade-off   between    the    computer    time    required    to     solve    the 
differential    equation    and    the    accuracy    of   the    solution.       The    FORTRAN 
computer   program  is  listed  below. 
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APPENDIX  B 
MATHEMATICAL  JUSTIFICATION  FOR  USING  THE  FOLDING  TECHNIQUE 


Consider  the  n'^p  junction  shown  in  Fig.  3.2  with  an  arbitrary 
volume  distribution  of  the  generation  rate.  The  right-  and  left- side 
grains  are  assumed  to  be  semi-infinite  in  the  positive  and  negative 
x-directions  respectively,  and  the  grain  boundary  is  situated  at  x=0. 
We  will  first  formulate  the  general  problem  without  using  the  folding 
technique,  i.e.,  without  considering  the  grain  boundary  as  a  surface, 
and  then  demonstrate  the  justification    for  using   the    folding   technique. 

In  the  right-side  grain,  the  minority  electron  continuity  equation, 
under  low-injection  conditions,  can  be  written  as 


a^N         a^N         ^\       H         g 

— 4"— 4"— T-^"n^=   0  (B.l) 

3x^         3y^         dz^       L^       ^n 
n 


where  L^   is  the  electron   diffusion  length,   D^   is  the  electron  diffusion 

coefficient,  g,^  is  the  generation  rate  (per  unit  volume   per   second),   and 

Nj^   is   the   electron   concentration    at  any    point   in    the   right-side  grain. 

The    boundary   conditions    fbr    this    problem    are  N|^(z  =  0)    =0  at   the   edge 

of  the  (shorted)   junction   space-charge  region;  N„(z  -v  »)   =  q     since  L     < 

K  n 

Wg,  where  Wg  is  the  width  of  the  base;  Nj^(x  ^  «)  .  o  since  the  grain  is 
semi-infinite  (»  L^)  in  the  x-direction;  and  at  the  edge  of  the  grain- 
boundary  space-charge  region  in  the  right-side  grain  (x=0'^), 
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°n-^=   Sn(eff)    ^^R       '  (B.2) 


PR 
which    defines  S^^g^^^    (see    Chapter    2).       The    boundary    condition    (B.2) 

accounts    for   only    those    carriers    that    flow   towards   the    grain    boundary 

from    the    right-side   grain    and    recombine    there.       Any   carriers   that   are 

trangnitted  through  the  grain  boundary  to  the  left- side  grain  as  well    as 

those    carriers    that    are    generated     there    can    be     accounted     for    by 

formulating   the   problem  in  the  left- side  grain. 

In     the    left-side    grain    we    can    write    the    differential     equation 

analogous  to  that  in   the  right-side  grain: 


8^N,         3^N  3^N  N,         g, 

8x^         3y^         az'^       L'^       ^n 
n 


where  N^  is  the  electron  concentration  and  g^  is  the  generation  rate 
(per  unit  volume  per  second)  at  any  point  in  the  left- side  grain.  We 
have  assumed  that  the  diffusion  lengths  in  the  two  grains  are  equal. 
The  boundary  conditions  for  this  problem  are  N.  (z=0)  =  0;  N  (z+")  =  O; 
N|^(x^)  =  0;  and  at  the  edge  of  the  grain-boundary  space-charge  region 
in  the  left- side  grain   (x=0"). 


^^\  PR 

\-^     =Cff)    \        •  (B.4) 


The    quasi-equilibrium    approximation    ensures    that    N,(0")    =    Nn(0'^)     and 
.GB 
'n(eff) 


PR 

hence     that  S  ,    .«    has    the     same    value    at    x=0^    and    at    x=0~.        Note 
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PR 

that  S^^g^^^  in  (B.2)  and  (B.4)  is  dependent  on  the  excitation  at  the 
grain  boundary  (see  Chapter  2).  The  total  current,  I^^iq,  collected  at 
the  junction  (z=0)   i  s 


"      0°       8Nn       3N, 
^EBIC  =  ^°n     /       /       t-17"^^     ,^^   dy        •  (B.5) 

-00  -00  2=0 


Hence  in  order  to  obtain  the  collected  current,   it  is  necessary  to   solve 
(B.l)    and  (B.3)    subject  to  the  respective  boundary  conditions. 

The  folding  technique  [26]  enables  us  to  simplify  this  problem. 
Provided  the  differential  equations  (B.l)  and  (B.3)  are  linear,  we  can 
add  them  to  get 

— 7+  — 2  +  ^  -  ^+  Q^=  0  (B.6) 

dx^       -b)!^       Zz-      \J-      ^n 


where  N  =  Nj^  +  N|_  i  s  the  folded  sun  of  the  electron  concentrations  in 
the  right-  and  in  the  left-side  grains  and  g  =  gn  +  g.  i  s  the  total 
generation  rate.  It  is  important  to  recognize  that  the  addition  of  the 
differential  equations  (B.l)  and  (B.3)  to  yield  (B.6)  is  possible  only 
if  Lp  is  the  same  in  the  two  grains.  The  boundary  conditions  are 
N(z=0)  =  0;   N(z  >  oo)   =  0;   N(x^)   =  0;   and  at  x=0, 

°n^M-^l      =  ^n(eff)t^R  ^  V       ' 


n     '^^       qGB 

°n  37-  ^n(eff)^       •  (B.7) 
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We    stress    that    in    order    to    obtain    (B.7)    we    have    assumed    that    quasi- 

equilibrium     prevails     across    the     grain-boundary     space-charge     region, 

PR 
which    ensures  S^^^^^^    has    the    same    value    on    either    side    of   the    grain 

boundary. 

We    have    thus    provided    a  mathematical    justification    for    using    the 

folding    technique.        The     solution    of   the    differential     equation    (B.6) 

yields   the    folded    sim    of  the   electron    concentrations    in    the   left-    and 

right-side     grains,      (N^     +     M^^),     and     is     equivalent     to     solving     the 

differential     equations    (B.l)    and    (B.3).      The    boundary   condition    (B.7) 

accounts    fbr    the    carriers    in    both    grains    that    recombine    at    the    grain 

boundary.      The    total    current,    I^biq,    given    by    (B.5)    can   be   calculated 

from  the   solution  of  (B.6)   as  described  in  Chapter  3. 


APPENDIX  C 

NUMERICAL  ALGORITHM  FOR  SOLVING  THE   TWO-DIMENSIONAL  CONTINUITY   EQUATION 

UNDERLYING  THE   EBIC   RESPONSE 

We  describe  here  the  basic  numerical  algorithn  used  to  solve  the 
two-dimensional  electron  continuity  equation  formulated  in  Chapter  3. 
It  should  be  pointed  out  that  our  algorithm  is  not  necessarily  the  most 
efficient  one.  The  method  was  found  to  be  successful,  and  hence  it  was 
utilized  without  additional   considerations. 

C-1     Mumerical    Solution  of  the  Steady-State  Continuity  Equation 

To    obtain     the    EBIC    response    we    must     first     solve    the    electron 

continuity   equation    (3,6)    under   low-injection   conditions   in   the    p-type 

base, 

a^N       3^N         N  ^     G       „ 
n 

subject  to  the  appropriate  boundary  conditions.  The  Gauss-Seidel 
elimination  method  [57,  58],  with  successive-overrelaxation,  is  used  to 
obtain  the   solution  of  (C.l),   N(x,z). 

A  rectangular  grid  system  was  employed  over  the  x-z  plane.  Two 
different  sets  of  grid  sizes  were  used.  The  finer  grid  system  was 
employed  in  the  region  around  the  generation  source.  Each  point  is 
labeled  (i,j),  which  indicates  its  coordinate  position.  I  f  ax  and  az 
are  the  coarse  mesh   sizes  in  the  x-  and  z-directions  respectively,   then 
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x^.  =   (i-l)Ax,       i   =   1,2,3,... 

(C.2) 
z.  =   (j-l)Az,       j  =   1,2,3,... 

are    the    coordinate     positions.        The     fine    mesh    sizes  ax'    and  az'    are 
calculated  as   follows: 


ax'    =   AX/(K-1) 


(C.3) 


Az/{L-1) 


where  K  and  L  are  the  ntnber  of  fine  mesh  points  in  the  x-  and 
z-directions  respectively.  The  values  of  K  and  L  are  chosen  by  trial- 
and-error  so  as  to  minimize  the  computer  time  and  maximize  the 
accuracy.     The  excess  electron  density  at  each  point  is  indicated  by 

N(x.,z.)   hN.^.        .  (C.4) 

A    finite   difference    form   of  (C.l)    can   be  written   at  each  interior 
point  in  the  coarse  mesh  region  using  the    following  approximations  [57]: 


8X  AX^ 

and 
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In  the  fine  mesh  reqion  similar  equations  can  be  written  by  replacing  ax 
and  Az  by  Ax'  and  Az'  respectively.  The  differential  equation  (C.l)  can 
now  be  written   in    finite-di  ff^rence    form   using   (C.5)   and  (C.6): 


«2  i,2  '     l2     'T 

n 


n 


(C.7) 


Rearranging,  we  obtain 


(2  *  2A  .  ax^/l2) '  =  •" 


"n 
2 


where  A  =  -^.      The    finite-difference    fom    in    the    fine-mesh  region   can 

also   be   rearranged   to   obtain   an   expression    for   N.   .    similar   to   that  in 

(C.8),   with  AX  and  Az  replaced  by  ax'    and  Az' . 

The    value    of  G^- ^    =   G(x^.,Zj)    in    (C.8)    is  evaluated    at   each   grid 

point  using  the  appropriate  model,  i.e.,  spherical  or  Gaussian,  to 
describe  the  generation  rate  distribution.  For  a  spherical -source 
model,  the  value  of  G  at  each  interior  point  is  given  by  a  constant 
value  (se4  Chapter  3),  and  for  a  Gaussian- source  model  G  is  evaluated 
using  (3.17).  The  total  generation  rate  G^  (electron- hole  pairs/ sec)  is 
chosen  such  that  the  theoretically  calculated  maximum  EBIC  in  the  grain 
is  equal    to  the  experimentally  observed  value. 

Equation   (C.8)    is  evaluated  at  each   interior  grid    point,    using   the 
dependent  and   independent  boundary  conditions  to  define   the  values  of  U 
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at     the    boundaries.        These    values    are    then    overrelaxed    by    using     a 
weighting    factor  w  [89]: 


N  =  wN         +    (l-w)N.    .  .  (C.9) 


Equation  (C.9)  is  used  in  conjunction  with  (C.8)  to  successively  assign 
new  values  to  the  excess  electron  density  at  each  interior  point.  The 
procedure  is  iterated  until  some  acceptable  difference  between  the  old 
and  new  values  of  electron  density  at  each  interior  point  is  achieved. 
The  nimber  of  iterations  required  to  obtain  a  given  accuracy  depends  on 
the  value  of  the  weighting  factor  used  as  well  as  the  coarse  and  fine 
mesh   sizes  used. 

The  sizes  of  the  coarse  grid  used  in  our  analysis  were  about  2  pm 
in  the  x-direction  and  3  ym  in  the  z-direction.  The  values  of  K  and  L 
were  49  and  13  respectively.  The  value  of  the  weighting  factor  used  was 
w  =  1.9  [89].  The  required  CPU  times  using  an  AMDAHL  470  computer  were 
typically  about  200  seconds  to  obtain  a  single  EBIC  line  scan,  i.e.,  to 
solve  (C.l)  and  to  calculate  the  EBIC  at  nine  different  positions  of  the 
beam. 

The  independent  boundary  conditions  at  z  =  0,  z  ^  «,  and  x  ^  =°  were 
incorporated  into  the  solution  by  simply  assigning  the  excess  electron 
density  values  to  the  respective  boundary  grid  points.  For  example,  at 
the  edge  of  the  junction   space-charge  region,   N^-    ^  =  0   for  i  =  1,2,3,... 

The  dependent  boundary  condition  at  x=0  is  modified  every  time  N  is 
evaluated.  The  dependent  boundary  condition  in  this  problem  involves 
the  normal  derivative  at  the  edge  of  the  grain-boundary  space-charge 
region,      and     can     be     numerically     evaluated     using     Newton's      forward 
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interpolation  formula  [58].  Differentiating  the  three- point  version  of 
this  formula,  we  obtain  the  following  one-sided  derivative 
approximation: 


9N^    .        -H         +   4N     i    -   3N 

-^--L2 ^ hi       for  j   =   1,2,3...        .  (C.iO) 


The   relation   (CIO)    can   be   used   to   evaluate,    after   each  iteration,   the 
electron  density  at  the  edge  of  the  grain-boundary   space-charge  region: 


-N     .  +  4N         -  3N 


PR 

In  (C.ll)  S^^g^^^  is,  in  general,  dependent  on  H  at  the  grain  boundary 
(see  Chapter  2),  and  is  evaluated  al^er  each  iteration.  Rewriting 
(C.ll) ,   we  obtain 


4N2  j   -  N3   . 


C.2     Numerical    Evaluation  of  the  Currents 
The  total   generation  current  in  the  base  is 

q  /T  G  dx   dz  =   qG^  =   iB^  .   iB^^  .   I^^^^        ,  (C.13) 


where  I^^   is  the  grain-boundary  recombination  current,    l^     is  the  quasi- 
neutral-base    recombination   current,   and    I^gj^  is  the  collected  junction 
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current.  The  values  of  the  current  on  the  right  side  of  (C.13)  are 
calculated  numerically  and  compared  with  the  total  generation  current 
qGq  to  serve  as  a  convergence  check.  The  numerical  procedure  is  stopped 
when  (C.13)   is  satisfied  to  within  5%  accuracy. 

We    now    discuss    the    nimerical     evaluation    of  ijr,,   I?.,,   and  Im-rr-- 

bo         [)N  ldIL 

The  grain-boundary  recombination  current  is  given  by 

PR 

Note  that  the  value  o^'S  ,  ^^v  at  any  point  along  the  edge  of  the  grain- 
boundary  space-charge  region  is  dependent  on  the  excitation  at  that 
point  (see  Chapter  2).  The  integral  in  (C.14)  is  evaluated  using  the 
trapezoidal  rule  [57].  The  qua  si -neutral -base  recombination  current  is 
obtained  by  integrating,  using  the  trapezoidal  rule  [57],.  the  electron 
density  divided  by  the  lifetime  in  the  base: 


^  0  0         n 


The  EBIC  collected  at  the  junction  is  given  by 


■eBIC'^^J       IiL^/x       •  '  =  •«' 


The  derivative  of  N  at  z=0  in  (C.16)  is  obtained  using  Newton's  forward 
interpolation  formula  [58]  in  the  z-direction.  Differentiating  the 
three-point  version  of  this  formula  yields  an  expression  similar  to 
(CIO): 
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aN.  -N         +  4N.        -   3N. 


The   integral    in  (C.16)    is  then  evaluated  using   the   trapezoidal    rule  [57] 
to  obtain   Ipgyr- 

The  boundaries  in  the  x-  and  z-directions  in  our  analysis  were 
fixed  at  about  three  diffusion  lengths  away  from  the  center  of  the 
generation  source.  The  electron  concentration  at  three  diffusion 
lengths  away  from  the  center  of  the  generation  source,  obtained  from  the 
numerical  analysis,  was  at  least  five  orders  of  magnitude  lower  than  the 
maximun  value,  which  justifies  our  choice  of  the  boundaries.  The  number 
of  coarse  and  fine  mesh  points  have  been  chosen  by  trial  and  error,  and 
the  method  yields  good  convergence.  The  FORTRAM  computer  program  to 
solve  the  continuity  equation  as  well  as  to  calculate  the  currents  is 
1  i  sted  below. 
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APPENDIX  D 
POLYSILICON  CELL  FABRICATION 


We  describe  below  the  processing  steps  involved  in  the  fabrication 
of  the  large-  and  gnall-area  poly  silicon  cells  used  in  bulk  and  thin- 
film  aluminum  diffusion  studies  (see  Chapter  4).  The  standard  cell 
cleaning  procedure,  which  was  common  to  all  the  processes,  is  described 
first.  The  photolithographic  process  (D.3)  was  used  in  the  fabrication 
of  the   small-area  cells. 

D.l     Standard  Cleaning  Process 

1.  Degrease   the  wafers  by    scrubbing   them  with  cotton   swabs  dipped 
in  hot  trichlorethylene   (TCE). 

2.  Immerse  the  wafers  in  boiling  TCE   for   five  minutes. 

3.  Immerse  in  a    second   solution  of  boiling  TCE    for    five  minutes. 

4.  Immerse  in  boiling  acetone    for   five  minutes. 

5.  Immerse  in  boiling  methanol    for    five  minutes. 

6.  Rinse  in  deionized  (DI)  water   for  ten  minutes. 

D.2     Bulk   Process  (Large-Area  Cells) 

1.  Clean   the  wafers  using   the   standard   process. 

2.  Polish  using  2:15:5  (HF:HN03:CH3C00H)  Si-etch  solution  for 
10-20  minutes. 
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3.  Clean  usinq  a  proxy  solution;  immerse  the  wafers  in  a  solution 
of  H2SO4  +  HgOgO:!)  for  about  10  minutes,  then  dip  in  HF  for  about  10 
seconds. 

4.  Predeposit  phosphorous  (n"*")   at  875°C    for  30  minutes. 

5.  Lap  the  back   surface. 

6.  Evaporate  1.5-ym-thick  aluminum  onto  the  back   surface. 

7.  Anneal    at  800°C   fbr  one   hour  in   the   presence  of  forming  gas. 

8.  Remove    the     phosphorous    glass     from    the     front    surface    using 
buffered-oxi de-etch  (BOE)    solution. 

9.  Evaporate  about  6000  A  of  aluminum  onto   the  back   surface. 

10.  Sinter,   if  necessary,   at  450°C    fbr  15  minutes  in    fonning  gas. 

11.  Etch  edge  using  6:1:1   (HN03:'rlF:CH3C00H)   Si-etch   solution. 

12.  Evaporate  aluminim  onto  the   front  surface  (~  300-A-thick) . 

D.3     Photo!  ithography 

1.  Bake  the  wafers  at  300°C   for  30  minutes. 

2.  Spin  on  negative  photoresist. 

3.  Prebake  the  wafers  at  65°C   fiDr  about  20  minutes. 

4.  Expose,   using  the  appropriate  mask,    for  about  15   seconds. 

5.  Develop    the    wafers    using    the    negative-resist    developer     fbr 
about  20  seconds. 

6.  Postbake  the  wafers  at  135°C   for  about  20  minutes. 

D.4     Bulk  Process  (Small-Area  Cells) 

1.  Clean   the  wafers  using   the   standard   process. 

2.  Polish  using  2:15:5  Si-etch   solution. 

3.  Clean   using   proxy   solution. 
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4.  Grow   a   thermal    oxide   about  5000-A-thick   at  1100°C:    10  minutes 
dry  Og,   22  minutes  wet  O2,   and  10  minutes  dry  O2. 

5.  Open    windows    for    the    n^    diffusion    using    the    photolithography 
process   (1st  mask) . 

5.       Etch  off  the  oxide  in  the  windows  using  BOE   solution. 

7.  Remove   photoresist  using   proxy   solution   (~  8  min). 

8.  Clean   using   proxy   solution. 

9.  Predeposit  phosphorous  at  950°C    for  40  minutes. 

10.  Clean   using   proxy   solution. 

11.  Drive   the   emitter   in   at  1000°C    for  27  minutes:    10  minutes  dry 
O2,   12  minutes  wet  O2,   5  minutes  dry  O2. 

12.  Lap  the  back   surface. 

13.  Evaporate  about  1.5  ym-thick  aluminum  onto   the  back   surface. 

14.  Anneal    at  800°C    fbr  one   hour  in    forming  gas. 

15.  Open  windows   for  the    front  contacts  using  the   photolithography 
process  (2nd  mask) . 

16.  Etch  off  the  oxide  in  the  windows  using  BOE   solution. 

17.  Remove   photoresist  using   proxy   solution   (~  8  min.) 

18.  Clean   using   proxy   solution. 

19.  Evaporate  aluminum  onto   the    front  surface  (~  5000-A-thick) . 

20.  Define        the        front       metallization        pattern        using        the 
photolithography  process. 

21.  Remove  aluminum  using  aluminum  etch. 

22.  Remove   photoresist  using   plasma  etch. 

23.  Evaporate  alumintm  onto   the  back   surface   (~  6000-A-thick) . 

24.  Sinter,   if  necessary,   at  450°C    fbr  15  minutes  in    forming  gas. 
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D.5     Thin-Fi1m  Process  (Large-Area  Cells) 

1.  Clean   the  wafers  using   the   standard   process. 

2.  Polish  using  2:15:5  Si-etch   solution. 

3.  Clean   using   proxy   solution. 

4.  Evaporate   1.5  urn  of  aluminun  onto  the    front   surface. 

5.  Anneal    at  450°C    for  3-6   hours  in    forming  gas. 

6.  Remove  aluminun  using  aluminim  etch. 

7.  Etch    the    front    surface    (0-2   pm)    using    1:2:9    (HF:CH3C00H:HM03) 
Si -etch  solution. 

8.  Clean   using   proxy   solution. 

9.  Predeposit   phosphorous  at  875°C    for  30  minutes. 

10.  Lap  the  back   surface. 

11.  Evaporate  aluminum  onto   the  back   surface  (~  6000-A-thick) . 

12.  Edge  etch  using  6:1:1   Si-etch   solution. 

13.  Evaporate  aluminum  onto   the    front   surface   (~  300A-thick) . 
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